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O ruchu planetoid typu *, (Thule). 


Wyznaczenie wyrazów typowych oraz przerwy. 


Wstęp. 


Wiadomo, że za charakterystykę danej grupy planetoid służy 
wymiernosé stosunku ruchu dziennego Jowisza, Marsa lub Saturna 
do ruchu dziennego planetoidy. Gdy dana planetoida posiada ruch 
dzienny taki, iż stosunek ruchu dziennego planety głównej do ru- 
chu dziennego planetoidy wyraża się w przybliżeniu stosunkiem 
dwóch małych liczb całkowitych, jak np. 4, 4, 2, 2, ... etc., to zabu- 
rzenia w ruchu tej planetoidy, spowodowane planetą główną, są 
tem większe, im bardziej się ów stosunek ruchów dziennych zbliża 
do stosunku takich dwóch liczb. Największe zakłócenia w ruchach 
planetoid powoduje Jowisz, dzięki swej olbrzymiej masie; niektóre 

tylko planetoidy, których ruch dzienny jest współmierny z ruchem 
Marsa, podlegają dość znacznym zaburzeniom z jego strony, np. 
Hebe (6), Vala (131), Asia (67), Nysa (44), Lutetia (21) o dość 
znacznych ruchach dziennych i mimośrodach. Zaburzenia Jowi- 
szowe wywierane na tę grupę planetoid nie są duże, a toz powodu 
znacznej odległości Jowisza od tej grupy, Dopiero planetoidy gru- 
py 3 (z ruchem dziennym około 900"), o małym mimośrodzie, pla- 
netoidy grupy i (z ruchem dziennym około 600") właściwie rozpo- 
czynają grupę planetoid silnie podlegających zaburzeniom ze stro- 
ny Jowisza. O ile stosunek ruchu planetoidy i planety głównej nie 
jest bardzo blizki stosunku dwóch małych liczb, to perturbacye 


O ruchu planetoid ] 


EZ i 


2 J. KRASSOWSKI 


nie będą zbyt wielkie, i wtedy zwykłe metody obliczania perturba- 
cyj stosują się bez trudności. 

Gdy jednak ów stosunek staje się coraz biiższy takiego pro- 
stego ułamka, jak np. dla przedstawicieli każdej grupy planetoid, 
a więc Hestii, Hekuby, Hildy etc., których ruchy dzienne najbar- 
dziej są zbliżone do pewnej niewielkiej wielokrotności ruchu Jo- 
wisza—trudności przy obliczaniu perturbacyj stają się coraz więk- 
sze; zwykłe metody zawodzą, gdyż w wyrażeniach na perturbacye 
występują bardzo małe dzielniki, które przy całkowaniu jeszcze 
bardziej się zmniejszają. 

Przy pierwszem przybliżonem obliczeniu perturbacyj takich 
planetoid należałoby już uwzględnić wyrazy rzędu kwadratu masy 
zakłócającej, co, jak wiemy, niepomiernie komplikuje rachunki 
i często nawet uniemożliwia ich wykonanie o ile chcemy postępo- 
wać według metod klasycznych. 

Gylden podał właściwie jednę z metod przybliżonych, która, 
jak paru uczniów jego wykazało (Brendel, Backlund, Harzer 
i їп.) те zawodzi nawet i wtedy, gdy stosunek ruchów jest bardzo 
blizki do prostego ułamka. Metodę Gyldén'a, dość skompliko- 
waną i niezbyt dobrze przystosowaną do rachunków praktycznych, 
przerobił i uzupełnił Brendel. Metodę tę uczniowie Gyldén'a 
zastosowali już do planetoid grupy Hestii, Hekuby i Hiidy — i za- 
wsze z dobrym rezultatem. Brendel i Kramer podali obszer- 
ne tablice, pozwalające w niedługim czasie obliczyć najważniejsze 
perturbacye planetoid typu Hestii i Hekuby i ich pozycye geocen- 
tryczne, tak że prawdopodobnie z łatwością będzie można zawsze 
odnaleźć te planetoidy. Wyniki dla ruchu Hekuby, otrzymane 
przez Harzer’a!), zostały do pewnego stopnia skontrolowane na 
innej drodze przez Simonin'a?). Przypadek planetoid typu Hil- 


dy (5) trudny z powodu bardzo blizkiego sąsiedztwa Jowisza z te- 


mi planetoidami, zbadał Buchholz (obliczył on charakterystycz- 
ne wyrazy, ale nie podał jeszcze wyrażeń samych perturbacyj). 
Planetoida Thule (279) należy do oddzielnej klasy planetoid: 


= i jest, jak dotąd, jedynym tej klasy przedstawicielem. Thule je 


1) Harzer P. Untersuchungen über einen speziellen Fall... etc. 


2) Simonin M. Sur l'orbite d'Hécube, Annales de l'observatoire de 
Nice, t. IV. 
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szcze bardziej niż planetoidy typu Hildy może się zbliżyć do Jowi- 
sza; takich zbliżeń, jak wspomina Bidschof!), bywa po trzy na sto 
lat. Podczas takiego zbliżenia z Jowiszem, ponieważ Thule przez 
dłuższy czas znajduje się pod jego siinym wpływem zakłócajacym, 
droga jej ulega bardzo znacznym zmianom; np. dzięki zbliżeniu 
do Jowisza, które miało miejsce w 1912 r., średni ruch dzienny 
Thule znacznie się zmniejszył, jak to widzimy z rachunków Vilje- 
va?) Kobb?) przypuszczał nawet na zasadzie wzorów S.H. 
Darwin'a*), odnoszących sie do orbit okresowych, że ruch Thule 
może z czasem stać się niestałym i droga jej może się nawet kiedyś 
przeciąć z drogą Jowisza. 


Dotąd nie badano ogólnie ruchu tej ciekawej planetoidy. 
O ile wiem, tylko Bidschofi Wedemeyer?) obliczyli lub po- 
prawili pierwiastki drogi Thule, a Viljev obliczył tylko zaburze- 
nia ruchu, odnoszące się do niedługiego okresu czasu, i to opie- 
rając się na danych orbity, otrzymanych przez Wedemeyer'a. 

Wydało mi się więc rzeczą inreresującą przystąpić na zasa- 
dzie metod, podanych przez H. Gylden'a, które z tak wielkiem 
powodzeniem dopełnił M. Brendel, do zbadania ruchu tej cie- 
kawej grupy planetoid typu Thule, których ruch średni dzienny 
iest bardzo blizki 4/, ruchu średniego dziennego Jowisza. Punk- 
tem wyjścia naszych rozważań będą dane liczbowe, odnoszące się 
do planetoidy Thule, podane'przez Wedemeyer'a. 

W niniejszej pracy zamierzam streścić ogólny zarys teoryi 
Gyldén'a, obliczenia perturbacyj małych planet w modyfikacyi, 
wprowadzonej przez Brendel'a, i następnie wyznaczyć charakte- 
rystyczne wyrazy, odnoszące się do ruchu planetoidy Thule, oraz 
szerokość przerwy *) w pierścieniu planetoidalnym, która dla typu 
Thule dotąd nie została jeszcze oznaczoną. 


1) Sitzungsberichte d. Akad. d. Wissenschaften in Wien. Bd. 100 
Matlı.-Naturwiss., p. 937. J. Bidschof, Bestimmung der Bahn des Planeten 
(279) Thule. 

7) Viljew AN Bd. 195 

3) Kobb G. Sur un cas d'instabilité possible, Bull. Astr. 1902. 

1) Darwin G. Н. On periodic orbits, Acta Mathematica t. XV. 

5 Wedemeyer. Archiv der Seewarte XXXI Jhrg. 1909. 

5$) Brendel. Theorie d. kl. PI. II, p. 5. 
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I. Wyprowadzenie równań ruchu zakłóconego w spół- 
rzędnych Gyldenowskich. 


1. Wiadomo, że ogólne równania ruchu trzech ciał, które 
charakteryzują ruch zakłócony jakiejś planety dokoła słońca, sa 
następujące: 


dex N ТМ, 20 
ap TA (1 +m) == Um) e 
L dy 2 o 29 
dt? Kë k (1-- m) y3 af (m 4-1) dY 


d’z 5 cui 90 
dt: al) r3 = (m1) ES 


gdzie, według przyjętego zwyczaju, przez k oznaczamy stałą Gaus- 
sa, przez m masę planety zakłóconej (przyjmując masę słońca za 
jedność), r promień wodzący tej planety; m, r” takie same wiel- 
kości dla planety zakłócającej; x, y, 2 і x', y”, 2 odpowiednio dla 
każdej są to spółrzędne prostokątne, odniesione do pewnego stałe- 
go kierunku w przestrzeni i których początek jest w środku słońca; 
t czas liczony w dniach średnich słonecznych, który upłynął od ja- 
kiejś epoki początkowej. Funkcya © jest to funkcya perturbacyjna: 
Papy Ada da yy sette e mi Й. r 
na ]4-m 5 ri que |= l+ m ^ „ORŁA | 
gdzie 
A? — (x — a+ (y — y} -- (e — z)? 2 r?-- 7r? —2rr' . cos Н 
zaś H jest to kąt w środku słońca, pomiędzy promieniami wodzą- 
cymi r i r' obydwóch planet. 

Równania | należy przedewszystkiem odpowiednio prze- 
kształcić. 

Н. Gyldén analogicznie do myśli Ha nse n’a rozpatruje ruch 
planetoidy zakłóconej na swej chwilowej drodze oddzielnie od ru- 
chu samej płaszczyzny drogi w przestrzeni. 

Hansen okazał), że istotnie można znaleźć takie wyrażenie 
analityczne dla ruchu zakłóconego planetoidy, iż elementy chara- 
kteryzujące położenie orbity w przestrzeni, a wiec 8:12 ИШ 
związane stałe, wystąpią oddzielnie od pozostałych elementów. 
Taki rozdział tych dwóch rodzajów charakterystycznych wielkości 


Hansen P. A. Auseinandersetzung etc. I. $ 8. 
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można osiągnąć przez wprowadzenie pewnego nowego układu 
spółrzędnych z,, Y,, 2, mającego — jak zresztą i poprzedni układ 
czyż — swój początek w środku ciężkości słońca. 

Całkowanie układu równań różniczkowych 1, w których pra- 
we strony byłyby zerami, dawałoby nam, jak wiemy, ruch po elip- 
sie Keplerowskiej. 

Równania te są następujące: 
dex 
dt? ab 
d*y 
dt? 
(^ PN NS 

dt? rr ya k2(1 + m) —0. 

Stałych całkowania wtym wypadku jest 6; są one dwóch 
rodzajów: jedne określają położenie drogi planetoidy w przestrze- 
ni, pozostałe zaś kształt i rozmiary drogi. Metoda waryacyi stałych 
dowolnych, jak wiemy, upoważnia nas, zamiast całek równania I 
(ruchu niekeplerowskiego, zakłóconego) wziąć całki równania (2) 
dla ruchu Keplerowskiego— pod warunkiem atoli, abyśmy stałe 
dowolne, które są w całkach równań ruchu keplerowskiego znajdu- 
ją, uważali jako funkcye czasu, które muszą spełniać równania |. 
dz du dz 
dt” dt ' dt 
obliczone na zasadzie układów równ. 2 1 I mieć będą i te sama po- 
stać analityczną. Wszelkie spółrzędne, spełniające powyższy waru- 
nek (t. j. mające taki sam kształt analityczny w ruchu Keplerow- 
skim i w ruchu zakłóconym) Hansen nazwał spółrzędnemi ideal- 
nemi ?). 

Niech te spółrzędne idealne będą z,, y,, 21; dostawy kierun- 
kowe osi x, oznaczymy przez о, B. y, OSi y, przez a,, Bi, 11 a OSI 2, 
przez à, Bo, Ya — wówczas, jak wiemy; zachodzić będą związki 
następujące: 


(8) zy — ax -by + 12; yi eut PY HTN 2, = art Fay 122 
1 odwrotnie 
(4) ж == ad, 1-04, + 22, y — Br. + 810, + 21; 

g = т tfi + 281 


x 
73 k?(1 + in) = 0 


(2) LI ка 4+m=0 


y3 


Wtedy wartości spółrzędnych x, y, 2 i ich pochodnych 


! Hansen P. A. Auseinandersetzung etc. l. $ 8. 
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| % + BB, + tt = 0 SP RUN. A M GÓR 
aas -- Bg + її; = 0 oi Bi r 11 = 1 
tt Bio "fa = 0 aż 4- 83 -|- 15 = 1. 
Aby nowe spółrzędne x,y,2, były „idealne“, wystarcza aby 
były spełnione związki następujące: 


(4a) 


do | dB É d'i LT ПАС ES 
Z dt TY dt yc xy mU mar mese sr Y 
do, , | dä | dto __ 
VETAT аа = 


Związki (5) redukują się do dwóch tylko warunków niezalez- 
nych; jakoż uwzględniając związki (3) i podstawiając w równ. (0): 


(9) APCE 


gdzie dla krótkości oznaczyliśmy : 


Adt = N ap=— Wadi 
Bdt = У ady=— Y yda 


Сй = > dao — Ż «dr 
(wielkości te są to dostawy kierunkowe chwilowej osi obrotu). 

Z równania (6) widzimy przeto, że we wszelkim układzie 
spółrzędnych idealnych, odniesionym do osi ruchomych, oś chwi- 
lowa obrotu układu schodzi się z promieniem wodzacym planetoi- 
dy. Jeśli założymy, że stała 2,=0, wtedy płaszczyzna z, y stale 
zawierać będzie w sobie promień wodzący planetoidy i nowy układ 
spółrzędnych ruchomych x,, y,, 2, będzie najzupełniej określonym 
w przestrzeni. 

Ponieważ wprowadziliśmy warunek !) 2, =0, mieć więc bę- 
dziemy 


(7) AT Bay -1- 122 — 0; 


1) Należy przypomnieć, że Hoene Wroński, podając swoje nowe metody 
obliczenia perturbacyj w „Reforme absolue du savoir humain*, rozpatruje per- 
turbacye, odnoszące się do płaszczyzny drogi planetoidy oddzielnie od pertur- 
bacyj elementów jej orbity. 
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to równanie, jak to widzimy od razu, jest równaniem płaszczyzny 
przechodzącej przez promień wodzący, a zmiennej z czasem ż; 
płaszczyzna, zawierająca w sobie w czasie ї nietylko promień wo- 
dzący planetoidy dla tej chwili, ale jeszcze i styczną do jej drogi, 
spełniać będzie równanie, w którym zamiast z, y, 2 podstawimy 
xz--dz, y + dy, z+- dz. 

Wówczas muszą być spełnione związki następujące: 


(8) x dz -|- Pady + 1402 = 0 
Oraz 
(9) ada, + YAB, + 201, = 0. 


Równania (7) i (8) dają nam warunek, który musi być speł- 
niony, aby płaszczyzna =, у; znajdowała się w płaszczyznie chwi- 
lowej drogi planetoidy; przy pomocy zaś równania : 

(10) ode, + 8d8, + ydy, = 0 
wyznaczymy na tej płaszczyznie położenie osi Oz. 

Pomnóżmy równania (4) raz przez da, dß, dy, a drugi raz 
przez da, dB, dr, dodajmy i uwzględnijmy związki: (4a), 2, =0, 
oraz (10); otrzymamy: 

ado. +ydf +12d = 0 

xda, + ydg, + 207, = 0. 
Równania (3) zróżniczkujmy i uwzględnijmy równ. (9) i (11), 
wtedy otrzymamy, dołączając równanie 2,=0, równania nastę- 
pujące: 


(11) 


| d= ade Hay + pz 
(12) | dy, = a dz + Bd y + 1 d2 
dz, = 0 = adx + 8d y + 12. 


Równania (12) pomnózmy przez а, Gu, œs, później drugi raz 
przez B. B,, Bs, i wreszcie trzeci raz przez Y, 11, Ta 1 dodajmy, to po 
uwzględnieniu związków (4a) otrzymamy : 

dz = ada, + ady 
dy = kde, + BP dy, 
de = tax, -+ Td, 


Z tych równań po uwzględnieniu (4a) otrzymamy: 
dada + азау + dydz = 0 
da dx + dB,dy + dr dz = 0 
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po uwzględnieniu których z (12) otrzymamy równania następujące: 
dix, =ad'x + Bdiu Lade, 

(13) d'y, = o d?z LB, d'y + dz, 

| d?2, 5050122: 
du = adx, + o,d*y, + dadx, + da, da, 
d'y = 8022, -- 8, dx, + 4802, + df, dy, 
de = Ма, + dy + d'tyz, + di dy, . 
Funkcya ice Q da nam dł następujące: 


H ` «28 
Qd, a $3 Aë 


08 a, 1, d 


14 2 af = EM " 
( ) Y Az a | B, dY | F az 
08 0 
B rese 


Jeżeli wreszcie równania I odpowiednio pomnożymy przez 
a, B, 1 i dodamy do siebie po uwzględnieniu dopiero co otrzyma- 
nych związków, to mieć będziemy !: 


dex, 


dg + k* (1 +m) 5 "E + m) = 
(1) 


dk +k (1 +m) =m (14m) 5 


Układ równań II jest najzupełniej m ei równania | 
i wyraża ruch planetoidy po chwilowej drodze”). W układzie rów- 
nań II zamiast spółrzędnych prostokątnych wprowadźmy spółrzę- 
dne biegunowe r iv. Ponieważ płaszczyzna drogi planetoidy leży 
w płaszczyznie x,y,, więc wskutek perturbacyj droga ta bynaj- 
mniej nie jest zamkniętą, i dla tego musimy liczyć kąty v nie tak, 
jak zwykle w elipsach od 0% do 360°, lecz od -- oo do -|- oo 

Spółrzędne biegunowe, które wprowadzamy do równań II, są 
nam dane przez równanie następujące: 


(15) L = MEOS V; yQm—r-sin ©; 


zaś rów. 11 przejdą w następujące: 


! Hansen. Auseinandersetzung l. p. 69. 
2) Harzer P. Untersuchungen etc. p. 8. 
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d?r dot Kl m). 2Q 
(HD) | dt? r (i) tech, 3 FRUT ee 
dv dr, dv 1 08 
d | =—> 2 = 
| dt dt ' | dt? A y du 


W tych równaniach v oznacza prawdziwą długość planety, 
zaś 7 jest to jej promień wodzący. Te równania są to równania 
zasadnicze Haesen'a, ale w spółrzędnych biegunowych. 

2. Równ. ШН. Gyldćn przekształca przez wprowadzenie pe- 
wnych nowych zmiennych zależnych p, n, S oraz zmiennej niezależ- 
nej v „prawdziwej długości*; te wielkości są to tak zwane spółrzęd- 
ne Gyldénowskie. Zanim jednak to przekształcenie równania 
III wykonamy, przypomnimy w ogólnych zarysach bieg myśli Gyl- 


30 3 


dćna. Jeżeli w równaniach III przypuścimy, że > I S= sa równe 
0, to wtedy te równania przedstawiać będą ruch Keplerowski 
i, jak to wiemy, otrzymamy całki następujące: 
bre e e 
ПУ = 1 Fe cos (от) 
oraz stałą prędkość wycinkową 
„du — — 
y AH Va (1— e?). 
Ponieważ masę planetoidy zawsze możemy zaniedbać wobec 
masy słońca, którą przyjmujemy za 1, więc w ruchu Keplerow- 
skim jest zawsze: 


(2) 


(3) a(l—e) <r < a(l --e). 


Jeżeli >" xd 2,0, to mieć będziemy do czynienia z ruchem 
zakłóconym, i wtedy prawdziwa droga planetoidy znacznie się róż- 
nić będzie od elipsy Keplerowskiej, jednak, jak to nam wyka- 
zują spostrzeżenia, dla czasu skończonego promienie wodzące pla- 
netoidy wahać się będą w pewnych skończonych granicach. Otóż 
Gyldén wprowadził tak zwane krzywe periplegmatyczne 1), zawar- 
te między dwiema współśrodkowemi kulami, skończenie od siebie 
oddalonemi. Krzywe te są to poszczególne drogi planety, które się 
tworzą wskutek ciągłych zmian wielkości połowy osi jej orbity. 


') Gyldén H. Orbites absolues г. 1. 1. 1; id. Undersökningar... 1. 
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Zmiany te, według Gyldćn'a są nieograniczone w czasie, jednak 
zawarte są w skończonych granicach; wielkość a Gyldén nazywa 
protometrem. Ponieważ w dalszych rozważaniach naszych 
a nie będzie absolutnym *) elementem, przeto, idąc za Brendel'em, 
wielkość a nazywać będziemy wprost połową osi drogi (nie 
jest to jednak półoś chwilowej elipsy Keplerowskiej). 

Przypuśćmy, że planeta zakłócona porusza się po drodze, le- 
żącej w pewnej stałej płaszczyznie; niech 7 będzie jej promień wo- 
dzący, v—prawdziwa długość, liczona od pewnego stałego kierunku. 
Wskutek perturbacyj elipsa Keplerowska stale się odkształca, 
irównanie drogi, opisanej przez planetę zakłóconą, napisać moż- 
na, według Gyldén'a, w postaci następującej: 


(4) ping 10) 


G ТЕЗ] 4 
gdzie a jest to pewna średnia wartość promienia wodzącego, a 1p 
są to funkcye prawdziwej długości, które przedstawić możemy 
przy pomocy szeregów trygonometrycznych; n jest to pewna wiel- 
kość, bardzo powoli zmienna z czasem i porównywalna z mimosro- 
dem eliptycznym. W ruchu Keplerowsk iin prędkość wycinkowa 
była stałą, w ruchu zakłóconym, według Gyldé n’a, już tonie ma 
miejsca, gdyż ta prędkość wycinkowa waha się dokoła pewnej Sre- 
dniej wartości. Przez analogię z ruchem Keplerowskim, mo- 
żemy więc napisać, że w ruchu zakłóconym zachodzi dla prędko- 
Sci wycinkowej równanie nastepujace: 
dv Foi — 13) 
Va u E RE 

(5) es 
gdzie S jest to pewna mała wielkość. Wprowadźmy funkcye S, p, 7 
w równania III. 

Drugie z równań Ш możemy napisać: 
d > dv — 29 
dt (7 LL AE ET 

Połóżmy m= 0; wtedy będziemy mogli napisać przez analo- 
gię, ze wzorem na ruch Keplerowski, równanie następujące: 


1) Absolutną droga Gyldén nazywa pewną orbitę, w której elementy są 
to pe wne zupelne stałe, i która różni się od rzeczywiście opisywanej przez da- 
ne ciało drogi, o wiełkości rzędu siły zakłócającej. 
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d JiVWa(i — 1?) |. ,,99. 
dt | 1--S Ј дф 


wykonajmy wskazane różniczkowanie, wtedy otrzymamy: 


I kVad—qp)| Mali) 48 1а 1 dP) H 
pases) d do x 2 1358; ac dv] 


A ` 
` 


czyli, jeśli wraz z Gyldén'em i Brendel’em') oznaczymy: 


y 20 


[V bis br pues NEA 


powyższy wzór przejdzie w równanie następujące, które nam okre- 
sli S: 
1 l T ` 


eh re EE 


Zwróćmy się teraz do wzorów (4) 1 (5); bez trudności znaj- 
dziemy: 


DA 
dr ` Kol “бу 
dt 1-8 dv 
oraz 
d? Í sl 2) d d 
AT dd Ad d TO IE 
ZONE od re ap ta S on el i 


Uważajmy teraz równania Ill; w tych równaniach podstawmy 
2 9 


à ' der ¡dry | 
znalezioną dopiero co wartość na oraz na Eé , to otrzymamy: 


dt? 


| d 


| | „29 
a) a(l — 1?) | шу + 0-870 | -ü += —n 0-087 


а(1—1 


2 
Równanie r=--.- ) da nam: 
ken 


) Brendel. Om Användningen p. 5; Gyldćn, Undersókningar II. 
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r | EM а ар 


со podstawione w równanie (а) da nam w rezultacie: 


5 Г? 
(09 $$ p > WE 14874 |g, 2848 — 
| ! 2 eg 
(ISP – _ es ha y Gen 4 
(1-8); 
ER | — v? pal IP) 


gdzie analogicznie do poprzedniego przez P oznaczylismy funkcye 
następującą: 

(Vets) Р — ү? 
Równanie (V) daje nam właśnie wyrażenie na p w funkcyi v. 


. Ponieważ P i Q są to pewne szeregi trygonometryczne, więc 
równania IV i V mieć będą postać następującą: 


CH) ES = N z, Sin (Anv — B,) 
de 4 | 
(2) Tp (1 — В)р = y b, COS (A„v — Bn). 


Zgodnie z oznaczeniami Brendel’a przyjmujemy, że w równ. 
(1) i (2) wielkości а,, On, № Bn, B sa to pewne stałe wielkości; 
w wielkościach а„, Dw, B, (co jest istotne) masa ciała, wywołujące- 
go zaburzenia (w tym przypadku masa Jowisza), występuje jako 
czynnik. Te trzy wieikości są przeto porównywalne z masą zakłó- 
cającego ciała. 

Całkowanie równania (1) da nam odrazu rezultat następujący: 


(3) S=a — Y cos (A,v В), 


gdzie przez a, oznaczylišmy stałą całkowania. 
Z równania zaś (2) otrzymamy: 


D. ¿Bremdici Dis do skl sPlanetengl: 


ag + nein © PL де ET 
dv? al y? lav: 1—12 do dv!" (1—p9 dv ln dw 
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? EU COS (A, — Ba). 

Wyrażenia (3) i (4) zawierają jako dzielnik wielkość zwykle 
małą A, rzędu masy zakłócającej; wskutek całkowania we wszyst- 
kich wyrazach zależnych od A, wystąpią wielkości, w których masa 
zakłócająca zachodzi w mianowniku w wyższej potędze, wskutek 
czego wyrazy zależne od A, zostaną bardzo powiększone. 

W równ. (1) istnieją wyrazy pomnożone przez masę ciała za- 
kłócającego w pierwszej potędze; te wyrazy w całce (3) już nie bę- 
dą zawierać tej masy jako czynnika. Podobne wyrazy Gylden 
nazywa elementarnymi; ponieważ w przypadku tych wyrazów 
elementarnych spółczynniki A, będą bardzo małe (masa Jowisza 
dix) przeto te wyrazy będą bardzo powolnie zmienne z czasem 

| i według określenia Harzer'a!) są to wyrazy elementarne 
| długookresowe. 
| Jeżeli ruch planetoidy i planety zakłócającej jest tego rodza- 


(4) pol cos [V1 —fu—I'] + > 


ju, że stosunek ich ruchów średnich dziennych jest prawie równy 


p 


prostemu ulamkowi Cé to wtedy też w równ. (3) wystąpią wyrazy 


będące pierwszego rzędu w stosunku do masy planety zakłócają- 
cej; w niektórych przypadkach specyalnych mogą nawet być rzędu 
niższego od pierwszego, a mianowicie takiego, jak pierwiastek 
trzeciego stopnia z masy zakłócającej. Wyrazy, których rząd jest 
zawarty pomiędzy 1 i 1—Gyldén nazywa charakterystyczny- 
mi; ponieważ jednak w tych wyrazach czynnik A, jest bardzo ma- 
ły, przeto i te wyrazy będą charakterystyczne długookre- 
| sowe. Wreszcie w równ. (1) mogą występować takie wyrazy, dla 
| których czynniki A, nie będą małemi wielkościami, wskutek czego 
po zcałkowaniu w (3) wyrazy te nie będą zwiększone; podobne 

| wyrazy Są to wyrazy zwykłe. 
Oznaczmy przez 5, spółczynniki przy v w wyrażeniach wyra- 
zów elementarnych długookresowych (o" jest rzędu masy zakłóca- 
| jacej); przez 9, tę samą wielkość w wyrazach charakterystycznych 
| długookresowych (9, jest wielkością małą, jednak znacznie więk- 
szą, niż masa zakłócająca); wtedy będziemy mogli napisać równ. 

(3) w postaci następującej: *) 


1) Harzer P. Untersuchugen, р. 3. 
2) Brendel. Theor. d. kl. Planoten I. p. 18 et sqq. 


1 4 J. KRASSOWSKI. 
vwd xa! 4 
(5) S= ¿a — m e cos (5,4 — An) > "COS (040 — Gu 
n On 


af ` ` 
ÑU cos (Av — Ba) 


p An 


i w tem wyrażeniu A, już nie jest małą wielkością. 

Zbadajmy obecnie w analogiczny sposób równ. (2). Przede- 
wszystkiem należy tutaj zauważyć, ze 1 — B == A, wskutek czego 
w (4) nigdy nie wystąpią dzielniki zerowe. Z pośród różnych wy- 
razów 2 А, wchodzących w równ. (4), musimy odróżnić przede- 
wszystkiem wyrazy, w których A, różni się od jedności tylko o wiel- 
kości rzędu masy ciała zakłócającego; w tych wyrazach A, będzie 
mieć kształt 1 —5,, gdzie 5n, jak i poprzednio, oznacza wielkość 
rzędu masy zakłócającej. W równ. (4) w odnośnych wyrazach wy- 
stąpią dzielniki tego samego rzędu, wskutek czego te wyrazy będą 
wyrazami elementarnymi; jednak będą one elementarnymi 
krótkookresowymi. Dalej w równ. (2) wystąpią także wyrazy, 
w których A, różni się od 1 o wielkości rzędu б,, tak iż możemy 
napisać A, = 1 — 2, w tych wyrazach występują dzielniki rzędu 0,; 
wyrazy te są to wyrazy charakterystyczne krótkookre- 
sowe. 

Wreszcie w równ. (2) Są i takie A,, które bardzo się różnią od 
jedności; wyrazy, zawierające te wielkości, przy całkowaniu nie 
uzyskują małych dzielników, nie są więc one powiększone: są to 
wyrazy zwykłe. 

Podobnie, jak to zrobiliśmy dla S, możemy napisać równanie 
na p w ten sposób, że każdy rodzaj wyrazów charakterystycznych 
występować będzie oddzielnie; p zatem mieć będzie postać nastę- 


pującą: 
| p= k cos [(1— <) v — T] + ХК, cos (сл) — Bal 
| +B, cos [Il —à,) v—D]+ Zb, cos [Anv — Bn) 


gdzie, idąc za Brendel'em, oznaczyliśmy : 


(6) 


hb b 
7 1 ae” ү T E n : M n ^ 
(7) ° an pda Pag] i5 Bu UHR 
B Poe On Б") Om , Я рп RZA: Dn” 
VM 20, - 0,? V. p T 2(0, FER c) UE (2,2 c?) zT l ед Anz — B | 


D Harzer. Untersuchungen, p. A 


| 
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Wyrazy charakterystyczne dtugookresowe i krótkookresowe 
wyrazy elementarne obydwu kategoryj odgrywają zasadniczą rolę 
w rozwinięciach wielkości /’i Q, które wchodzą w wyrażenia samej 
funkcyi perturbacyjnej. 

Harzer !) przejrzyście oznaczył w następujący sposób te 
wszystkie wyrazy, wypisując dla każdego rodzaju charakterystycz- 
ny kąt, którego sinus lub cosinus zachodzi w naszych wzorach. 
Harzer rozróżnia cztery postacie A, B, C, D wyrazów, charaktery- 
zujące się następującymi argumentami: 


(8) | A) cw— An; C) Bai Ca; 
EBY (E DA 


Wyrazy A)i B) są to wyrazy elementarne; C) i D) — chara- 
kterystyczne. Wyrazy kształtu A) i C) są długookresowe; zaś B) 
i D) krótkookresowe; wszelkie wyrazy, których argumenty nie są 
kształtu (8), są wyrazami zwykłymi. Należy zauważyć, że najniż- 
szy rząd względem inimośrodów drogi ciała zakłóconego i zakłóca- 
jącego dla wyrazów kształtu (A) jest rząd drugi, zaś dla wyra- 
zów (B) już pierwszy. Najniższy rząd wyrazów kształtu (C) i (D) 
względem mimośrodów, zależy od stosunku ruchów średnich y. pla- 
nety zakłóconej i planety zakłócającej. Niech stosunek ten mało 
się różni od stosunku dwóch liczb całkowitych pi g, wtedy ów 
rząd najniższy będzie | p ~ q | albo į | (p—q)—1 |l; w przypad- 


ku planetoid typu Thule, dla których 5 mało się różni od %/,, wy- 


razy (B), (C) i (D) będą już pierwszego rzędu, a więc bardzo znaczne. 
Przy całkowaniu, które dało początek spółczynnikom b,, mo- 


że się zdarzyć, że niektóre z tych b„ wypadną niezupełnie elemen- 
tarne, t. j. takie, które nie są rzędu 0 względem m’, lecz wyższego. 
Te wyrazy Gyldćn nazywa wyrazami podelementarnymi. 
Wyrazy elementarne zawierają w sobie perturbacye wiekowe, w któ- 
rych według klasycznych teoryj czas występuje przed znakami 
funkcyj Sinus i cosinus. Gyldćn przez wprowadzenie wyrazów 
elementarnych uniknął tej niedogodności. 

Wielkością stopnia n-go nazywać będziemy wyrażenie, które 
zawiera jako czynniki niektóre z wielkości y, 4 kn, k, przyczem 
suma wykładników przy tych czynnikach jest równa n. Brendel 


1 Harzer. Untersuchungen, p. 3. 
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do teoryi małych planet wprowadza jeszcze określenie następujące. 
Wiadomo, że wielkością rzędu л-ро nazywamy wielkość, która za- 
wiera jako czynnik masę zakłócającą w potędze n-tej. Może się je- 
dnak zdarzyć, że podobny wyraz oprócz masy ciała zakłócającego 
zawiera jeszcze mały dzielnik rzędu 0, otóż Brendel!) nazywa 
wyrazami rzędu istotn ¡e n-go takie wyrazy, których wartość bez- 
względna może być porównywalna z n-tą potęgą masy m'. Np. 
w równ. (7) zwykłe wyrazy (a także ich współczynniki b,) są isto- 
tnie 1-go rzędu; charakterystyczne wyrazy (także i spółczynniki P) 
nie są istotnie pierwszego rzędu; natomiast iloczyn 5, X B, jest 
formalnie drugiego rzędu, ale tylko istotnie pierwszego rzędu; 
idąc za Brendel'em oznaczać będziemy przez litery łacińskiego 
alfabetu wielkości istotnie pierwszego rzędu, zaś przez litery grec- 
kiego alfabetu wszystkie te wielkości, w których występują dziel- 
niki rzędu б. 

3. Przejdziemy teraz do zbadania ruchu planetoidy w jej 
chwilowej płaszczyźnie drogi. Gylden zakłada, że protometr Or- 
bity a (w naszym przypadku jest to półoś orbity) jest zdefiniowany 
przez takie samo równanie, które daje nam związek w ruchu Ke- 
plerowskim pomiędzy a i n, a więc, 


(1) ан 


a? 


w ruchu Keplerowskim. Według Gyldén'a a musi też takie 
samo równanie formalnie spełnić, tylko w tem równaniu n nie od- 
powiada już ruchowi średniemu po elipsie Keplerowskiej, jest 
raczej pewną stałą całkowania, którą Brendel nazywa stałą ruchu 
planetoidy, i która, jak zobaczymy, pozostaje w związku z wiel- 
kością n, , oznaczającą ruch dzienny średni. 


ү de ADI o A 07 Maik 
oniewaz równ. „= 118 do "Vad mm 
więc możemy napisać: | 

dt (1 — m: 
(2) n dv a" (1 -|- p? (1 + Š) 


równanie to będziemy musieli całkować. 
Wielkość p Gyldćn rozbija na dwie części (p) i R: 


(3) eh 


') Brendel. Theorie d. LIMPIEZAS 


<= 
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w ten sposób, że w (p) są objęte wszystkie wyrazy typu (B), a P 
jest funkcyą pierwszego rzędu i zawiera wyrazy charakte rysty- 
czne oraz najwieksze zwyczajne. Równanie na (p) możemy 
napisać w postaci następującej: 


(p) = К cos [(—s)v—T]-- Xk, cos [(1 —0,)°— 14] 


albo 
(4) (p) = К cos (v — o) -|- X cos (v — w,) 
gdzie 


w == Г + pv; ln + Dell i). 
Wprowadźmy do równ. na (p) funkcye 1 ill przy pomocy 
związków następujących: 
1 cos [I = k cos w + X kn cos o, 
y Sin I= k sin w + XE, sin w, 
z tych związków z łatwością znajdziemy 
(5) (p) =qcos (v — ID), 
gdzie q i II są to funkcye zawierające wyrazy elementarne diugo- 
okresowe. Gdy funkcye y, П, R są już znane, to możemy przy po- 
mocy następujących związków obliczyć promień wodzący 7, mia- 
nowicie równanie str. 4 jest: 
_a(l—?)_ а(1— 3) а (1 — 1) 
e  1+0)+R E EECH 
gdzie v=v— II. 
Z równania ostatniego wydzielmy znów pewną część 7, ana- 
logicznie do tego cośmy zrobili dla p, mianowicie 


(4) 


(6) (3) = Pd | : 


Wielko$é (r) w ten sposób zdefiniowana jestto według Gylden'a 
„absolutny promień wodzacy“. 

Röwnanie (2) z poprzedniego paragrafu daje nam związek 
pomiędzy czasem ź i miejscem plauetoidy na jej drodze; to równanie 
przekształca Gyldćn, posługując się analogią wzorów dla zmien- 


1) p jest rzędu masy zakłócającej i charakteryzuje ruch apsyd. 


O ruchu planetold. o 


Шымы с 


18 J. KRASSOWSKI. 


nych ruchu eliptycznego i tych nowych zmiennych, które on wpro- 
wadził jako charakteryzujące ruch planetoidy. W tym celu Gyl- 
den wprowadził nową zmienną e, analogiczną do anomalii mimo- 
środowej w ruchu eliptycznym, która jednakże występuje nie w wy- 
rażeniach promienia wodzacego r, lecz absolutnego promienia (r). 
Posługując się wspomnianą analogią kształtu wzorów dla ruchu 
Keplerowskiego i zakłóconego, znajdziemy 


| (7) (r) —a (1 — 4 cos e) 
| stad 


1-1 cos V 
ln cos v 


— "|? 
l+ncosv 


NS /1 --1 | € 
c) шо =] Eam. tg. 

Podobnie oznaczając przez M zmienną analogiczną do ano- 
malii średniej w ruchu eliptycznym, mieć będziemy: 
| (9) M=e—nsine. 


Posługując się formalnie wzorami ruchu eliptycznego, rozwi- 
niemy anomalie M w szereg ?): 


a) sin e = sin v; b) COS === oraz 


(8) 


В, =— 21; B=in ++... В, = — {1 — р>; 


5 EZ BL. U : | 


М = — qsin e = v D sinn v, gdzie 
(10) | 
n 


Różniczkujemy pierwsze równanie z (10), uważając, iż 1 jest 
stałą : 
1 (de—qcosede =dv--ZnBncosnvdv, czyli 


1 | 
(11) | (1—1 COS €) 3 =1+%, B,cos nv; 


] — y 
PCOS e Торе 
różniczkujmy ten związek względem у: 
| ln  . V1 — ai 
sin € GR 1 sin V, lecz sin є == - br 


dv (i-+ncosv)? I+ncosv 


| na. Br e wde SET ed kl ZB Sp +30: 
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podstawiając te wartość otrzymamy: 
7008 Pd ds op И1 — 1 
ау  (l--«cosv) dv 1--qcosv' 


Obie strony (11) (а) pomnóżmy stronami przez otrzymane równa- 


ABE 


ie dl ge trzymamy: 
nie dla š otrzy y: 


3 
de — PN! — 12)? 
dv (I-+- cos v?) | 
co znów podstawiając w równanie (11) da nam rozwinięcie nastę- 
puiące: 


(1 — cos e) 


(12) Са) ayia 1 + En B, cos nv 

(1-]- «cos vy? 3 | 
Równanie (12) jest bardzo ważne a to dlatego, że dzięki samej defi- 
пісуі Gylden osiągnął to, iż po prawej stronie (12) niema żadne- 
go wyrazu długookresowego, podczas gdy z lewej strony znajduje 
się funkcya o długim peryodzie. 

Zwróćmy się obecnie do związku pomiędzy anomalią v i cza- 
sem L mamy równanie (2), р. 16; w tem równaniu n jest to tak 
zwana stała ruchu planetoidy. Niech Z oznacza długość średnią, 
analogicznie do definicyi w ruchu eliptycznym, będzie ona okre- 
śloną w sposób następujący: 


(13) L=ni--A, 
gdzie A oznacza długość średnią w chwili gdy t=0; wtedy będzie- 


my mogli zdefiniować nową zmienną czas zredukowan у (t). 
„Czas zredukowany“ spełnia związek następujący: 


(14) п) + A=v+ Z В, sin nv. 


Na zasadzie znanego wzoru z trygonometryi możemy napisać 
uwzględniając wartość v: 


B, Sin nv = B,sin nvcosn II — PB, cos nv sin П, 
skad otrzymamy: 


dt) _ 1 dE 
(14a) n q, 7 1-2 nB, cos n v + zę 
de  sd(B,ecosnH) . ` a (B, sinn П) 
gdzie > До HUE ~y  COSnV 
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jest to zresztą mała wielkość, którą przeważnie można zaniedbać 
(w naszym przypadku nie możemy tego zrobić). Ostatecznie więc, 
uwzględniając poprzednie związki, otrzymamy: 


(15) „0. U ux ny 
SE ао (1-- «cos v? ' dp 

Niech W oznacza róznice pomiedzy katem, odpowiadajacym 
prawdziwemu czasowi nt i kątem odpowiadającym czasowi zredu- 
kowanemu n(£)—mieć będziemy : 


(16) W = ті —n(t), 
a ponieważ 
L=nt+A 
zaś 
A=v—n(l) ++ XB, sin nv, 
więc 
L=nt— n(t)4- v 4- ХВ, sin nv, 
czyli 


L=v+XB, sin nv + W. 
Wynika z (16) 


aW d At) 
do w ао 
ale dla n T. mamy wyrażenie (2), zaś dla n wyrażenie (14a). 
co uwzględniając otrzymamy: 
3 ] Е 
Aga elo dd cogo | W x rd 2 B 
= {л RB. | 
do sry) 5929) | | | PD ere | dv 
Rozwijając wyrażenie (17) według potęg £ 1 S, otrzymamy : 
aw 2 à 
E — S— 2R—2RS + ЗЕ? IL 3SR*—4R* +... 
+ (62 — 28 — 1242? - 6 ES 4...) cos v 
VI 


-RHS S—6R+... j cos 2v +... 


9 | 
-- 6Ry° cos v +- e R — s) 7 COS ЗУ 


gdzie v=%—Il. 
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Cała metoda Gylden’a jest zawarta we wzorach IV, IVbis V, 
Visi VI. Zajmiemy się teraz zcałkowaniem tych równań i wyzna- 
czeniem wyrazów charakterystych i długookresowych w przypad- 
ku ruchu planetoidy Thule. Nim jednak do tego przystąpimy, 
musimy Sie zapoznać z rozwinięciem funkcyi perturbacyjnej 
w spółrzędnych Gyldénowskich. 


II. Rozwinięcie funkcyi perturbacyjnej. 


W równaniach otrzymanych poprzednio dla Sip występnją 
wielkości P i Q, które zawieraja pochodne czastkowe funkcyi per- 
turbacyjnej Q. W funkcyi perturbacyjnej występuje część 


1 
5 


] ną 
к= (r?r? — 2rr' cos H) ° , 


wskutek tego nie możemy całkować powyższych wyrazów bezpo- 
średnio. Aby całkowanie wykonać jesteśmy zmuszeni rozwinąć 


A na taki szereg, którego wyrazy można całkować. Gyldćn roz- 
wija funkcyę perturbacyjną na szereg trygonometryczny według 
wielokrotności kąta H i potęg wielkości r i y”. Gyldén zakłada że: 


(1) nn o-F 2E, cos H|-2R, cos 2H+-... 


W tem rozwinięciu spółczynnik №, wyrazi się przy pomocy 
wzorów następujących: 


! nu da f | 
dire E cos nd dd > | cos nydy 


| a 
E VyhLy? 9pp'cosQ т ly (ryż „т | 
0 (RN 1445) —2|-,|cosę. 
| Y Y 
a . | TER 
Oznaczając a = т! wprowadzając nową wielkość zmienną A, 


bardzo małą, określoną przez równanie: 


a 2 y 2 
(3) (=; 
otrzymamy 
теү - 


i związek (2) przejdzie w następujący: 
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RN f cos фа № 
rm y) TE. E OI. . 
w V 1 +- k? — 2kcosy 
Jeżeli do tego ostatniego wzoru zastosujemy znane prze- 


kształcenie z teoryi całek eliptycznych, to w rezultacie otrzy- 
mamy: 


мя 


| 2 (y 2n 
(5) R, -— 2 : = on" (12)? sin pa 
И] — a? sin 20 Loi sin? p 


0 


W całce (5) wyrażenie pod pierwiastkiem można napisać 
w sposób następujący: 


2 Ama f. is: sin^" ed | 
T VI —a?sin?p + а? A ѕіпф — 
(6) 
sin?" p d ç 
анн 


|-ażsinżę 


gr! 


4 
т 


Gyldén, mając już wartość każdego R z wyrażeniem (5) 
rozwija је według potęg A. W ten sposób znajdziemy : 
2а XS f. sin "edo | ti 
% Su man Nr {n.1 À Тн ^^.. 
T J V1— o? sin?p(1—A 


gdzie za Gyldén'em?!), oznaczając przez 


т 


P] 
2 (6,2 SAW (10 
8.0) — EC ER 
(8) P 1 | {1 —a?sin? p] 5 
otrzymamy 
(1) ] (3) DAS 
| fn.0 = an+! В»; йк ү 9 атту Вы ИЕ 8 art B.A 
d J, SEE 1) (20-41) 
5 v DE = El 
А ZE : n+29+1 
| [15079 ТЕ: NOS ps: Pete 


1, Gylden H. Undersökningar etc. H art. 38. 
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Spółczynniki 4,. zależą więc tylko od а. Wprowadzmy 
otrzymane wyrażenia z (7) i (9) do wzoru (5), wtedy znajdziemy: 


i a R 
(10) Fin = pi (1— А)? {о — fai A -H Yn. тиз? H.. 3! 
Rozwiniecie (10) jest jeszcze zbiezne i w przypadku planeto- 
idy Thule. Podstawiajac wzór (10) w równ (1) otrzymamy żądane 
rozwinięcie tylko części ^ funkcyi perturbacyjnej— wiadomo bo- 


wiem, Ze cala funkcya perturbacyjna!) ma kształt nastepujacy: 


— | 9 Of [Es | KO ra р ЖЬ | 
(11) ai — m! [^ E cos H}= m |+ Sat À)? cos Н 


Rozwińmy w szereg wzór (11); w tym celu uwzględniając 


WEIT a 
rozwinięcie dla A otrzymamy ze wzoru (11) wzór następujący : 


a? 
9 


(12) a(Q) = т!к, +2 |, — E (1— А) | cos H+ 2R, cos 284-..] 


możemy ogólnie więc napisać, że szukane przez nas rozwinięcie 
funkcyi perturbacyjnej mieć będzie kształt następujący: 


(13) a(%) =Q, +29, cos H+ 28, cos 2HL,..20, cos nH — 
= 0, +49, cos n H, 


gdzie 
(14) OQ, — m S (1 TEM A)? Ion — m1 À Ta 2 À? — зде; .. | 


Two Oznacza, że dla wszystkich wartości n, z wyjątkiem n=1, mieć 
będziemy 
2 
o: 
Otrzymaliśmy więc rozwinięcie funkcyi perturbarcyjnej w po- 
staci, podanej przez Gyldćn'a. Jest to jednak jeszcze postać 
nie ostateczna, gdyż w rozwinięcie (13) zamiast wielkości À musi- 
my wprowadzić wielkości р, p”, n°, 2€. 
Związki pomiędzy v, r’ i 1, d Sa następujące: 


Yno = Yno, Zas dla n=l то 4,9— 


„Um, A т). 
agii eg. Le. 


') Z pominięciem zaburzeń samej płaszczyzny drogi. 
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a więc 
y y2 42 1922 112 
de Pzd ph s EE 
| al \r 1+p 
Podnosząc A do potęgi od 1-ej do 3-ej i uwzględniając wyra- 
zy trzeciego rzędu włącznie (a więc w krórych mamy pš, p", p”, 
gin np, etc.) znajdziemy: 
À = 2p — 2р'— 3p2+ App! — p' 24-29 — 224 4p! — 6p?p' + 
+ 2pp'? — pv? + 4p^n* -t 4р1? — 40 1. 
= 4р?— 8pp'+ 4p? — 12p°-4- 28p?p'— 20pp *-1- 4p" 1- 8P + 
+ pn? — Spin? — Spin”. 
A3= 8p3 — 24p3p' — 24pp'?— 8p'3. 


(16) 


skąd na > (1 Ba otrzymamy związek następujący: 


a " (1 + p^)! ( ] 2 т ү" W 
E À = — = | — 
(17) (1 ) AENA wën? lI—np+(n+ lp + 


n(n-- 1) , | |, NANT D), À " 
T (ąz ) *— n (n + 1)pp T T p?—nq + (n+1)7?- 


| 9 2 2 ] 
er ) e O "POLA d dole 


UO — Dni T 1) 
| 6 
-- n (n -F 1)’. 
Otrzymany zwiazek (17) wprowadzmy do wzoru (13), wtedy 
otrzymamy ogólnie wzór następujący: 


p'3-- п?рт?— n(nJ-1)p^r?—n(n Lef? 


(18) a (Q) =2m' » Q (n, S, 8 рер" 127 27 cos n H. 

W tym wzorze należy pamietać, że dla » =0 mamy opuścić 
czynnik 2. W ten sposób ostatecznie otrzymaliśmy rozwinięcie 
postępujące według potęg v’, 7/?, рі р dla całej funkcyi perturba- 
cyjuej. 

Spółczynniki ©, wyrażają się przy pomocy związku następu- 
jącego: 


v D ani] 
(19) 0 = m' y mas р? А я 
kd 
n, 8, 8! v, v' 
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Te spółczynniki są to funkcye liczbowe znanej nam wielko- 
$ci a, a więc nie trudne do obliczenia. Harzer!) obliczył licz- 
bowe spółczynniki w wyrażeniu Q,, zaś Gyldén?) podał warto- 
ści tych spółczynników dla wielkich planet aż do wyrażeń zawiera- 
jących a’ a także i liczbowe wartości stałych niezależnych od a, 
które wchodzą w wyrażenia na ©. 

О... уи W funkcyi spólczynników liczbowych ү. wyrażają 
się w sposób następujący *): 


(0 оо = Yn.0 
Олло == — NYn0— MT? 
n (n--1 
Q,20 = GE {not (20|-3)1n1-h 4124 


0,4111 (N SI 1) (n + 2) Va+1.1 — 2(2т-|-5)т, ал z BY +-2.2 
Oui Z—(R—1) mn 41,349 — Z(2N-L)Tn 211 — 871 1.2 


(06.1.0 = — (n -- Dena gs DIAM 
Q4 110— — (n 1)їта--1л — 218—312 
(20) zu za = Зтол + 41o2 
n + 1) (n + 2)? CR 
О ойе QE 2 RA Tata D An — See + 
+4(3n +1 1)Yn41.2 -|- 401.3 
(n+ 1) (n + 2 
Q „+1.2.0 = | d ks "(n41171- (20-5) Teil + AYn+12 
n? (n —1 
Mn AH io On! Hn-H2)1. a + 
Së 4(3n T O)Tn—12 -}- 2470-13 
= 
Oe оз ms a -|-(2n kn dps HAE: etc. 


( 
W równ. różniczkowych па Si р wchodzą wyrażenia pocho- 
dnych cząstkowych funkcyi perturbacyjnej względem r i v--a nie 
bezpośrednio rozwinięcie samej funkcyi przeto obecnie, zna- 
jąc wyrażenia rozwiniętej już funkcyi perturbacyjnej możemy zna- 


1) Harzer P. Untersuchungen pg. 22. 

*) Gyldén. Orbites absolues t. I: 

3, Brendel. Th. d. kl. PL I, pg. 52; Harzer.’ Untersych., pg. 20; 
Kramer. Th. d. kl. Planeten, pg. 9. 
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leźć i związek pomiędzy P i Q oraz funkcyami Q,,,/5; w ten spo- 
sób wielkości P i (Y staną się znanemi funkcyami liczbowemi 7,., 
co zatem idzie a, a więc ich obliczenie, choć kłopotliwe, dla danej 
planetoidy nie przedstawia zasadniczych trudności. 


Mieliśmy równ. następujące: 


EI y? EI 
шило”. gx | | 
Te =} ży oraz © EIER 
Ponieważ 
@(1— 1°) 
Y == — 
leri n 
to znajdziemy: 
Q e? 29 ] — | 90 


al) 3v ^ (1--p) än 
W tem wyrażeniu pierwszy czynnik rozwińmy w szereg po 
tęgowy, potem zróżniczkujmy względem v; otrzymamy wtedy 
Q= — (1 — 2p + 3p?—...) 2Èn0 sn PE 2 y 2" sin m HOS 
(bo wiemy, ze H zależy od v) co możemy napisać: 
(21) qum A E PEN sino HAS 


(dla n = 0 należy opuścić czynnik 2) — gdzie przez Oe. ozna- 
czyliśmy pewną funkcye Qss, a więc funkcyę wielkości 7,:, któ- 
rą można otrzymać na drodze rachunkowej; postać tej funkcyi jest: 


] 
(22) m u LÉI, Mai - 9, TE 2 20. -1.8'. v. y! d- oM, я— AR, vv 
/ 


Zupelnie analogicznie znajdziemy wyrazenie 


| ,90 „ 90 Q 
Ps dr = — (1 n°) әр ' 


Po wykonaniu działań i podstawieñ wskazanych, jak to 
zrobiliśmy dla Q, otrzymamy rezultat następujący: 
(23) ОРООД ap OS MOM 


gdzie przez P,,,,, oznaczylismy wyrażenie następujące: 


l | | 
(24) ТТ, ia AS xcv m: (S - y 1), 8-18, лу" - (S zi ` IPS e: 
D 
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Wielkości P są tylko funkcyami liczbowemi wielkości Man 
a міес лс, przez co także te wielkości są znanemi funkcyami a. 


Brendel'), Harzer?) podali wartości Pns i Qus: 


i NOGA OFE EE S ligne, Mega 
Pa.1.0 = — 29,20 (Олло = Олан — 2 Ооо 
Pazo = — 39,30 Омо = Qazo — 281.042 Mann 
Jo MERI: | = — Qapi Qn+1.01 => Qu 1.04 
P. уол = — HA 111 On—1.0.1 — H. 104 
Pu-100= — Q, von (n--1.0.0 ZE (0, 1.0.0 
P o dee 200 (1.0.1 = Qio 

(24) d Dion v mer an = Mine 

| Dun = — Quo Qs 1 — CERE, 33 29, 1.0. 
Pea = 0 Qatı 2.1 ОИ ЕТ ZA 1.01 
Pur Z — 294842 Onyi 1.0 Z Фал — 2241.00 
Pn. = — 2120 In-110 = Q, 140 — 294 100 
den PEB AA age Quia A MEL e 29104 
Pozo =—3%.30 Olo ein 281.00 
Pins =— 28a. 

| Dun Z— 2 Q120 


W wyrażeniach na P i Q występuje kat H. Musimy obecnie 
ten kąt wyrazić przy pomocy prawdziwych długości v i v’ plane- 
toidy i planety zakłócającej; jak dotąd tak iteraz rozpatrujemy 
tylko ruch planetoidy w płaszczyznie drogi chwilowej, pozostawia- 
jąc do dalszej części niniejszych rozważań, wyznaczenie wyrazów 
funkcyi perturbacyjnej, zależnych od nachylenia i zmian położenia 
płaszczyzny samej drogi. W ten sposób znajdziemy po zaniedba- 
niu drugiej potęgi wstawy kąta nachylenia pomiędzy drogą plane- 
toidy i planety zakłócającej: 

(25) ?) Bzy 
Po wprowadzemiu tej wartości H w wyrażenia na PiQ 


otrzymamy wzory, w których występować będą obecnie długości 
prawdziwe v iv. W równaniach, z których możemy otrzymać 


1) Brendel, I. cit., pg. 56. Harzer, I. cit., pg. 30, 31. 
э) Masal H. Formeln u. Tafeln etc., pg. о. 
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S, p, W—zmienną niezależną jest prawdziwa długość v — więc aby 
módz zcałkować równ. dla P iQ) musimy wyrazić v przez v. 
W tym celu zwróćmy sie, idąc za Brendel'em, do równania (13) 
pg. 19, z którego otrzymamy związek pomiędzy położeniem plane- 
toidy na jej drodze i czasem. Znane wzory ruchu Keplerow- 
skiego dadzą nam dla planetoidy: 


(a) nt 4- A =v — 24 sin v + in sin 2v —...--W 

i analogicznie dla planety zakłócającej: 

(b) n't -- A'—v'— 21! sin v'--5/,1' sin 2v'— ...--W" 
gdzie v=v— П, zas v'—v' — I. 


Mnożąc wyrażenie (a) przez „=> otrzymamy: 
(26) nt=uv—pA—2py sin v ¿pon? sin 2v4- ... Wu. 
Z równań (b) i (26) przez porównanie znajdziemy: 
v'= (po — pA — Jun sin v --?/, ux? sin 2 v 4- ... Wy) 
+ A+ 29 sin v'4-...2-W'. 
" Wprowadzmy oznaczenie następujące: 
(07) B=A—uA i G=u2B,sin пу — ÈB, sin nv 


to otrzymamy: 


(28) v=w4B+@-+eW—W' 
a wiec: 
(29) v—v—= (1 —p)0— B—G — VW =w —@. 


(Możemy opuścić W' jako małą wielkość, gdyż w tym wyra- 
zie jest zawarty wpływ Saturna na ruch Jowisza, zresztą ta wiel- 
kość jest nam znaną z teoryi Jowisza), gdzie: 


(30) w, =(1—n)v— bM | 


Połóżmy dalej | 
(31) v' = — w, + r+ V, 1) 
wtedy 


(32)  G — — 2unsinv—+ 27 sin v,'+ 1&2? sin 2v-- 22" sin 2v, + ... 


) Brendel. Th. d. kl. Pl. I, pg. 63. 
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Aby przedstawić cos n(v—v') jako funkcya v, rozwińmy 
przedewszystkiem, uwzględniając związek (31), wyrażenie d sin v, 
w szereg według potęg 1': 

(33) w sin vy =—y' sin [w, — G — vj =—1/ Sin (w,--V,)+ 
-|-q'G cos (W, — Vi) —... 
gdy w (32) uwzględnimy wyrazy tylko do pierwszego stopnia 
włącznie, to otrzymamy: 
G= — Am sin v — 21) sin (W, — Vj). 

Podstawiając tę wartość w (33) znajdziemy: 

y sin v,'— —* sin (w, — V,) — 2enn Sin v cos (w,— dry wa 

— Ou ai? sin (W, —V,) cos (w, — V) =—1 Sin (W, — Vi) — 

— en sin (V— Vj + Wy) 9097 Sin (Wy — re Ya! 
RY sin (Z Wy Ve : 

Analogicznie postępując otrzymany dalej: 
EI, IS AW v 9) - 13? sin (wi — Vi) +... 

Podstawiajac te wyrazenia we wzór na G, otrzymamy w re- 
zultacie: 

| G=—2u4 sin v—2 sin (w, — Vi) + ¿pon? sin 2v — 

(34)! —21 (1j sin (W,-- v —v) 2-2w* sin (w,—v — Vi) — 
ox sin (2Ww E 2v) re 
Poniewaz ze wzoru (29) mamy 
v—v=w,—(, 
więc 
| n? G2 
cos n (v — v') = cos n w, + nG sin n w, — a COS RW, +... 


Mozemy wiec, uwzgledniajac wyrazy do Ш stopnia włącznie, 
wstawić wyrażenie (34) w to ostatnie rozwinięcie. Ostatecznie 
otrzymujemy cos nv—v) w postaci szeregu, w którym zmienną 
jest tylko v. Dla sin n(v —v') analogiczny wzór mieć będziemy 
z tą różnicą, że zamiast cos wystąpi sin. 
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Będzie więc: 
( COS Tee av =. OS 
E. 110—0 yes nen W- LUN QU Уг у) NOZ 
, cos cos 
qw c fm—Dwid-v] +k Ae [@-Е1)эл— 
+ Кт? COS zw + k > (mw. 12у) 
(SNU M Li sin 
, COS 
(85) 4-е т? (mw, — 2v) + km LED wi vil 
ho mm a [n4w vv, Lien! WER wi —v--vi] 
cos 
т: CES L)w,— vie? ^W, 


cos 
ETE ми [(n — 2)w, + 2v,] gg, ` T > [(n + 2) ¡—2v,] 


gdzie: | 
k =u; ką =—nu; kę4==— mh, =--h; kę ==—n u"; | 
MAS na: ale obla OUS PONI | 
к= 9 grep ERR 9 | g "^P | 
ke —n(n-—1)p; + n(n-- 1); ku =+nn DÉI | 
(36) Fun 
k = —nqin--1)p; k,,——mn* k E lg п |: 
11 о, ) "13 E 8 J 
ULP M 
KEE 


Gdy już wykonaliśmy rozwinięcie z n H według potęg v, 


to wtedy wyrażenia (21) i (23) pg. 26 napisać będzie można w po- 
staci następującej: 


z 

(36a) P == 2 > Be: s')vy' p° as E nas COS n (v = v^, | 
Oraz y | 
(37) (9 — — 28 Qin sw pp m2 2" sin n(v —v!). 


* 
(Y oznacza, że dla 2—0 należy opuścić czynnik 2). 
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W tych wzorach wydzielmy z wielkości p część elementarną 
kształtu b oraz zwykłą wraz z charakterystyczną; pomiędzy temi 
róznemi częściami p mamy, na zasadzie poprzedniego, związki na- 
stępujące: 
(38) p—=ncosv4R. oraz р (р) +R, 
gdzie, jak zawsze, (p) oznacza część elementarną p postaci B; zaś 
w À są zawarte wyrazy zwykłe i charakterystyczne. 

Ponieważ: 


Pl) + RI) ESO RE (o) - R’+.. Re, 


wiec gdy w tym wyrazeniu na p uwzglednimy tylko czesé charakte- 
rystyczną i podstawimy w (36a), to otrzymamy rozwinięcie postaci 
następującej: 

| P= P, + PRP, R+ P, R+... P EH... 
(39), gdzie A 

| P,=2Y Р, ыу, (9) . (Pen cos (v — v). 

Dla Q, mieć będziemy analogiczne wzory, tylko w tych wzo- 

rach zamiast cos (v—v') wszędzie znajdować się będzie sin (v—v'). 
Zatem według drugiego wzoru ze związku (39) mieć będziemy: 


+ x 
A ES Р (n.0.0)0.0 P ,.1.0)0.0 A Де (n.0.1)0.0 


-4 П 
AG ()(.0.0).0.0 +5 Om .1.0).0.0 (P) + — a ning ‘Ing 


п.о. 7 ёр .1)0. x: d 
+ y P ( dh Gg ` ( ka up S (n.0.2)0.0 e 1 


Oe 2.0). À — Сз 1.10.0 — Q(n.0.2)0.0 


Lis S ODU ud y Р (n.0:0)0.1 biles | i ld) 


dmi lang pum wma 


ut —29 S P («100.0 a ju" Pit "n Vide 1.1)0.0 x 


(0; a S oed Ys SA e TA Q(.1.1)0.0 


X RC | "nie it. d. 


Aby otrzymać ostateczną postać rozwinięcia funkcyi P i Q — 
musimy wykonać jeszcze jedno przekształcenie. Każda funkcya 
S, R, W, zawiera bowiem część zwykłą, krótkookresową, oraz dłu- 
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oookresowa, wyrazy długo i krótkookresowe elementarne, takież 
wyrazy charakterystyczne i wreszcie wyrazy zwykłe. 

W ten sposób funkcyę W możemy napisać w postaci nastę- 
pującej: 
(40) ` W то ЕМ, М-М. 

Część yu jest to część wiekowa. Wa oznacza wyrazy długo- 


okresowe, W;—wyrazy krótkookresowe, W. zaś wyrazy zwykłe. 
Według Brendel'a oznaczmy 


(41) Wı-- М, = K. 
to wtedy 
(49) wun Wa+ K, 


ale we wzorze (30) pg. 28 mamy 
(43) nw,=n(1—u)v—nB—nuW=n(l —&)o—nB-—nywW ,1—nyK, 
gdzie dla któtkości oznaczyli$my 
| RT ipte 
Jezeli jeszcze oznaczymy 
(1 — pa) v — B—pW¿=0, 


to wtedy !) 
w, =w EX K. 
Wtedy dla 
sin NW, 
otrzymamy wzór: 
; ' n*u? КЫ | 
(44) sin mao, — sin nw — ny. K соз nw — 5! sin nw... 


Jeżeli teraz uwzględnimy wszystkie wartości ШИ weas, 
i złączymy razem wyrazy zależne od tego samego argumentu, 
oraz uwzględnimy związki (44), to otrzymamy w rezultacie rozwi- 


nięcie w szereg nieskończony według potęg 1, т, R, К czastko- 


wej pochodnej P funkcyi perturbacyjnej, w którym występować 
będą tylko argumenty V, Vi, Oraz W, i pewne liczbowe spółczynni- 


ki (m.s.s')v.v' zależne jedynie od = — a, które dla każdego po- 


szczególnego przypadku zawsze będziemy mogli obliczyć. 


5 Brendel, I c. pg. 69. 
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W rezultacie więc otrzymamy dla P rozwinięcie!) podane 
przez Brendel'a w postaci następującej: 


(45) P=ZB, go cos mot ptt? ncos(nw+v)+ DES cos(nw -- v) 
+ ХВ 0 у cos(nw-rv,) + X ВС!) w cos (nw — v) 
+ > B, 2.017? COS nw 
FEB.) 1 cos (nw--2v) + SBE? q?cos(nw —2v) 


+ SB mm’cos(nw+v+v, HEBEN gy cos(nw-Lv—v, ; 


—1 


+ LR, Inn’ cos(nw—v-L vi) HLB "mcos (nw=v—v,) 
+12 B».0.2 1 ?COSR w 
+2 В, о '2с05 (nw-F-2 v,)H- X BC) Sege (nw. 2 v) 
TP ME 

 RAEXBL,cosmw- XB. T cos (nw Lu 
JB. O? 1 cos (nw — v) + > reg d cos (nw -|- v) 
+B. cos (mo —vi) +...) 

Luk À >n B ao Sin NWA Хт JENA y Sin (nw--v) 
+ Xn B; jn sin (nw —v) 4-2mB чү sin (nw-v,) 
+En BU) чү sin (nw—v,) +... ) 

+ В2 { Ў Bi, cos nw-|-...) +e RK Zn Во sin nw...) 


2 2 
Doy eH m. B, оо COS nW Y ) 


1) Brendel, |. c. pg. 70. Harzer. Untersuchungen etc., pg. 30. 

2) Wskaźniki przy spółczynnikach Ai B są tak dobrane, że u dołu 
pierwszy wskaźnik odpowiada czynnikowi przy w, następny daje nam wykła- 
dnik w, trzeci zaś wykładnik y. — górne wskaźniki dają czynniki przy р oraz 
wykładnik przy E (conir. Brendel I. c, pg. 71). Dla wartości n=0 niektó- 
re ze spółczynników A i B są zerami (contr. Brendel, 1. c., pg. 71 i 79). 
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Dla Q zaś mieć będziemy analogicznie: 


(46) Q=} Anoo sin nw4-X AC? 4 sin(nw+v) +42 АСУ) ysin(nw v) 


n.1.0 n.1 0 


val . ` — . 
kozy а a sin (nw + v) LAT a’ sin (nw — vi) 


п 0.1 N 
+YA ,, "^ Sin NW 
-- X40 y sin (mw--2v) HE ANS n° sin (nw — 2v) 


n.2.0 n.2.0 


LA my sin (ишу +v) T 2 А, na sin (nw 4-V — vi) 


A, 
SACD sat si (nw— ASAS m! sin (m: 
+ ZA, т sin (mw— v +V) r24,,, m SIN (nu v — Vi) 
+ ЖА ,, 1^ sin nw 


sA? тї? sin (nw-F-2v,) -ЕЎА у, 0? sin (m0—2v,)-|-... 


m.0.2 


L RAE ALS, sin nw 4-2 ARO sin (nwy) +2 А ту nsin(nw—v) 


n.1.0 


+ SAT" w! sin (mw + v4) + A ad cos (nw — v) +... \ 


n.0.1 


—pK | Xn 4,00 COS nu + Xn AGV a cos (nw + V) 


n 1.0 


+ yn AC y cos (nw—-v)-- In dt vd cos (n 104-v,) 


+ Xn АС? cos! nw—vwv,)4-...j 


m.0.1 


+ RZ { XAM sin nw...) — ВЕК { Уп А0 cos nw | 


2 1 
= 4 > o Anoo sinn W | > 

Wzory (45) i (46 rozwinięte do 3-go rzędu włącznie posłużą 
nam do wynalezienia wyrażeń analitycznych dla ruchu planetoidy 
Thule (279). Wyrażenia (45) 1 (46) będą wystarczające dla nasze- 
go celu, ponieważ nachylenie orbity Thule względem orbity Jowi- 
sza jest bardzo małe. 

We wzorach (45) i (46) występują spółczynniki A i B opa- 
trzone wskażnikami; te spółczynniki są w dość prostym związku 
z wielkościami Pass i Qu», które, jak wiemy, są to funkcye Quss w 
przeto spółczynniki Ai B ze wskaźnikami dadzą się z łatwością 


1) 1. с. pg. 3. 
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obliczyć jako funkcye liczbowe n.s, czyli wielkości o., znanej nam 
z elementów drogi uważanej planetoidy i Jowisza. 
Związki pomiędzy Ai B oraz Pi Q są następujące !): 


Ánoo = — 2n on 

Ao = — 9n Qus 

Ag! = — 27 Qu.2.0 

Año = — "4 Quio п Que o! 
An! Lp Qul 

Zen Sone ët Genen Bert DÉI uge 
Aż; SE — (n—1)19. Aan +2 (02-1) Q; o 


HEN "---4 120,20 + 27р. ol 


EE Ee 
dean = DEI naa — Е) Оо 


A2 TR (n— 1)i Qu 1341 + 2(n — 1) oam 
БВ» о.о — 280 0 


ВТЕ gp 

аар ЫБ, 

Bao YYY | 
za = Puio — 270 Paan 


Bin Banz An+1)Parıo 0, 
De = Pa-10.1 +2(n—1)P;_100 
Bi, = Рачо P, uj 
Bio =2Pn20— п Paso 


тети — пр 2(n-|-1) SARAT 
Eug Eeer 


n.0.1 


:) Brendel. Th. d. kl. Pl. I kap. V sparsim. 
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III. Wyznaczenie wyrazów charakterystycznych 
i elementarnych w rozwinięciach odnoszących się 
do ruchu planetoidy Thule. 


Nastr. 13i 14 podaliśmy ogólny kształt równań różniczkowych 
wyznaczających wielkości S i p; równanie dające W jest analogicz- 
ne do równ. na S, mianowicie W — C, -- X 2. sin (Anv Ball 

W powyższych równaniach wielkości a, b, с, sa to bardzo 
małe stałe wielkości; À,, jak to wiemy, jest mały dzielnik. Mo- 
że się zdarzyć, iż pomimo, że wielkości a, b, c, sa małe, je- 
dnak wskutek całkowania wystąpią jako dzielniki A,, które w nie- 
których przypadkach mogą mieć takie wartości, że dany wyraz, 
w którym wystąpią po całkowaniu czynniki = albo Ss moie 
się stać albo bardzo małym blizkim zera, albo też bardzo duzym 1). 

Gdy A, jest bardzo małe, blizkie zera; wtedy argument cosi- 
musa jest bardzo małym kątem, z czego wynika, że okres takiego 
wyrazu, zawierającego w sobie jako dzielnik bardzo małe A, jest 
bardzo znacznym. -— Wyrazy zależne od takiego A, będą wyra- 
zami długookresowemi. 

Gdy zaś A, jest blizkie 1, wtedy okres takiego wyrazu jest 
niewiele różny od okresu zmienności kąta v, a więc od czasu obie- 
gu planetoidy. Te wyrazy nazywać będziemy wyrazami krótko- 
okresowemi. 

Jest rzeczą oczywistą, iż okresy krótkookresowych wyrazów 
są blizkie okresu obiegu planetoidy. Jak to już poprzednio wspo- 
minaliśmy w wyrażeniach odnoszących sie do S, wskutek całko- 
wania tylko wyrazy długookresowe zostaną powiększone, podczas 
gdy wyrazy o krótkim okresie nie ulegną zmianie; odwrotnie zaś 
dziać się będzie, gdy całkować będziemy wyrażenia dające nam p— 
tam długookresowe wyrazy nie zmienią się, natomiast krótkookre- 
sowe będą powiększone. Brendel zauwazyl?), że w przypadku 
współmierności bardzo blizkiej pomiędzy ruchem planetoidy i pla- 
nety zakłócającej, wyrazy nawet drugiego rzędu bywają jeszcze 
bardzo znaczne. -- Z drugiej strony zasadnicza myśl Gylden’a 


) Brendel. Om amvándningen pg. 6. Ludendorff. Die Jupiter- 
stórungen d. kl. Planeten von Hecuba-Typus pg. 4. 
2) Brendel. Th. d. kl. Pl. I Kap. VII sparsim. 
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polega właśnie na tem, że w wyrażeniach przybliżonych, jakie 
otrzymujemy na P i Q porządkujemy wyrazy nie według rzędu ani 
stopnia, lecz jedynie według rzeczywistej wartości tych wyrażeń, 
które one osiągają w wyrażeniach całkowych. 

Dla tego też przedewszystkiem musimy w rozwinięciach na 
S, p, R, W zająć się wyznaczeniem krótko i długookresowych wy- 
razów dla specyalnego przypadku ruchu planetoidy Thule (279). 
W wielkościach S, p, W, R, zgodnie z Brendel’em, oddzielmy 
część długookresową, część zwykłą i część krótkookresową w ten 
sposób, że pod odnośną liczbą dopisujemy wskaźniki d, 2, k; przez 
(p), zgodnie zresztą z tem co poprzednio pisaliśmy, oznaczymy 
w wielkości p wszystkie wyrazy elementarne postaci D; te wyrazy 
zresztą, jak to już zauważył Brendel, są krótkookresowe ') (są 
one oznaczone przez (p)). Mieć więc będziemy: 


(1) S= Sa + Š + 5,; R= R+ Ri R. 
p= (р) + E; w=y0+ Wa4- K, 


gdzie K= И, + W,, zaś yv oznacza wyraz wiekowy. 

Przejdźmy do naszego przypadku specyalnego; stosunek ru- 
chu dziennego n’ Jowisza do ruchu dziennego średniego n, plane- 
toidy (279) Thule wynosi: 


PIESZY 30, 
(2) md 7398 — 4 
gdzie 0, jest bardzo mała wielkość. 
Zatem 
(3) Las „== 


a ponieważ ogólna postać argumentu, wchodzącego w nasze roz- 
winięcie jest 

(4) nw=n(l —u)v—nB—pnW;i=n(1l —i,)v—nB—nuWa, 
to w wyrazach krótkookresowych spółczynnik przy v w tym argu- 


mencie jest blizkim 1, w długookresowych zaś różni się on mało 
od zera. Jeżeli więc chcemy znaleźć w rozwinięciach na Рі Q 


» Brendel, І.с. 
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wyrazy krótkookresowe i długookresowe, to musimy tak dobrać 
czynnik 2, mnożnik przy v, aby argument cosinusa, względnie si- 
nusa odpowiedniego kąta był w pierwszym przypadku blizki je- 
dności, w drugim zaś blizkim zera. Uwazajmy zatem nasze rozwi- 
nięcie (45) i (46). 

Ponieważ w obu tych wzorach wskaźniki odnośnych spół- 
czynników Są te same, zaś same spółczynniki są inaczej oznaczo- 
ne, przeto chwilowo wystarczy nam w zupełności zbadać tylko 
spółczynniki zawarte w rozwinięciu dla P, gdyż spółczynniki roz- 
winięcia dla Q otrzymamy ze spółczynników P, zamieniając tylko 
symbol A na symbol B. 

W rozwinięciu (45) następujące wyrazy Leo rzędy 0,112 
stopnia będą długookresowe. W 1,6, 13 wyrazie składnik argu- 
mentu zależny od v jest blizki 0 dla nm =0, a więc długookresowe 
funkcye będą przy Bn Baas: Br B,,,; dla n=4 3-ci і 5-ty 
wyrazy są długookresowe. Odnośne tedy spółczynniki dla tych 
wyrazów będą: B,,,, B,,,. Podobnie 8, 12, 15 wyraz dla n= 8 
będzie blizki zera; spółczynniki tych wyrazów będą odpowiednio 
Bii Bari» Baoa. Wyrazy krótkookresowe tego samego rzędu 
co poprzednie, a więc I-go, są te, dla których argument zależny od 
v będzie blizki 1; te wyrazy są następujące: drugi i czwarty wy- 
raz dlan=0. Spółczynniki tych wyrazów będą; BJ i В! ; da- 
lej З i 5-ty które nam dadzą Всі і Ву. Biorąc n=4 otrzyma: 
my krótkookresowe wyrażenia dla I, 6, 9, 10, 12 wyrazu, dla któ- 
rych“ spótazynniki! bed BAY, "B, вив 3 BER Er 
Bi» Bios; dla n=8 argumenty powodujące wyrazy krótkookre- 
sowe występują w wyrazach 3 i 5-ym, dla których spółczynniki be- 


dą В." і B,,,; wreszcie dla n=12 takie same wyrazy występują 


w 8, 12, 15 wyrazie rozwinięcia, skąd otrzymamy: Rie Bus ; 
Bi Wyczerpalismy tedy wśród wyrazów Leo rzędu już 
wszystkie, które dla Thule (279) będą krótko i długookresowe. 
Wypiszmy te wyrazy: 
Wyrazy o długim okresie: Wyrazy o krótkim okresie: 


+i 
B 0.0.0 Boro т COS V 


2 av Pi; 
Baan? Doo, W COS v, 
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al), 


po 


cm qq cos (v — vi) 


(Gei т cos (v, — v) 


—1 
Epices 


\ 


(221) 
В 1 COS V, 
В, oo COS 4 w 


B,, 1° cos 4 w 


39 


ВС) y cos (4w — v) 


4.1.0 


Bios! cos (4w —v;) B, o2 1 COS Aw 


4.0.1 
` BR) qm! cos (404 v — vi) 


| 
| 
| (2: 
| Bzy 1° cos (8w--2V) 
BR nn! cos (8w — v — v) В! nm’ cos (4w — v + v) 


1.1 
Bar; n’? cos (8ш 2v,) Brad, 1° cos (4w — 2v) 


BT} 4? cos (4w — 2v) 


Ems m’ cos (4w —v — Vi) 


ВС" y cos (8w — v) 


8.1.0 


BC w cos (80 — v) 


ВСЭ 4? cos (19140 — 2v) 


12.2.0 


DER m’ cos (12w --v —v,) 


2.1.1 


ВС 9? cos (12w — 2v,) 


Analogiczne wyrażenia znaleźlibyśmy dla spółczynników A. 
Zajmijmy się teraz wyszukaniem wyrazów o długim 1 o krótkim 
okresach wśród wyrazów drugiego rzędu. W tym celu, jak to już 
zrobiliśmy poprzednio, wielkość R rozbijmy na dwie części, dłu- 
gookresową R. i krótkookresowa R,, tak ze R= Ra+ Rr. 


Ale R jest to szereg trygonometryczny, musimy więc rozwi- 
nąć oddzielnie cześć wyrazów drugiego rzędu, zależną od długo- 
okresowych argumentów wchodzących w E, t.j. zawartych w R, 
i część w której występuje E,. Postać tedy tego rozwinięcia, 
w którym zawarte będą wyrazy drugiego rzędu, będzie następująca: 


Ry, = (В,4- Ва) | Bos, cos nw-+ Brio n cos (mb ut, 
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Wyrażenie w nawiasie |... oczywista zawiera wyrazy o krót- 
kim i długim okresie, a także wyrazy zwykłe. Jak to powiedzie- 
liśmy wielkość Е jest to szereg trygonometryczny, w którym wy- 
stępować będą argumenty kątów, analogiczne do tych, które są 
w wyrażeniu wi |. Przez wykonanie mnożenia wskazanego 
powstaną nowe argumenty, będące sumami lub różnicami po- 
szczególnych pierwotnych argumentów. Buchholz podał 
praktyczne wskazówki, kiedy w wyniku ostatecznym, przy mno- 
żeniu dwóch funkcyi różnych (krótko lub długookresowych) otrzy - 
mamy wyrazy o długim lub krótkim okresie lub wreszcie wyrazy 
zwykłe. Oznaczymy według Buchholz'a przez x i 6 funkcye 
długo lub krótkookresowe, występujące w wyrazach szeregu, które 
mnożymy przez siebie i niech rezultat otrzymany będzie W, wtedy 
mieć będziemy: 


(2) | Ха. Ë= Va; Xa. = Vis ya Fat Fs; PEN Ee 
Dy. 6 — it; y, = Wa V EV. 


Zastosujmy te wzory Buchholza do naszego przypadku 
i wyszukajmy wyrazy drugiego rzędu w rozwinięciach (45). 

Przedewszystkiem rozważajmy tę część iloczynu R,,, w któ- 
rym wchodzi jako czynnik wielkość Ra. W tej części wynajdziemy 
następujące wyrazy Il rzędu: 


Długookresowe Krótkookresowe 
B. B^, cos 4 w 
COS CON AP veda ir u 
Ву т} cos (410 — vi) Bi" «cos v 
(3) A COS 
BA Т/ COS Vi 


B, 11? y cos (Bw v) 


8.1.0 


RY y! cos (Bw — v,) 


8.01 


| 
Dalej w iloczynie R,, musimy uwzględnić także wyrazy kry- 
tyczne powstałe wskutek mnożenia przez R,. Wyrazy długo 
¡ krótkookresowe będą przedewszystkiem wszystkie te, które już 


wyznaczyliśmy dla poprzedniej części (w której jako czynnik wy- 
stępowało Rig)»; oprócz tych wyrazów wystąpią jeszcze wyrazy po- 
chodne od wyrazów zwykłych, szeregu | | dla których argument 
jest równy 2v. Te wyrazy musimy uwzględnić, gdyż po pomno- 
żeniu przez №, one dadzą nam nowe, wśród których będą wy- 


razy o krótkim okresie. Oprócz więc powyżej już wypisanych 
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wyrazów (3) w tej drugiej części rozwinięcia dla R,, wystąpią je- 


szcze wyrazy zależne od następujących zwyczajnych wyrazów, mia- 


nowicie: 


1.0 Q 
В, Cos 840 


gro 


4,0 T COS (40 —+ v) 


Bio, w cos (4w v,) 


Bib: y cos (1210 — v) 


,1.0 


e q cos (1210 — v!) 


Wyrazy drugiego rzedu jeszcze powstana z iloczynu 
LK iEn Bros sin л w- ...[. 


Rozłóżmy znów K па części składowe K — К-К; na za- 
sadzie poprzedniego wiemy, że K¿=0, a więc pozostanie rozpa- 


trzeć tylko część 


„КУ Воо Sin w--...|. 
| następujące wyrazy będą te, które na- 


W wyrażeniu w | 
leży uwzględnić: 


Wyrazy długookresowe: 
[ 4. В) y sin (40 —- v) 


An В, ол d sin (40 — vj) 


(4) 1 


5.1.0 


8, ВЕ) q sin (Bw — v) 


d 
Wyrazy zwykłe: 
Зр B, o, Sin 810 
dAn Bo m sin (4w v) 
An BR W sin (4w + v) 
128% y Sin (12w—v) 


12 Bi, W sin (12w—v) 


Wyrazy o krótkim okresie: 
4 p В, оо Sin 4w 


81 BĘ.) y sin (8w— v) 


3* 


49 J. KRASBOWSKI. 


W ten sposób wyczerpaliśmy już wszystkie wyrazy w rozwi- 
nięciu funkcyi P aż do drugiego rzędu włącznie. Wyrazy trzeciego 
rzędu zwykle są już tak małe, że można ich nie uwzględniać — je- 
dnak w przypadku ruchu niektórych planet krytycznych bliskich 
Jowisza będziemy musieli uwzględnić jeszcze wyższe wyrazy 
(trzeciego rzędu) w tych rozwinięciach, a to dla tego, że one nie 
będą o wiele mniejsze od wyrazów niższych rzędów, które uwzglę- 
dniamy. Z pośród wyrazów trzeciego rzędu należy uwzględnić 
oczywista tylko te, które są zerowego stopnia. 


cos 4w; AR sin 4w; 8 Ba, COS 410 


ге 4.0.0 


B соѕ 12w; 128" 


12.0.0 we cos 12w 


adio $ 
d 


sin 12w; 62B 


12.0.6 


te wyrazy przemnożone przez Ri K dadzą nam wyrazy o długim 
io krótkim okresie. 


W rezultacie będziemy mogli ostatecznie napisać następujące 
rozwinięcie dla P. W tym rozwinięciu zostały uwzględnione 
wszystkie wyrazy szczególne aż do Ill-go rzędu, zachodzące 
w przypadku ruchu małej planety typu Thule — t. j. takiej, której 


ruch dzienny średni jest bliski 3 ruchu średniego dziennego Jo- 
wisza. 
PB Bioocos:4w + B n cosiv Bo, т соз Му 


+ BOP y cos (4w —v)+ Big) q! cos (4w—v,)+ 


4.1.0 4.0.) 
+ BD q cos (8 w —v) + By, 1! cos (w — v:)+ 


p Zo Wz ce В? qq! cos (v—v4) + B ao т? cos 4ш + 


+B,,, n° cos 4w + Ву n° cos (40 — 2v) + 
HBT? ny? cos (4w— DUCES n cos (8, — v) + 


+ ВС) w cos (80 — v) ВО) nd cos (dw+v— V,) + 
8.0.1 U 1 .1 


= —2) 
+ В? np! cos (4w — v v) + B, 


nm’ cos (4w — v — v) + 


АЁ as nn cos (ab tva qa 1 — Ú: 1m cos (1226 —v ah + 
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+ BE) n? cos (8w—2v)+ Bios 1/2 cos (8w — 2v,) + 


+ Binos 11? cos (12 — 2v,) + BG? т? cos (12w — 2v) + 


12.2.0 


+4 Bao) cos 81 4- В}? y cos (4w + v) + Bł 7 cos (4шю + v) + 


+В € cos (120 — у) + Bzy w! cos (12w—v,) У Bi 


+4 В + Bio cos 4w + BI q cos v+ BHO y! cos v + 


1.1.0 


—1.1.0 49771.1.0 y— | 
TB 108 V Ra 1 cos Vi+ B. a cos (8w — v) 4- 


+В cos (80 — v) ) (В: Ru) + 


+ (4B,,, sin 4w + 8 B,,, sin 8w + 4ВС y sin (4w—v) + 


4.1.0 


+ ART af sin (4 w — v) +4 BC y Sin (dw + v) + 


11.0 


+AB тү sin (4w + v,) +8 BC a sin (80 — v)+ 


5.1.0 


1- 8 В. m sin (8w — v)+12 BI” q sin (124 — v) + 


12.1.0 
+12815) у sin (1245 -— v) ) u (K, But) 


20 2. 
++ Bi), cos 4w+- RI 


12-0.0 


cos 12w ЎР? + 42870 


по Sin 410 + 
SUE sin 12w y y POPE { 8B, oo COS 4w+ 


ROAD OSMA ae а 


12.0.0 


Analogiczne wyrażenie otrzymamy dla wielkości Q, tylko od- 
powiednio zamieniając spółczynnik liczbowy B ze wskaźnikami na 
odpowiadający spółczynnik A. 

Otrzymamy tedy: 


P> Sin An ACC) n sin Y БАС" y sin (4w—v)+ 
Ато nsin (8w— v) + AC) 3 sin (40—у,) HAC y sin (8w—v,)-|- 
2 аи т sin (v —v,)--4,,, 9? sin 4w A n? sin (4w -- 2v)4- 
HACH 12 sin (8w — 2v) + АС т? sin (12w — 2v) 4- 


HA) m sin (404 — v ҺА nn! sin (dan v — v) 4- 
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FAC mn sin (Bu — v — vy) TAS m sin (120 — v — vy) + 
HAG? my sin (4w—v+-v,)-+ AC? +? sin (445 — 2v,) + 
А52 1? sin (81,5—2v,)H- AG, q? sin (12w — 2v,) + 

+ À 459, Sin BWAR. n sin (4wtv)t A, у sin (12w—v)+- 
+ lat T sin (4w + IRA 4! sin (12w — v) ? Ri+ | 
+A 4490 Sin 404519" «sin (4w—v)-L A; a sin (8w—v)+ 


АШ e da 9 Sin (40 — v.) + 


sr A, ew q’ sin (8w — vj) y (Ra + hy) — { ДА, оо COS 4W + 


EA o COS 80 + 4 А) y cos (41w-|-v) 3-44 C Pw! cos (4w-|-V,)-|- 
+H4A т т cos (dw — v) + 8 АС" y cos (Bw — v) + 

+ A cos (12w — +44, n cos (410 — v) + 

+845? n cos (Bw—v,)-|-12 40 w cos (120—v,) } p. Rut KAL 
+ 144, Sin 40445, Sin 12w) R2— {4410 cos 4w + 

+ 124), cos 12w ) y R, K — 4 84,,, sin 440 + 

+ 624,,,, sin 12wj KAL Á, oo COS Sue LAT y cos (4w-|-v) + 
+ AC, gege (1230 v) E ACE? yl cos (410 oa kel | 
+ Año, 9 cos (12w — vj) ). 


Abyśmy jednak mogli ostatecznie obliczyć wielkości P i Q 
musimy mieć możność obliczenia 
(5) Sa, Sk, Ra, Rk, Kr. 


Wielkości te są to znów szeregi trygonometryczne, w których 
występować będą wyrazy o któtkich lub długich okresach. Musimy 
więc obecnie, mając rozwinięcie dla P 1 (Y, wynaleźć jeszcze wartości 
dla wielkości (*), które podstawimy w nasze wzory. Przedewszy- 
stkiein należy zauważyć, że we wszystkich wielkościach (*) wystę - 
pować będą wyrazy krótko i długookresowe, zależne od tych sa- 
mych argumentów, które już występowały w rozwinięciach dla P 
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iQ. Mieć więc będziemy w rozwinieciach dla S, R, K wyrażenia 


następujące: 


Długookresowe. 
Wyraz stały 


y 
n 
nn’ cos (v — vi) 


y cos (4w — v) 


n cos (4t — v4) 


n? cos (8ш— 2v) 


Krótkookresowe. 


» COSV 
т! COS V, 
cos tw 
n cos (8w— v) 
a cos (8w .— v,) 
1? cos 4w 


т! cos (4 wHv—v;) 


(9) . т? cos (810 — 2v,) 
nn cos (8w —v—v,) 


nn’ cos (4 w—v--v,) 
n? cos 4w 


n? cos (4w— 2v) 
1" cos (4w—2v,) 
nn cos (4w—v— v4) 
1 COS (12w-—2v) 
a? cos (12w—2v,) 


i nn? cos (12w—v—v,) 


W tych wyrażeniach (9) są zawarte wszystkie kategorye argu- 
mentów Gyldén'owskich, o których wspominaliśmy na str. 15 
wzór (8) mianowicie argumenty typu 4, B, Ci D. 

Z tego, cośmy teraz powiedzieli widzimy, że do typu 4 na- 
leżeć będą 4 pierwsze wyrazy długookresowe, będą one elementar- 
ne, gdyż te wyrazy są drugiego stopnia względem n — wszystkie 
zaś pozostałe długookresowe wyrazy należeć będą do wyrażeń 
charakterystycznych kształtu C. Podobnie, do wyrazów elemen- 
tarnych krótkookresowych należeć będą dwa pierwsze wyrazy 
o krótkim okresie, reszta zaś należeć będzie do wyrazów kształ- 
tu D 

Aby móc lepiej zbadać istotę wyrazów (9) postarajmy się 
wyrazić wszystkie argumenty w postaci tej, jaką Harzer otrzy- 
mał dla każdego typu charakterystycznego i który wypisaliśmy już 
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na str. 15 (8). W tym celu przypomnijmy sobie, że mieliśmy 
związki następujące: 


3 —0 
== 4 : 
skąd: 
4(1-— u) = 1 + ô; 4w =(1-|-08,)0— 4 B— 4uW% 
Oraz 
у 0—1; 1) v,—v-—Il,, 
czyli 


v ==v(1— <) — Il; V =v(1 — s) ME 


(gdzie II,, stała, zaś Ilo znamy z teoryi Jowisza, s ruch apsyd). 


Wprowadźmy te wielkości do wzorów (9) wtedy otrzymamy: 


Const. 
kia 
(ay) q KI cos [(1 — p )v П, | 
| m’ cos [cv + (П, — П,)) In cos [(1 — 0) v- Hal 


y cos [3,7 — 4B — 4. W, + Il] 
| y cos [ 8,0 —4B-—4pW¿+1T] 
(C)! n° cos [28,0 — 8B — 8 W ¿+ 21] 
| 1? cos [22,0 — 8B -—8y4W¿-+ 211] 
nn cos [26,0 —8B — 84 Wa +III] 


cos [(1 + 9,)u —4 B -- 4. W.,] 
10); y cos [(1 +28, v —8 B — Baum, II] 
D cos [(1-H 20) — 8 B — 8uW;--Il] 
т cos [(1-- 2,) 0—4 B--4&W,| 
тт’ cos [(1+ 0) 0—48 —4y W,—11 — П, | 
E cos [I + 85v —4 B — 4 W, --II—IT,] 
(DIE i? cos [I + &)v—4B-—-4uW,]| 
y? cos [(1+ ôv + 4 B4-4pW¿— 211] 
1/? cos [(1— ô) v 4- 4 B+4pW,— 211] 
gn! cos [l — 3)v--4B-I-4&W4— II—IL] 
n cos [(1 4- 385 c— 12 B — 124 W, + 211] 
n° cos [(1 |-30,)v — 12 B— 124 W-- 211, | 
um? cos [(1 + 38,) v — 12 B— 124 W ,--II2- TI, ] 


га Brendel. 1. c., pg. 20. 
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W powyższem zestawieniu wyrazów tylko 5 nie zawiera 
wielkości 6,, jak wiemy, związanej z ruchem pewnej określonej 
planetoidy; te wyrazy odnosić się więc będą do ruchu dowolnej 
planety 1), podczas, gdy reszta wyrazów wypisanych w (10) musi 
być obliczona oddzielnie w każdym przypadku szczególnym, dla- 
tego też Gyldén nazwał te wyrazy wyrazami charakterystyczne- 
mi, jako zawierające w sobie cechy ruchu pewnej badanej plane- 
toidy. Zwróćmy się do zestawienia tych szczególnych postaci 
wyrazów podanych na str. 15, wtedy z łatwością w (10) będziemy 
mogli wydzielić wyrazy, odnoszące się do każdego rodzaju, co już 
zresztą uczyniliśmy, wydzielając klamrami każdą grupę argu- 
mentów. 

W wyrażeniach kształtu (C) i (D) występuje Wu, które jest 
bardzo mało zmienne z czasem, gdyż jest to długookresowa część 
wielkości W —redukcyi czasu. Obecność tej wielkości utrudnia 
znacznie całkowanie naszych równań różniczkowych. W wyraże- 
niach kształtu A wskutek całkowania wystąpi mały dzielnik ç któ- 
ry jest rzędu masy zakłócającej. W wyrażeniach kształtu Di C 
przez całkowanie w mianownikach wystąpi wielkość 6, która, acz- 
kolwiek mała, nie jest jednak rzędu masy zakłócającej, lecz zawie- 
ra masę z wykładnikiem mniejszym od jedności: wyrazy więc te- 
go kształtu stają się zerem wraz z masą. Zgodnie z tym, co już 
poprzednio mówiliśmy, musimy zauważyć, że gdy całkujemy rów- 
nanie dające Š. to przez całkowanie wyrazy kształtu A i C powie- 
kszają się, podczas gdy przez całkowanie równania dającego p — 
są zwiększone tylko wyrazy typu Bi D. Wskutek tej uwagi skon- 
statujemy, że Sz będzie rzędu m’, zaś S, rzędu 5 ; Ra, Ry, W; są 


/ 
rzedu o Wa i K,są tegoż samego rzędu. Na tej samej zasa- 


dzie można z łatwością stwierdzić, że zwykłe wyrazy w 8, R, K sa 
rzędu masy zakłócającej ?). 

Po tych uwagach, będziemy mogli napisać następujące wy- 
rażenia *) dla S (w S wyrażenia kształtu (A) będą conajmniej dru- 
giego stopnia). 


) Brendel. | c. pg. 106. 

2) Kramer J. Theorie der kl. Planeten von Hecuba-Typus, pg. 22. 

3) Ludendorff. Jupiterstórungen etc. kap. I u. Il sparsim.; 
Kramer J., l. c., pg. 43. 
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( S= Sa + Sk = + a, COS 4w -|- а cos v 
4-a, q cos (4W=V) ta, q cos (Bw—v) 
+ a," COS v, 
Las "| cos (4wv,) +45 w cos (80 —v,) 
+ an? cos 4w + а, т> cos (4w — 2v) 
La, cos (8w—2v) +a, n? cos (12w— 2v) 
DE 
| + aym cos (41--v — v4) 4- 9,971! cos(4$0— v — v4) 
Lo, qm cos (80—V—v;) + a ml cos (UN v4) 
ann cos (4w0—v+v,) Has? cos (4w - 2v) 
Lo 1? cos (8w—2v,)+a,, 1? cos (12w 2%;) 
+a n cos 4w + 
i aro. d 
w tem rozwinięciu jest: 
Sa=a, ayq cos (40 — у) ж? cos (80 — 2v) 
agn! cos (4w — v4) + a, 1m cos (810 — v —V,) 
-|- aet? cos (810 — 2v,) 
zaś S, Oznacza wszystkie wyrazy w S kształtu A. 
Analogicznie mieć będziemy rozwinięcie (2) dla R, zważy- 
wszy, że wyrazy kształtu (B) są już zawarte w (p) : 
(R=Db,+8, cos 4w + far cos (4w-—- v) T В, COS (80 — v) + 
-H Ba cos (4w — у) +B n° cos (880 —v) + 
| Ban? cos 4w J- B, ,?cos(4t —2v)-|-B,, 7^ cos(8&w - 2v) + 
-- 81? cos (12w— 2v) 
(12) 1 + gm! cos (4w-|- v—v,) + Pio si cos (40 — v — v) + 
4- В, т" cos (8w—v—v,) Ba тй cos (12)0—v— vy 
Lg qw. cos (tw— vv) Bt" cos (4w—2v)+ 
+B, cos (8ш —2v,) --Bis 11? cos (12w — 2v,) 


| Pio” cos 4w LR 
gdzie Ra Oznacza część E, w której sa ugrupowane wyrazy kształ: 
tu (4). 


у) Kramer J, l. с., pg. 28. 
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Wyrazy długookresowe w R czyli Rz są: 

Ra= bo LB cos (40 —v )+B,n? cos (8w — 2v) 
3- B, »' cos (dw — v,) + Bis тр cos (8w — v — v4) 
+B, 2 cos (8w —2v,). 

Jak to już poprzednio zaznaczyliśmy, po zcałkowaniu w tych 
wyrazach nie będzie występować mały dzielnik 8. 

Ponieważ K¿=0, przeto wielkość K = K;4-K,. Niektóre wy- 
razy zwyczajne będziemy musieli uwzględnić, gdyż one w naszym 


rozwinięciu będą już nawet pierwszego rzędu. Mieć więc bę- 
dziemy: 


( K=x, sin 4w + y, 4 Sin (8 w— v) 7, 1 sin (4w v) 
+7, лү sin (8 w— v) 
+7, 4? sin 4w +" "sin (4w — 2v) 


+z 1 sin (12w — 2v) 
- Y m’ sin (4w+v—v,) +7, m sin (4w — v — ү) 
Hp sin (12:0 — v — vi) 
ts m sin (duty) Hr, sin (du — 2v) 
+ 9 sin (12w — 2v,) 
-l- Yio Oe sin 4w + Yao 1° sin (dw + 2v) 
Ha 9 sin 8w 


F Yaa sin (8W-|- v — vi) 
gdzie jak widać: 


(13): 


Ken sin (4w — v) 

+ tao” sin (4w + 2v) Yan? sin 8w Lt qq sin (80 + v — v). 
dw dw) 
de ZIE (del: 

Dla tej samej wielkości mieć będziemy rozwinięcie nastę - 
pujące: 


Wreszcie część o długim okresie wyrażenia 


ale 1,1 cos (4w — v) + 1,41? cos (8ш — 2v) 


(14) | + Тз 1 cos (40 —vi) + Yis m’ cos (8w—v — v) 


“!2 an ae dw 
+ 116 1 cos (8w Seel 


O ruchu planetold. 
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Е lw 
gdzie znów przez Ka oznaczymy wyrazy kształtu (4). 
A 
Należy zauważyć, że spółczynniki a są rzędu masy zakłóca- 
bn b BAG ZE: m En 
jącej, zaś spölczynniki a i B są rzędu > (masa zakłócająca po- 
dzielona przez mały dzielnik). 


Z wyrażeń (12) (13) (14) odrazu widzimy, które rodzaje wy- 
razów zachodzą w każdej z uważanych wielkości. Napisane szere- 
gi musimy wstawić w wyrażenia wyżej otrzymane dla P i Q i wy- 
konać mnożenie, grupując razem wyrazy zależne od jednego argu- 
mentu. Przy mnożeniu musimy uważać aby w rezultacie, zgodnie 
z zasadniczą myślą Gyldćna uwzględniać tylko wyrazy o dłu- 
gim lub krótkim okresie, zaś opuszczać wszystkie zwykłe wyrazy, 
które później w otrzymanym rozwinięciu na P 1 (Y sprowadzimy, 
przez dodanie sumarycznego symbolu. Jednakże, ponieważ w ilo- 
czynie zawierającym K, wystąpią, jako rezultat mnożenia wyrazów 
zwykłych w K przez wyrazy krótkookresowe w nawiasie, nietylko 
wyrazy zwykłe, które zwykle pomijamy, lecz także wyrazy krótko- 
okresowe, więc musieliśmy także te wyrazy zwykłe uwzględnić 
w rozwinięciu K. 

Wykonajmy więc wskazane podstawienie i mnożenie, korzy- 
stając z elementarnych związków goniometrycznych dwuch kątów 
À ig 

l 


| 
sin À . sin p= — y cos Mrs cos (A —- u.) 


re ; 
(а) sin À . cos p = 9 Sin (À +) - 5 sin (À — u.) 


1 
COS À . COS y = > cos (A+ u.) + y COS (A— u). 


Otrzymujemy wzór nastepujacy: 
P= Bono + B400 COS 4w + Ba cos v+ B; 1, cos (4w—v) + 
+ Ве! cos (Bw — v) + В} үлү cos v, + Bio, 11 cos (4w —v,) + 
+ RS gd cos (Bw-—v,)+ Bazo +B, 9,0 + EN y am cos (v— v) + 
F B a? cos 404 ВУ чп cos Awt v—v,)+ 
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+ Bree 1? cos (4w—2v) + DE. 12 cos (8 w— 2v) + 

+ Lé, cos (12w— 2v) + B. qq cos (4w — v + vi) +- 

+ Ba, cos (40 — v — v) + Ву m’ cos (Bin —v — v,) -+ 
RE m m cos (12w—v- v,) + Bios N” cos 4w + 

+ Din. 1'* cos (4w—2v,)+B,,, n? cos (80 — 2v,) + 

+ Bei" сөз 120 Sie) AL 

TA, cos 81 + BIE cos (4w--v) + By; qcos(129w у) 
+ Bio W cos (4w+-v,) 4- Bos тү cos (12w — vj) ) 

1 B, cos 4w--B, y cos (810 —v) +B, n° cos 4w4- Bi, n° cos (4w—2v)- 
+ Ёла 12 cos (12)0—2v)-- 8, 1! cos (8w—v,)-E8, jn cos (4w--v-- v) 
+ Bra! cos (4w — v — v) + Dem cos (12w — v — v4) + 
+ Bann’ cos (dw — v -+ vi) + Bis 42 cos (4w — 2v,) + 

+ fie"? cos (12w — 2v) + Bor? cos 400 ) 

+4 Bic + B оо cos 4w + Bus 1 cos v + 

+ Bio 1 cos (40—v)-- BZ" 4 cos (Bi— v) + 

+ By чү cos vi HB q! cos (4w—v,)+ BZ 1’ cos (8w—v,)) 
À ba + $, cos 4w + 85,1 cos (dän v) -+ 5, 1 cos (8w—v) + 

+ Bs cos (4w—v,) + 6, 1' cos (Bw — v.) -+ f»? cos 4w -|- 

- Bui? cos (4w—2v)-|-f,,4? cos (810—2v) +64, 4? cos (120—2v)-- 
+ Bam cos (40 + v — v.) + Bau! COS (40 — у —v) + 

+ Bis тй cos (Bio — v — v4) + Bra ria cos (1990 — v — vp + 

+ Ps mj cos (dw — vv) + 5,312 cos (5t — 2v) + 

+ Se" cos (8w— 2v,) + 6,, 4/2 cos (12w— 2vi) Bal? cos 4w | 


+- t. T B, ou sin 4w + 8B, ,, Sin Su + 4 D T| Sin (4) + v)4- 
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AR: 


4.0.1 


ad sin (4w + v) + 4B,,, 4 sin (4w—v) + 
+8B,,, 4 sin (8w—v) + 12 By, „ sin (120 — v) + 
-- 4B,, 1 sin (4w — v,)+8B;, q' sin (8w—v,) + 
+ 12B;,,, «/ sin (12w—v;)} 
| y, sin 4w-- 1, ү sin (80 — v) + y, 4 sin (4w-|-V) + т, -+ 
+75 1 sin (8w — v,) + 177% sin 4w + т 12 sin (4w — 2v) -|- 
+ yir sin (12w - 2v) + Ysi sin (40 + v — v4) + 
y, M sin (40 — v — v4) - Tig t! Sin (120 — v + v4) + 
+ yan! sin (4w— v +) + 7111? sin (40 — 2v,) + 
+ 1,9! sin (120 — 2v4) + 110 1/* sin 4w -|- 
+7,07? sin (4w-|- 2v) + ta Y sin 810 -+ Tas tn! sin (8802 + v — v,) Í 
+ | Bij, cos 4w + В, cos 12014 (81? +3 B? cos Sw! 
+ [1 4B sin 4w + 12 Bizo sin 12w | 
148,7, sin 8w} — p {8В, оо cos 47-62 B12.0.0 cos 12w, 
142—415? cos 8wj]y.. | 
Z powyższego wyrażenia otrzymujemy ostatecznie: 


(15) Р=р, + р, cos 4w—+ pin cos v p,q cos (4w — v )-|- 
+ Den cos (80 — v) + po? cos (12w—2v) 
+ pz cos v, + p; 1 cos (40 — v.) -p:n cos (Bw — vi) + 
+ Pay M cos (12w — v — vi) 


TP 7? 4-01, N° cos (4w—2v)+- 
+2,79? cos (8w--2v)1-pa, N? cos (12w—2v,) 
Fpa q cos 4w Dis m COS (4Ww—V—V,) + 


Pie vm. COS (Bw—v--v,) 4-03 W COS (v— vi) 
À Dio "n cos (4tw--v— vi) +-p,, n°? cos (4w—-2V,)+ 
—-p,, N? cos (8w—2v;) 
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Fu «m cos (4w—v-Fvi) 
Pis 9° cos 4w 
+Pıs WE dub 


Przez Œ, oznaczyliśmy sumę wyrazów zwykłych. 
Dla krótkości oznaczyliśmy: 
p= B, +; Bio B, + 24 B, ee 2 Bien pl 
py— B, „+ | Boo + $ В|, tB... F 
аш КЕ сай н e ER 
u. Bio З Biz 1, nl 22,09 Bizoo) 
PB + Bach + Bo Bi +4 Boo Bah Bo Bd 
Län Ban НВ YZ Gu SC, hb 2€ en hi 
DB Hi Be HABI B + 4.81008. + ÄR h + 
Lau, t + Że Born + 28 Buts 
E a BE eH S ci EE Baer f 
21. RB pilak ŻE Beefert: Á t. BA 
рь =B Bono Ps + 4 Boo Bs--4 Bano bs +3 Bios Ва. f 
TB, EE Ré 
Bee et Boo + Broo ba + Во + I Bano Beck Bann T 
eco Баб, 
Be Bad FAB бї Baas + à Boo Bs + Bo f 
+61 Bst — Że Bout 
Poeci d Boo Pr KAPA, Hai Be B Ж Boy T: 
+4p В, „та + Be + 41 Bao Y 
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+ AGS, sin (12w—2v,) + | | 


+ {410 sin 8w+ А0119 qsin(4w+v)+ Ai, q sin (120—v) + 
+ Age W sin (12% — vi) }. 
| 3, cos 4w -- 6," cos (80 — v) + 8, y! cos (810 — v,) + 
+ 8,4? cos 4w В,, 3? cos (4w — 2v) + E, , 1? cos (12w — 2v) + 
+ By cos (4w + v — v H- Bann cos (4Ww—v — v4) + 
4-8,6 тр’ cos (12w — v — v4) + әлү cos (dw —v t vi) + | 
+B, 1? cos (4w— 2v,) 4- 8, cos (12w — 2v,) Bue" cos 4у{ + 
| 
| 
| 
| 


A 2а sin TE Azją" y sin (4w — v) + Арт sin (80 — v) + 


1.1.0 


FAT Ww sin v, + Ашу sin (40—v,)+ 

-„А ут sin (Bw — vj) j 
| ba +8, cos (4w—v) + B, 2 cos (8w—2v)+ Ban cos (4tw—v,)-|- 
-H fs um cos (8ш — v — v,) + Bie? cos (8: — 2v,)-|- 
LB, nn cos (4w — v — vi) + Bra cos (12w — v — V,) +- 
+ 8, ml cos (4w — v + v) 6459 * cos (4w — 2v,) + 
+ 4,1? cos (12w — 2v) + 101? cos 4w + B, cos 4w + 
+B, cos (Bw — у) 4-6, 1 cos (Bw — v,) +8, 1? cos 4w + 
+ Burn" cos (4w-- v — vj) + By, q cos (40 — v — v) + 
+ f» cos (40 — 2%) | — 
— {44,99 COS 40 + 84,4, COS 80 + 4А n cos (4w-v) + 
+44 w cos (40-Е) +441, y cos (40 — v) + 


+84! m cos (80 — у) +12 А5) n cos (12w —v)+ 


5.1.0 
+ 44,5 y cos (4w — v) + 8 AE n’ cos (dw — v) +- 
+12 Anaon cos (1290 — vi) pp 
| x, sin 4w+ y, y sin (80 — v) + 7,74 sin (4w + v) + 


y +17 sin (8w — у) +14? sin 440 +T sin (4w — 2v) + 
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11797 Sin (12w — 2v) +7, тр’ sin (4406 4- v — v) + 

+ 112701 sin (410 — v —vj)-E m sin (120 — v+ v,) + 

+ Т9 117 sin (4w— v - vj) J- 3,47? sin (4w — 2v,) + 

06190 Sin (12w —2v,)+ 3,41? sin 4w- тт sin (4w+- 2v) + 
Tar" sin 80 + yaa sin (810 - v — v,)} + 

+4 42° sin Ane-LA e sin 12w) (382 +4 Bl, cos 8w 3 — 


— 4440, cos 4w + 124% cos 12w } n (18,0, sin 8w Y — 


12.0.0 


— À 8 A, оо Sin 40462 Aan sin 12w j Ay? 11,? cos 8w | 
laan COS 8w + AŻ т cos (4w--v) + Aja, cos (12w —v)+ 
FAR, ү’ cos (4w + v HAGA тү cos (12w —v,) . 


12.0.1 


Te wyrażenia uporządkujemy, uwzględniając związek gonio- 
metryczny (a) w ten sposób, abyśmy otrzymali wyrazy jednako- 
wych argumentów dla wyrazów charakterystycznych oraz elemen- 
tarnych. Zachować będziemy musieli także i niektóre wyrazy 
zwykłe, zależne od argumentu 2v, gdyż podobne wyrazy wystąpią 
póżniej w wyrażeniach dla W. 

Gdy wykonamy wskazane podstawienie oraz przekształcenie, 
otrzymamy wtedy ostatecznie wzór następujący: 


| Q=q, sin 4w + 9,7 sin v + q, + sin (4w — v) 
+ 61 sin (Bu — v) 
+ Qs Sin v, -Hq54' sin (4w — v,) 
qua! sin (8w — v.) 
+ 95 12 sin 4w + 9447? sin (410 —2v)4-q, у? sin (8w—2v) 
+9,89? sin (12:10 —2v) 
(16) + gm sin (Aw v —v,) + q, mm sin (4w — v — v, + 
+ дү! sin (8w—v—v,)+,9 mn sin (12:0 —v—v,) 
+ on Sin (4w—v--v,) + Qu 4? sin (446 —2v,) + 
+ 9,9” sin (Bio — 2v,) + 95," sin (12w —- 2v ,) 


дул? Sin 4w + Gm! Sin (v — v,) 
g, sin8w +g, y sin (4w+v) --g," sin (12w—v) 
| + 9,4 sin (40--V,) + g, W! sin (12w- v,)+ Go. 
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Przez q oznaczyliśmy spółczynniki wyrazów charakterystycz- 
nych i elementarnych, zaś g oznacza spółczynniki największych 
wyrazów zwykłych. 


Wartości dla poszczególnych spółczynników są następujące : 
Ф.А, 450 Py + 4 Ast Awer 14400 ЕРА вов" 
—plAgoo Зоо ! 8114-82 | —4 44,59 2400 7 Aro ry 
о AR By + Agoo A io Pr — 284300 s E 
"Ab dza l ZA Ws Ze, T 
ФА ы А Ды B + 008 + Arno — E Fr 
+44 AgopTs + 26 Ajos li — 21430, 
MEA ¡o— $ Алое вА В, — 2p Алоо, + 4t AsooTe + 
+ Ap Az от 
Aa — 4 Algo Po — 44001 Bb A Bi — 28 Ago a + 
sn Aw Aga, Ti 
geed 4-44 8, +3 Aeon Ba А Bıt O Arno Ta => 
— 4n Aio 
E 1 D АВ, TE А 86 401 8 FOLA hi — 
— 2p AA a 
Aa + A Ba A An Bad Ao Pa — 
— 2p Anota + Ag aos; + 60 Ajo о — 29 Ао 
A 4-4 ASSI A B А 6,4-4 Аа" Ba — 
— 24 В, 00 122 F 48 Aaoots бе Aso. T4 — "E 
й»=А ы bdzie Fa tł oo Po Fän Ba — Aion 8: А 


LEAST Ba F 4p Asoo Ye "r GAS IE — 24 AY s. T4 
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014, ort b det Ps +4 Або 1-4 А BAR, + 
+4 Año Ber A 16 б ASA ys au SCH ai 05 
due aa à AS š Saar Ae 0.0 бм -4 4450 Pia 4p Aggo in + 
+ 4t А „Тв — 2 Ад at 4 А т, 
Qu Aqq, — MA BA B Ф410 Be — 4 Ara, B — 
— HA B,— 4 p Aano Yis — ŻE Aura = 2p A, Lota dp. Az o Te 
WA Арол Bs — Ам ә — HA 8, — ARO, — 
Audio 20 Año, Ts 
Aso Ato ®А Bit + 44,5 8 Le A „ти — 
Ahr б А; „To 
MAPA o Bia 144908 А BA BH 
HA aic eg nr 2p A oo 118 op SE ues Ts 
q= A АД Bis — 4.4490 бо Ао" Bs Ed Ago] Ba — 
AO A sro A bil» 
Ф455, + ł Ash Ао ti Ano Ba #45 Be — 
4А от: 28 B 0 T4 
Фо = Ant Akta Bo Ast &4 р 6 3 Boy B + 
А ef Zë = Mesue B,— 4р4, oo? — 2t. Zä la. le Ys 


Фо = А0: 346 Bis оо oi Bst da Ps — 
— 4p Agoo Yis — РТУ Az, Ts 
da = Aora Axao Be ААД Bo Ао Pa А-6 Aa Br — 
A BF Aes BO 264, Ta 2U ae 
d'Ee Tes" dg + 4p dts 
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Jı = Ао} Ai oo B + Aa 26 60 17 68 Aan ft 
ga == A; o t b 400 Be FALLO BA В o A брот, 
Оз = Aqo t t Aooa Pr + à Ages + 4.420) 6, +4 42,5; + 
oe LU Año, Tam OU ANTE 

J= Anot Aach, +4 Agra ñ 4-5 45008284, oot 

w A eh SPA la r SU eg EE E 
gs = Azoti 40085 + YA, Bi +4 Аво оз 20A, oo Ts 

— 4 u Años ti + 88 Ajos т, + E aen 


IV. Całkowanie równań różniczkowych, 
odnoszących się do ruchu planetoidy (279) Thule. 


W rozdziale I, równ. IV, V, VI dają nam wyrażenia różnicz- 
kowe spółrzędnych Gyldćnowskich; w poprzednim rozdziale 
zdołaliśmy już oddzielić wyrazy typowe odnoszące się do ruchu pla- 
netoid. Brendel!) i Kramer podzielili planetoidy wedle cech 
ich ruchu na pewne klasy. Według tego podziału planetoidę Thu- 
le (279) zalicza się do planetoid pierwszej klasy. W ruchu tej kla- 
sy planetoid wyrazy stopnia 0 w Ri K są znaczne w porównaniu 
z masą Jowisza. W tej klasie planetoid wystąpi jeden bardzo znacz- 
ny wyraz z w rozwinięciu funkcyi À, i ten wyraz, jak to wykazał 
Brendel?) będzie kształtu D, według oznaczeń Harzera czyli 
charakterystycznym krótkookresowym. Do tej samej 
klasy co Thule należy także i planetoida Hekuba, której ruch, jak 
wiemy, był niejednokrotnie badany przy pomocy metod Gylde- 
nowskich?), lub ich modyfikacyi. Zasadnicza jednak różnica po- 
między temi dwiema planetoidami polega na tem, że Thule przybli- 
żając się znacznie do drogi Jowisza, z jednej strony ulega silnym 
perturbacyoim, z drugiej zaś dla tej planetoidy wskutek znacznej 


1 Brendel Th. d. kl. PL I, pg. 92; Kramer. Th, d. kl. PL, pg. 20, 
str. 199. 

° Brendel. Om amvandningen, str. 26. 

3) Patrz prace Harzer'a, Kramera, Backlund'a. 
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blizkości drogi Jowisza i drogi Thule rozwinięcia wedle potęg sto- 
sunku osi lub promieni wodzących są bardzo powoli zbieżne. Aby 
otrzymać wystarczająco zbieżne rozwinięcia perturbacyi Thule 
(tak, aby położenie planetoidy było określone z dokładnością oko- 
ło 1' na przeciąg jednego stulecia) możemy w rozwinięciach spół- 
rzędnych i w równ. IV, V, IV wyżej wymienionego rozdziału —- po- 
minąć wyrażenia rzędu trzeciej potęgi masy Jowisza, zachować 
zaś w wyrażeniach O stopnia wielkości rzędu trzeciej potęgi masy 
Jowisza, podzielonej przez mały czynnik <. Natomiast w wyraże- 
niach I-go i Il-go stopnia wystarczy uwzględnić wyrażenia rzędu 
kwadratu masy zakłócającej, podzielonej przez mały czynnik ç, 
a to dlatego, że w tych wyrazach występuje jeszcze bardzo mały 
czynnik y lub y! w pierwszej i drugiej potędze. W rozdziale VII 
pierwszej części swojej „Theorie der kleinen Planeten* Brendel 
М Ew uet E dr . ds 

podat równania rózniczkowe, wyrazajace do? mu 

Ponieważ w zadaniu naszem postanowiliśmy się ograniczyć 
do wielkości rzędu do trzeciej potęgi masy ciała zakłócającego, 
dlatego też musimy rozwinięcie Brendel'a jeszcze dopełnić 
w sposób analogiczny, jak to zrobił Kramer dla Hekuby ^. Gdy 
dodamy do wzorów Brendel'a, odnoszących się do ruchów He- 
kuby, wyrazy dopełniające, otrzymamy wtedy, z wystarczającą zu- 
pełnie dla naszego celu dokładnością, wzory różniczkowe na związ- 
ki pomiędzy spółrzędnemi Gyldénowskiemi dla planetoidy 
I klasy. 

Ogólnie każdą funkcyę zależną od v to znaczy Q, P, R, S, 
K, Wit. p. napiszemy w ten sposób, że wydzielimy w oddzielną 
grupę wyrazy 0, 1, 2 it.d. stopnia i każdą grupę opatrzymy 
wskaźnikiem, który oznaczać nam będzie stopień uważanego wy- 


razu, np. 
= Q, +0 HO, H 


znaczy, że funkcya () składa się z wyrazów stopnia 0, 1, 2it d. 
Na stopień i rząd dwu wyrazów musimy zwracać baczną uwagę 
przy mnożeniu odnośnych wielkości, aby w wartościach P, Q, R... 


nie otrzymać błędu większego od s Zajmijmy sie wyznacze- 


niem tych wyrazów dopetniajacych. 


') Kramer. Th. d. kl. Planeten, I, pg. 44, 46 sqq. sparsim; Bren- 
del, l c. pg. 139. 


62 J. KRABSOWSKI. 


Mieliśmy 
p=(p) + Ë, 
gdzie (p) zawiera w sobie wszystkie wyrazy elementarne kształtu D, 


zaś w li są zgrupowane wszystkie wyrazy elementarne kształtu 4, 
oraz charakterystyczne kształtu Ci D. Stąd otrzymamy: 


dy_d(p) , dk 
dv dv ' dv 
zm możemy napisać w postaci 
iR — (2) š Ë , [A Е\ 
dv ` del" du ЖЕТ, 
gdzie wskaźniki О, 1, 2 oznaczają stopień odnośnych części w roz- 
łożeniu 22 
dv 


Analogicznie będziemy mogli napisać: 
Q=Q,+Q, +Q i S=% +81 +83. 


W rozwinięciu R występują spółczynniki ß, które, jak wiemy, 

m! ' | : LR 

są rzędu ^ (pg. 50) a więc wyrazy zależne od (6, +4, + Qa) ` x 
0 

wejdą w nasze rozwinięcie. Wskutek różniczkowania wyrażeń 


długookresowych R, i R,, gdy tworzyć będziemy wielkości E 


1 

„di Р n? 

l ży, spólczynniki ß w tych wyrazach będą pomnożone przez mały 
2 


czynnik ô, przez co te wyrażenia staną sie rzędu istotnie masy 
zakłócającej. Gdy zaś otrzymane wyrażenia zostaną pomnożone 
przez funkcye Q, która jest rzędu masy zakłócającej, to w rezulta- 
cie otrzymamy wielkości rzędu kwadratu masy zakłócającej, które 
będziemy musieli jeszcze uwzględnić. Zatem mieć będziemy: 


(0-0 (5) 


gdzie wskaźnik p oznacza, że w uwzględnimy tylko niektóre | 


wyrazy, mianowicie te, które są kształtu C. Podobnie przybędzie 
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1 Р 
jeszcze i wyraz — Q, (62). Ponieważ S jest rzędu m’, to iloczy- 
ny następujące będą musiały być uwzględnione: 
jdR 
2 S, Q, rk —28.P, i —28,P,; —28,8,40, 
+28,49Q,nsinv, + (S, 9%, 


gdzie przez wskaźnik (A, C) zaznaczam, że tylko tego kształtu wy- 
razy w odnośnem rozwinięciu są uwzględnione. Ponieważ dalej 


dn? . 
zl jest stopnia 2-go ale rzędu masy zakłócającej, przeto będziemy 
E ai "s ата R 
musieli uwzględnić jeszcze wyrazy — 2 | | będące jeszcze 
dv Ado lo 


drugiego stopnia. W ten sposób w wyrażeniu rózniczkowem da- 
jącem p jako funkcyę v zostaną uwzględnione wszystkie wyrazy, 
aż do rzędu trzeciej potęgi masy zakłócającej włącznie. 

Gdy do wyrazów podanych przez Brendel'a (!. c. kap. VII) 
dodamy jeszcze wyżej wymienione wyrazy (porówn. Kramer I, 
pg. 42 et sqq.), to otrzymamy poszukiwane rozwinięcie postaci 
następującej: 


dh 


(lv 


( d? p 


dpi 1 9 7 49 Y P,— Q| do), + 


| = 25 Q (9 


+28, — P,— 284 0 P, 128,840 + 


ei ~ dh 
1p 


— si LI 
+ Qan sin v des A qp „St 


+28, — P, — 28,40 P, — 28,40 P, + 


+28, (8,409 + 8,49 )?-- Q, y sin v 


p dR ¡AR | ¡ak 
dea Qi | dvi, Qo \ ag XM 


-— [bGO.—————— M— MÀ 


+ 28,40 Q,n sin v— 2 


dn = 
dv E 


dv 
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Ograniczając się do takiego samego stopnia dokładności, 
otrzymamy analogicznie drugie równanie następujące: ') 


45 d 1 dy? 
(9) 42 =— 0.— Q — 38,00, — 38,9 0—38, 0,— 5 фу: 


Równ. (1) i (2) są właśnie równaniami szukanemi. Zajmie- 
my się obecnie całkowaniem. Zanim jednak to zrobimy, musimy 
równ. (1) i (2) jeszcze cokolwiek przekształcić. Wielkości Рі Q 
składają się z szeregów, w których występują spółczynniki nie- 
oznaczone. Np. w R mieliśmy różne spółczynniki B, w drugich 
zaś funkcyach spólczynnikl o iy. Aby sobie ułatwić zadanie, po- 
staramy zgodnie z myślą przewodnią Brendel'a*) wyrazić wszy- 
stkie spółczynniki a i y przez spółczynniki B; w ten sposób uzy- 
skamy wielkie uproszczenie naszych wzorów. Ponieważ w tej 
pracy ograniczamy się do rozwinięć liczbowych nie idących dalej, 
niż drugi stopień (Thule posiada bardzo mały mimośród — przeto 
w pierwszym przybliżeniu; pomijając wyrazy wyższych rzędów, 
niż pierwszy względem mimośrodu, nie popełniamy wielkiej nie- 
dokładności), możemy napisać następujące przybliżone wyrażenia, 
które zawierać będą wyrazy do l-go stopnia włącznie, mianowicie: 


(3) S= а,1 cos (4w — v)+a,n cos (4w —v,) 
(w tej funkcyi S są znaczne wyrazy typu С) 


oraz 


(4) R—], cos 404-8, y cos (4w — v)+ В; ү cos (4t — v, + 
+B, x cos (810 — v) + B, »' cos (8 w — v). 
(Widzimy, że w R występują wyrazy charakterystyczne typu 
Ci D). Ale z (1) uwzględniając definicyę dla R, otrzymamy: 
d? R 
dv? + R=25.. 


a ponieważ 
d?R 
dv? 


1) Ludendorff, 1. c., Kap. II sparsim.; Kramer, |. c, pg. 43, 51. 
2 Brendel. Th. d. kl Planeten I. p. 134; Ludendorff, I. c., pg. 8; 
Kramer id. p. 24; Brendel: Om amvändningen, pg. 29. 
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jest co najwyżej rzędu iloczynu masy zakłócającej przez mały czyn- 
nik, więc możemy je opuścić, wtedy otrzymamy 


Tess 2 S, 
czyli 
(5) R = 20) cos (4w — v) 920, y cos (4w — v). 
Uwzględniając (4) i (5) jednocześnie otrzymamy z łatwością !) 


PF | 
mon. B, = 2 Ps: 


Ag 


Przystąpmy teraz do wyznaczenia spółczynników 


їз MUS JĄ AD) 


Mieliśmy na str. 20 wyrażenie VI, dające nam o w tun- 


kcyi Ri S mianowicie: 
d W 


(6) dp =—2R,+8,=2R, 4-6 Ry cos v+-... 


Z drugiej zaś strony mieliśmy na str. 37: 
W—y0+ Wi+K, 
gdzie yv oznacza wyrazy wiekowe. 
Po zróżniczkowaniu otrzymujemy: 
dw — (Wa, dk 
(a) CN dv 
ale 
(b K=y, sin 4w + т, т sin (80 — v) + 7,1 Sin (4w + v) 
+ 7% 7 sin (846 — y), 


gdzie 
йз, По E 
oraz część długookresową 
aw | 
(C) dy Ub 4 COS (41£ — V) + ү; 1) cos (40—V)--4,... 


') Brendel. Th. d. kl. Planeten, p. 134. 
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Wyznaczymy te spółczynniki wyrazów z dokładnością więk- 
szą, niż a i ß, mianowicie będziemy musieli pominąć dopiero 
wyrażenia 3-go rzędu, a z wyrażeń 2-go rzędu te, które Są istot- 
nie rzędu pierwszego lub wyższego, przeto wyrażenie (6) może- 
my także napisać w postaci następującej: | 


dw | 
dv y—2f, cos 4w + (а, — 25 + 38,8, in cos (4w—v)+ 


+ (х, — 284 + 38,6.) cos (4w—v,) — 28,9 cos (80 — v) -- 
— 26,1 cos (8w—v,) + 3B, cos (4w + v). 


(7) 


Z drugiej zas strony uwzgledniajac (a), (b) 1 (c) oraz (10) 
otrzymamy: 


dW 
(8) dy = 1 + (+81) 1, Cos 4w +y cos (4w— v) + 


+ (1+20,) 7,97 cos (Bw —v 4-1, cos (4w—v,) + 
+ (1+20,) y; n cos (80 —v +(2+8,) Ys n cos (4w+ v). 
Z (7) i (8) wynika, że: 
ka = +) cos 4w -- 7, 1 COS (4w—v)-+ 
+ (1-28) т cos (Bw — v) + 14 cos (4w — v) + 
J- (1 4- 204) Ys q cos (8w — v) +-(2+-0,) т; т cos (4w + v) = 
(9) = —2f, cos 4w -+ (а, — 28, + 38,8.) cos (4w — v) + 


+ (0, — 2534-366.) 4 cos (4w — v) — 2f, y cos (8w — v) — 
— 28. 1 cos (810 — v) 4- 3f, n cos (4w + v). 


Z (9) z łatwością otrzymujemy związki następujące !): 
2 
D =: п=38 HR +R; 


2 | 
Ts — he T 33185; = — rds › | | 
1 
DEE ЖН ZB, | y^ 9 B, 
A ALO E 


H) Brendel l, pg. 142; Ludendorff, l. c. š 4. 
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Gdybyśmy w wyrażeniu +, ograniczyli się do takiego samego 
stopnia dokładności, z jaką wyznaczyliśmy a,, to otrzymalibyś- 
my wzory następujące : 


(11) A dk: Wie Bas 
GRE Palac ЧЕТЕ аву: с 
Тв ="/2f,. 
Na str. 57 mieliśmy wzór dla Q, w którym występowały 
spółczynniki o, Q, ...; napiszmy ten wzór, ograniczając się do 


wyrazów pierwszego rzędu tylko, mieć wtedy będziemy: 

(12) (Q — q, sin 4w-- 0,7 sin v Lon sin (4w — v) + 
+ Q 1 Sin (Bw — v) 4- Q, q Sin v, +s 7/ sin (4w —v,) + 
+9 1" sin (8w — у). 


Ale na zasadzie związków, wypisanych na str. 58, 59, widzi- 
my, że wszystkie spółczynniki q sa także funkcyami spółczynni- 
ków 7 i B. 

Zamieńmy w odnośnych spółczynnikach wszystkie y i B, po- 
sługując się związkami (10) i (11) 1 uporządkujmy rezultat tak, aby 


spółczynnik ß występował jako czynnik, wtedy otrzymamy związki 
następujące: 


93 —QUP TOR P > 
Q = 92" +009 B, + 9,9 8, + 9,9 B, 
qa = q |) + 9309 2, 9,09 B, + 9408, 
9, — 940) + 940! &, Long, 
4; — 50? -F 9,0 В, + 9,0 6, 
d, 9:9 A + 9,09 B, 
Q. —= Qu 9:08, 4- 906, 


Analogicznie mieliśmy także w owem wyrażeniu dla () część 
wyrazów zwykłych w postaci następującej: 


Q= 9, sin 8w-+ 9,7 sin (4w -|- v) + 9.1 sin (12w — v) 4- 
9:4 Sin (41 $ v,) + 9, 1 sin (1210 — Vi) 
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będziemy więc mogli napisać dla spółczynników g wyrażenia na- 
stępujące : 
9, =g,0--9, 965 
| g, 05®) NM, 
Zach E 
| g = g) SÉ 8. 9.26.49. 8, 
PEI, Ер BET GER 
Z (13 i (14) oraz ze związków na str. 58 159 widzimy, że 


spółczynniki q” sa funkcyami wielkości Ai p, a więc sa to wiel- 
kości liczbowe, które będziemy mogli bez trudności obliczyć. 


(14) 


Zatem mieć będziemy: 
(15) Te orki ZES) Pa 
ne AA FA A) 
— d? A ad —84,,, + E: EM 
I i MA éd А0081 ДЫ ал Дш 08 


02000 = А; QU — d А5 — S 4,, 
049 == Azor ; Фо — AŻ: — 86 ( Años — 45) | 
q |) = 5 а dt Asoo — 8k A oo 


OA 0; ПЕ dodi o U AA; 
qi = i Ao EC 


gz E= Frs — Aha — 16m 455, jr 


8.0.1 0.0.1 
9,9 = — 4.440, + Aso 
O= А1; QW "KA T А, I АЗА у 247): 
q,U == becht | 


ALA q E E Eet ata ik dy. (ЗА Ek 


= PA : 


4.0.0 ? 
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Oraz: 


(16) Y” = Agua; дү“ ) = (4, ula 2 o) zi SUC Me SAN 00); 
gue 4 i tg == p АД — 8, AN, O A 


8.1.0 8.0.0? 
QUEM Ayo 12 А „е; 
gy "d (А ыл. AR AA SES Kan g |= Ao; 
g |) =} A — 1284,00; 
g UMP g, =} (A "yd GE AR T2» (Ав s 
+3 4600 — BAR To ): 
g, 9 = 4 dg, eo 94” —4 че TH Ayo): 
EA gu (AR ds ANA 
+8 р. CARA — DA; 


g, 9 == (А, + 84 4,00): 


Zupełnie analogicznie postępując będziemy mogli napisać 
funkcyę P w postaci następującej: 


P= p, + p, cos 4w-|- p,» cos v + p,n cos (4w — v) + 
+ 1] cos (80 + v) -I- py Y cos v, + py"! cos (4v — v) + 
4-р; cos (810 — v) 
gdzie spółczynniki p mają wartości następujące: 
| Pa P: Po" BPO BA? 
D, —p, OHP 5,2 - p, Up? 
р;==р;'?--р,\8,-Ер„»®8,-Ер„®, 
! Deeg BP, 
| 


p =p, -- p OB, J-p,C 8, ]-p, 18, 
P=: p V Hp 8, -p; GG. 
a ta a LS e 8, 
PZP” Mell A O, 0, 


10 
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| s^ 1:0 : 0 
(18) p,” —— 129и DAS Š u Bio ba Au ZSZ Po d = j RS a 


р. =B; o0 Dy = Весо 2 Bean — 84 Bzy; 
py”) = Š Ве" Ed Ed B Bb) ) (В, 3B 00) + 


1-21? (31 Biamo— A 
рз Bios 1р0 A (RT + Bue?) — (4B, 
qu 4В 7 y == 3з RUE 
п M zy et. wą” = 1 Bi — 8B, „4: 


pU = By; pU (В -EBM-- (BR, + 
"I RE ZPB Py 5) = 4 Boy — SL BN 
Pa VU DER UA NT =F Б) )— 2 (4 В, š 


ЗВ, alt P= $ Bas — 41. В 


400 4.0.0 ? 


DUB; p военна АШУ ЫСЫ 
DOTE Воо В ру = yB - dAn D, oo; 
p, = ELA A prag Find Boo: r zy eh SC Bio 
nisi de NO ię ZE 
p, = By 019} @ E ak Bi d) 4р. (3 Big Г 

a AB JO PNA IE CS 


W ten sposób mamy więc już w granicach przybliżenia, któ- 


re zamierzaliśmy osiągnąć, wyznaczone wszystkie spółczynniki 
rozwinięć P i Q, jako funkcye pewnych arytmetycznych wielkości, 
które, jak to później zobaczymy, możemy zawsze bez wielkich 
trudności obliczyć. 


Przejdziemy teraz do właściwego całkowania naszych rów- 


nań różniczkowych zasadniczych. Całkowanie tych równań wy- 
konamy zgodnie z metodą wskazaną przez Brendel'a 1), którą 


') Brendel. Th. d. kl. Pl. I kap. VII; Il kap. Ш. 


ey l 
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z wielkim powodzeniem zastosowali Ludendorff, Buchholz 
i Kramer. W niniejszej pracy ograniczymy się, przy całkowaniu 
naszych równań, do wielkości zerowego stopnia aż do rzędu trze- 
ciego włącznie. 


Zacznijmy od całkowania równania dającego nam 25, gdy się 
ograniczymy do wyrazów rzędu O, wtedy ze wzoru (2) otrzymamy 


dS 
dv ra EM 


ale na zasadzie wzoru (12) jest 


6% = 49, sin 4w, 


a zatem | 

dS і 
(19) en 40, sin 4w 
a ponieważ 


w=(1 — р) — B—wW; 


zawiera w sobie wyrazy okresowe rzędu wyższego, niż zerowy 
(a także i wyraz wiekowy, jak wiemy), to po całkowaniu w S wy- 
stąpią dzięki tej zmienności argumentu w, wyrazy rzędu conaj- 
mniej pierwszego, ponieważ w S, niema wyrazów zerowego rzędu 
kształtu A i C. 

Calkujac (19) otrzymamy: 


(20) So = as + a, cos 4w-|-wyrazy zwykłe. 
Zastanówmy się bliżej nad otrzymanym rezultatem (20). 


Przedewszystkiem wyrazy zwykłe, które oznaczymy Soo otrzy- 
mamy z następującego wzoru: 


= = — ЎА, по SIN nw 
gdzie n= 8, 12, ...; po scalkowaniu wyrazów zwykłych mieć be- 
dziemy: 
(21) se a cos nw; 


n L 
д E 


W wyrażeniu (20) wyraz a, jest to pewna stała"), którą póź: 


l) Brendel. A. N., 3346. 
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niej określimy. Aby wyznaczyć spółczynnik a, zwróćmy się do 
wzoru dającego nam Q w funkcyi spólczynników À oraz argumen- 
tu w (pg. 55) uwzględniając związki (3) pg. 37 (10) pg. 46 oraz pg. 
98 et sqq. przez różniczkowanie otrzymamy: 

4 (1 VL Ui) A dh nir 2 (Lu) 40, ri 16р. (1 Ai ш) Agoo Ti an 


1.0 


dO p) dran + (14200 daer 
- 12 (1 — p) Amon 6, 121—240 — 1) 4, 4,17 62(1 ` IE 


po scałkowaniu otrzymamy: a, = 


A PER A a= 
4.0.0 18.0.0 8.0.0 ez 
MLE AN 5 AGP ESTE SKS 116 (1-1) 
| H | | Dis 3 A RU | 
à m) . KN 4.0.0 12 0.0 d 
v Eh El u th LE EE 
"up ent A 12.00 Aso | x 
2154 Me) Ве I? 
ale 12-6 
| — pa = 2 i 
oraz ROWE dit 
episc | 


Uwzglednijmy te związki i uporządkujmy (a) według potęg 
rosnących 8, wtedy otrzymamy: 


(22) a, =4'14'Bi o" 8... 
gdzie dla krótkości napisaliśmy: 
Aa. 0.0 
VES, 
Ag š 0 8j. A, оо 
aufs) (1-8)? 
(22015) 0 20 | | 
Ai au е T Sr, 


E A IE Te 
Gun A" 16? А 8р? Aan , Di Aan 


12.0,0 4.0.0 


п): ars» 59 QF)" 
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Wyrazenie (20) po uwzglednieniu (22) przejdzie w naste- 
pujace: 


(23) % =, + (9 +q*8+0""B,?) cos 4w Bang, 


W wyrażeniu (23) możemy już wszystkie wyrazy obliczyć, 
z wyjątkiem tylko stałej całkowania ag, którą na końcu wyznaczy- 
my po scałkowaniu równań, odnoszących się do pozostałych spół- 
rzędnych Gyldénowskich. 

Przejdźmy do scałkowania równania różniczkowego, odno- 
szącego się do spółrzędnej R. Ponieważ obecnie chcemy pomi- 
nąć wielkości stopnia pierwszego, więc przedewszystkiem musimy 
zauważyć, że (p), —0, gdyż w (p), jak wiemy”), występują tylko 
"wyrazy krótkookresowe kształtu B, których najniższych stopień jest 
pierwszy. Mieć więc będziemy ze wzoru (1) str. 632): 

а? n dR 
(24) | U R—28, — P,- o do). 
gdyż 
= (p) +R. 

Analogicznie do tego, jak poprzednio pisaliśmy, wielkość (2, 
możemy napisać: 

P,=7,--P, cos 44. 


S, znamy już ze wzoru (23), zatem w (24) pozostaje nam obli- 
czyć tylko wielkość 
dR : 
Lo ER | , ) 


otóż wzór (46) str. 34 daje nam: 


(25) Q= ЎА, aa SIN nu + BR, > Ay o Sin nW—YpKENA, a Cos mu, 


zaś z (12) str. 48 znajdziemy : 


dR dw 
= |= w^ ABY dw 


czyli uwzględniając (10) str. 46 znajdziemy: 


(26) 


»Brendel, T^c ро мо9: 
3) Brendel, |. c. pg. 90. 
3) Brendel, 1. c. pg. 123. 
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(27) do | = — PO +à) sin 4w. 


Pomnózmy;(25) przez (26), KERN 


dR\_ aa 
CINEO, ES | =24,5, sin m | (5, x ВХА"? sin nw E | — 
—u Kin SS A, COS- 210 
0 
ale 
Do 

a wiec 

29) R, Es de) = ê? 0-6) cos 4w „sin 4w=— 484? (1-73) sinu 


KOZA ZY 
mieć więc będziemy: 


123] 


o | do = — ү, On SIN 4w.cos 4w=(cos 810—1) ub (1+ 


(30) K 

W (28) wyrazy krótkookresowe będą te, dla których wskaZ- 
nik » przyjmować będzie wartość 4, 8 lub 12 a więc mieć bę- 
dziemy: 


ак | 


dH 8) 
do do |, 


(31) (А, о T 4800 cos 410) Е; — 


ie Ö 


li Sn Ó HA; ao + Aizos Ii cos 4w + 
+ „(1 i Aso — 34300 (But COS 410. 


Wskutek tego (24) przejdzie w następujące: 


(32) zt A: R=20,+P 410), + (0), 4-0), + (99, 62 + 


+ (W); £,*1 + (Vs 1,2 cos 4w 
gdzie dla krótkości napisaliśmy: 


A ) 1.4 Л 
(38) Do = га В оо 5 7 BE Bi x ° EA Ж а, SZ deg? 


2.0 
m a 2 
2 Бо . 
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75 
Mr 0.0 
(54) (P) =— mak Tr 
AMT TEE 
ç pr! 1.0 8.0.0 au 4 
(3), = m BS А By, "FR d GT 9 Z ы 
Su A 
() = — 4e Boot нт, не 
(Ф), Tal (В — 4 Kiem T Bo? dr 
pe ‚2.0 „a 
ub „OMA fen ON le mu VILI de AM os. 
9(1-1-8,) ' 2(14- A | 4 Lann T 41209 /! 
n 6147 
: Le 10 4. 2.0. 
— (E -- Get 400 Solais 
12? A 36p* À. 
iy — MB à , " | E. | ° 12.0.0 
(Pe H mp 364. Boo Fa | 14-0, 


Ale ze wzoru (12) str. 48 wynika !): 


dR у š 
(35) dv + R=b,—[1—(1+0,)?]f8, cos 4w, 
czyli 
(36) d*R 


dy? + R=b,— (20, + 2,2) B, cos 4w, 

gdzie dla krótkosci piszemy: 
do = 249 90; 
(а) 20, —-0, = — | (4), + (0), 8; + Dan + Q9, 8,? + (0); Br H- 
CLR: 
zamiast y, napiszmy: 
_ 2% 
1-8, 


wtedy analogicznie do równania (23) str. 73 będziemy mogli napi 


UR Dere n del, Le 


‚ pg. 120. Ludendorff, l. c., pg. 17. 
pg. 45. 


Kramer I, 


76 J. KRASBOWSKI. 
sać spólczynnik przy cos 4w, jako rozłożony według potęg py, w po- 
staci następującej : 


pi 2" b Бр"? 


przez porównanie z (a) wynika: 


(37) (20, +) ar Bi --P By 
gdzie dla krótkości napisaliśmy: 
24 
NY WESA Т NE 
(38) p 4.0.0 js 
4° 14-6 Bso 
0 8.0.0 Ly М 
р = Воо 1 Г 3 Boo — ne qe 9 | Aso кыта» ДЕ 
16y. Agoc 
MALTE 


p" = B Hr Boo tt B, Aio lA oo 3 Alano — 
PLAT RA 772 Bio O ei Arne 
+ (ën Ain — 2% А, оо + б 400 — 882 Bos + 
+ 72»? B ARA + 48р? “с 14442 A 0: 


Uwzględniając (36) i (38) otrzymamy: 


, Sg Bt 6? 
(39) R= 20,-- Po — RER cos 4w, 


gdzie znów dla krótkości napisaliśmy: 
/ ) r, |- i ei 
(40) Po =— Bang — Booth + i А оо р Bioota 
di AR : 

Otrzymaliśmy więc wszystkie wyrazy typowe funkcyi #—po- | 
zostaje jeszcze obliczenie tych z wyrazów zwyklych, których | 
w naszym przypadku pominąć nie możemy. 

Możemy bowiem napisać: 


= A 


„2 -+ R=? S, — P,= X (2 Bang cos nw — 13,00 COS NW), 


albo, j oznaczymy 28,00 — Bnoo przez У оо, to mieć będziemy: 
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(41) aus + R = SW „os cos nw. 


Całka równań (41) będzie: 
(42) R, = 9, sin v — g, cos v. H 
gdzie g, i g, sa pewne funkcye, które należy wyznaczyć. 


Z łatwością otrzymamy związki następujące: 


d k | 
do = 5 EV oo cos (nw +) + 1 ХУ, оо cos (nw — v 
Р 

de = ) £V „ое sin (mw + 0) —3 ZWnoosin (nw — v). 


Scałkujmy te równania według metody podanej przez G y l- 
dén'a Brendel'a 2). 
Jest: 


| cos (nw + v)dv =- sin (nw+0) -4- 


l 
n(l— y, )+1 
п р [ dwa 


cos (mw + vidv 
n(1—p)>+>1) do Wii 1 4 


+ 


j sin( n w + v) dv = cos (nw + 0) -- 


all 
n(1—u,) +1 


n AO 20 a oa 
| * sin (nw -3-v)dv 


Tad — u)--1J dv 


a ponieważ W; zawiera tylko wyrazy długookresowe kształtu A 

lub C, a wyrazy kształtu € mogą występować dopiero gdy uwzglę- 

dnimy wielkości pierwszego stopnia?), dla tego też w naszych wy- 

rażeniach w całkach zależnych od W; nie będzie żadnych wyrazów 
| stopnia zerowego. Wstawiając ostatnie dwie równości w dwie po- 
| przednie, znajdziemy: 


1) Kramer, L c. 1, pg. 48. 
? Brendel. Om Användningen $. 11, 
Brendel‘ Th dad. КІРІ. Lupp 120; 
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| 


(43) 


YP o0 
METRES 
P „oo 
n(1—p)= 


b 
pe 


cos (nw + v) + 


rr 


YI 


+32 


q COS (nw— v), 


gdy (43) wstawimy w (42), otrzymamy: 


R, — >O Kon COS nw, 
gdzie | 
V n.0.0 Р | 
Коо = 1— PE | 
ale 


Ш Son — Boo a ln a (I Tal 


wiec ostatecznie otrzymamy: 


2 KIT ca ann 


44 AGO 
(44) 0.0 t +8) y 
| 4 | 


dla wartości n = 4 gdzie S,oo jest określone przez (21). 


Pozostaje nam jeszcze całkowanie równania dającego zwią- 
zek pomiędzy W i v Zajmijmy się obecnie tem równaniem. 
Mieliśmy na str. 20 równ. VI, w którem obecnie uważajmy tylko 
wyrazy stopnia zerowego, mieć tedy będziemy: 


dW 


Go) т 


OST or SLI MS GR ONE EC 


0 > 


Wielkości S, i À, sa nam już znane, jakoż uwzględniając 
związki na str. 48 mieć będziemy: 


So = а, + a, cos 4w; R,=b,+f, cos 4w, 


z których odrazu znajdziemy : 


O RUCHU PLANETOID. 79 


| S Ro =a ba + a, 0, cos 4w-l-a,8, cos 4w + а, cos? 4w; 
R, =D, + 20,0, cos 4w + 8,” cos? 4w; 

(46) ; R3 =b,* + 30,8, cos 4w + 3b,[f3,? cos? Au LBB cos? 4w; 
SoRo = Qo bo + 2a,b,8, cos 4w + a,ßı? cos? 4w + 

+ a, 0,2 cos 4w + 25,0, 8, cos? 4w+ 3,2а, cost 4w. 
Zbadajmy wyrażenia (46). Wielkości: a,, by, a,, b,. sa 
istotnie rzędu masy zaklócajacej m’, a więc Oe, 0,2, a, Da, a b?, 
a,b, sa istotnie rzędu m’? — należy je pominąć!); 3, jest rzędu d 


gdzie à oznacza małą wielkość; zatem mieć będziemy: 


R? = 2b,$, cos 4w + B,2 cos? 4w 
ale ponieważ 


cos? 4w = cos 4w . cos 4w =4 cos 8w LA, 
to 


(47) R = 2b, B, cos 4w + 4 8,? +- B,? cos 8w 

SR, = a B, cos 4w + 1a,8, + $ a,b, cos 8w 

S, R,* == 4 a,B,? + af,” cos 8w +  a,8,?* cos 4w --4 a, B,? cos 12w 
3 — š babi? + š DB,” cos 8w+ 28,? cos 420 + 164? cos 12w. 


Gdy uwzglednimy powyzsze zwiazki, otrzymamy ostateczne 


7 


równ. różniczkowe dające e 


(a 5 = e, tr (a, — 20,8, — 28, 66,8, + 


+ 2 a, 6,2 — 3f, cos 4w — (G, b, — š 8,? -. 
+ $ a Bi? — 6 b, B1?) cos 8w + (3 a, 2,? — 8,3) cos 12w, 


gdzie dla krótkości oznaczyliśmy wielkości następujące: 


1) Brendel w przypadku Hekuby i Hildy pomijał wszystkie wyrazy 
rzędu wyższego, niż drugi, a z wyrazów drugiego rzędu odrzucał te, które są 
rzędu istotnie pierwszego lub wyższego. Ponieważ w moich rachunkach pasu- 
wam dokładność o jeden krok dalej, przeto pomijam dopiero wyrazy rzędu 
wyższego, niż 1m', oraz te, które są rzędu istotnie drugiego lub wyższego. 
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Fir | 
Co = do — a,B, + š Go,” — 60,68, 20. 

W tej wielkości c, występuje jako wspólny czynnik masa m', pod- 

czas gdy w wyrazie $ß,? występuje masa m’ podzielona przez ma- 


łą wielkość ò. Wykonajmy całkowanie równania (47). Mieliśmy 
przedewszystkiem na str. 32 równanie (42): 


W = w 4- WZ + K, 
różniczkujmy ten związek względem v, otrzymamy: 

aw а: , ак 

т te) = i 
NONI dy 


( oznacza część wiekową, a ponieważ dtugookresowe wyrazy 0 
stopnia nie istnieją!) w funkcyi W; ani też tembardziej w K, prze- 
to będziemy mogli napisać: 

dW  —, dK 

(39) do T dv 
Porównajmy (47) i (48). Otóż widzimy, że 


(49) += Co +71» у, 


z drugiej zaś strony, uwzględniając związek (13) str. 49, oraz (10) 
str. 46, przez różniczkowanie znajdziemy: 


dK 


(49) E że aL 


(1 -+ 6,) cos 4t» + 2 (1 + 4,) K^; cos 8w + 
+ 3 (1 + à) К cos 12w, 
gdzie spólczynniki K^, i KE: oznaczają wyrazy zwykłe. 
Podstawmy w (48) związek (49) znajdziemy : 


(50) rt (1 +8,) cos 4w + 2 (1 + 6j) К", cos 8w + 
+ 3 (1 + 8,) K', cos 12w. 


Porównajmy wzór (50) ze wzorem (47) to otrzymamy : 


!) Patrz wzór 14, str. 49, Brendel. Th. d. kl. PL. I, p. 68. 
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| =6+1 
| ES 20,4, — 26, + 65,0; +4 Gd — 96 
| L ] 4-3, 
Шш sze ira à Say + DA 
| : (1 +0,) 
E RE 
| 12 (1 +8) 


Do wyznaczenia psp calkowania jeszcze powrócimy, tym- 
czasem ograniczymy się do uwagi, że ponieważ uważana przez 
nas planetoida Thule należy do krytycznych planetoid pierwszej 
klasy, według podziału Brendel’a!), to dla takiej planetoidy 
mieć będziemy: 

(42m 0, ale 7 =0, 
a więc w naszym przypadku y=y *) i równ. (48) przejdzie po 
scałkowaniu w następujące: 
(92) W = TU + K. 

Przypuśćmy, że stałą a, już wyznaczyliśmy, wtedy uwzgled- 
niając związek b,=2a,--p' na str. 70 i zw. (40) str. 76 będzie- 
my mogli napisać po uporządkowaniu rezultatu według potęg B, : 


b= 2 a4 +9 + 9181 + 030,7, 


gdzie 
Qo ~= — Dooo; 
(53) 1 à 4 
| yl. ++ 0, . H — 
езе Deos 9 Aut Tote 24.00 
Qa = — $ Bio: 


Zatem, gdy znamy stałą a,, to na zasadzie równ. (51) i (53) 
możemy z łatwością obliczyć już poszczególne części równania (52) 
i tym samym otrzymać funkcyę W. ograniczając się ciągle do wiel- 
kości stopnia zerowego. Postępując zupełnie analogicznie do te- 
go, cośmy zrobili gdyśmy szukali wyrazów zwykłych w funkcyi R 
(str. 73), znajdziemy z łatwością podobne wyrazy i dla W, miano- 
wicie będą one: 


) Brendel, I. c., pg. 92, 100. 
2 Brendel, 1. c; Kramer l. $ 1. 
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(54) У „оо = dla 2224. 


Przejdźmy obecnie do wyznaczenia stałych całkowania. Po- 
nieważ obecnie ograniczymy się tylko do badania ruchu planetoi- 
dy w samej płaszczyznie jej drogi, więe mieć będziemy tylko 4 sta- 
łe, które powinny występować w naszych równaniach po ich scał- 
kowaniu. Te stałe całkowania są, jak to wiemy, x moduł mimo- 
środowy, T długość perihelium, a połowa osi drogi planetoidy 
oraz tak zwana stała ruchu n. Jak widzieliśmy po wykonaniu cał- 
kowania wyrażeń różniczkowych spółrzędnych Gyldén'owskich 
występuje jeszcze jedna nowa stała a,. którą, dzięki temu, że jest 
nadliczbową, możemy uważać jako stałą dowolną; tę stałą a, 
określmy za pomocą warunku, że bœ część stała, w wyrażeniu na R 
ma być rzędu istotnie pierwszego— wtedy także i b, będzie takiego 
samego rzędu (str. 70). Ta stała a, będzie tak dobrana, aby część 
stała c, w wyrażeniu (47) przyjmowała wartość zero. Ponieważ 
W =W,-|-K, a K składa się z funkcyi R i 5, to stała część w W 
składać się będzie ze stałych części funkcyi R oraz S; wskutek te- 
go mieć będziemy 
(55) Co 1 — 4,200. 

Ale Brendel!) wykazał, że o ile mamy do czynienia z plane- 
toidą charakterystyczną. to musi zachodzić związek następujący : 

\ 
y 2, = O0= c, + ү; 

dla krytycznych planetoid, a taką jest Thule, y nie może być ze- 
гет ?); a równość 

а — 2 0, — C, zi T 
pozostaje ważną i dla planetoid krytycznych, będziemy więc mieć 
związki następujące : 
(56) 0 =0; 1+0. 

Z pierwszego związku (53) widzimy, że b, da się wyrazić ja- 
ko funkcya liczbowa wielkości a, oraz potęg f,, a więc w przypad- 


ku ruchu planetoidy typu Thule, należącej do klasy pierwszej, 
równania (53) i (56) najzupełniej określą nam wartość stałej a,. 


i) Brendel, |. c., rozdz. V sparsim. 
2) Brendel, Astr. Nachr. 3346; id. Th. d. kl. Plan. I rozdz. VII spars. 
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Wielkość W (str. 37) możemy napisać w postaci następu- 
jacej: 
(57) W —(C,--4)v + wyrazy okresowe, 


czyli 
W —=(C, 1-1 +- fo) + wyrazy okresowe, 


gdzie jak zwykle 
= + Y 
ale y, jest przynajmniej stopnia drugiego, więc obecnie możemy 
je opuścić, zatem 
(58) 1=1- 
Mieliśmy związek następujący: 


фу MU efe 
i 
w,=(1—)0— B — x. W, 
a wiec 
v=w -LK=(1—u)v— B=pW+yk, 
ale 
W —W, +T v 
w= (1— р) 0 туо + wyrazy okresowe = 
= (1 — к н) + wyrazy okresowe, 
gdzie 


po =p (l +T); 
zatem część wiekowa W jest następująca: 
(59) pars sec w= (1 + ko) v. 


Przypomnijmy sobie określenie wielkości y i p.;. 
Mieliśmy związek następujący: 
s т’ п 
60) RTR "WELT 
(60) DO Labo E 
We wzorach (60) jak wiemy z poprzedniego, przez n ozna- 
czamy stałą całkowania; wielkość, którą Brendel nazwał „stałą 
ruchu“; m, zaś oznacza rzeczywisty ruch średni dzienny. 
Wprowadźmy dalej wielkości à і à zdefiniowane w sposób 


następujący: 


aT ы RA 


pe: 
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Z wymienionych przez nas małych wielkości à i ô, występu - 
је w naszych wyrażeniach, jako dzielnik tylko ö,. 

Brendel wykazał, że 0, nie może być zerem w przypadku 
Ścisłej wymierności, jest ono rzędu m o wykładniku < š przeto za- 
wsze ono jest znacznie większe od masy zakłócającej. Co zaś do 
wielkości à, to ona jest rzędu m’ z wykładnikiem x 1 a w przy- 
padku wymierności staje się nawet równą zeru. Otóż zasadnicza 
myśl metody podanej przez Brendel'a polega właśnie na tem, 
że w wyrażeniach perturbacyi występują, jako małe dzielniki nie 
à lecz ô. Brendel wykazał, że stosunek ścisłej wymierności mo- 
że zachodzić tylko pomiędzy wielkościami m i n' (a nie pomiędzy 
n, i n), wtedy 0=0. 

Na zasadzie wzorów poprzedzających mieć będziemy *): 


i 10 =p dsl Pep 
3—8,—(8 — 8) (1 Lake 
A=". Ж 46746 — 3 Co — 8%; 


ale 


ponieważ mamy do czynienia z planetoidą krytyczną, dla której 
delt ban TIE" 
à, —:01-2-6 — Зу, 


4y — 3 — 6, | 


to mieć będziemy: 


ale ponieważ 


to ostatecznie otrzymamy : 
(61) 0, =0—47u. 
Zupełnie analogicznie możemy napisać, że 
оз e ОА 
"Wr deser ТЗ, 
gdy uwzględnimy (59), to n, możemy uważać jak ruch średni pla- 
netoidy w długości. Postępując zupełnie analogicznie jak poprze- 
dnio robiliśmy, otrzymamy: 
(62) 0,—0 — 4 u y == —4 Ly, 


Gdy 8, staje się bardzo małem lub nawet zerem, wtedy za- 


Y Brendel, A. N. 3346, 
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chodzi przypadek libracyi, który w teoryi Gyldén'a jest łatwy do 
uwzględnienia. m, właśnie określa nam ten ruch średni dzienny, 
który powoduje libracyę; musimy też znać dokładnie tę wielkość т, 
aby módz wyznaczyć czas obiegu planetoidy. Należy zauwa- 
żyć, Zen, im, różnią się tylko o stałą yy. Wielkość y, jest dru- 
giego stopnia, a więc obecnie ją opuszczamy, podczas gdy y = 1 
jest stale dodatnią i w przypadku planet charakterystycznych 
pierwszej klasy, do jakich zresztą należy Tliule, jest już wielko- 
Ścią stopnia zerowego, a więc w obecnych warunkach musi być 
ona uwzględnioną. 


V. Zastosowanie liczbowe otrzymanych wzorów 
do ruchu planetoidy (279) Thule. 


Planetoidę Thule (279) jako gwiazdę 13.5 odkrył w Wiedniu 
J. Palisa dnia 25 października 1888 г. Po kilkunastu dniach 
obserwacyi, ten astronom odrazu się zoryentował po wstępnym ra- 
chunku, że droga tej planetoidy leży w bardzo znacznej odległo- 
ści od słońca, mianowicie prawie 4 razy dalej niż średnia odległość 
ziemi od słońca. Przypuszczenia Palisy były potwierdzone 
przez rachunki Langego z Berlina i Halma w Kiel. Nieby- 
wałe dotąd wśród małych płanet rozmiary drogi nowo odkrytej 
planetoidy wskazywały, że ta planetoida, znajdująca się na ze- 
wnętrznej części pierścienia planetoidalnego będzie bardzo cieka- 
wym objektem dla badań teoretycznych. Można było bowiem 
przewidywać, że planetoida Thule z powodu swego znacznego od- 
dalenia od Słońca oraz z powodu znacznego zbliżenia się do Jowi- 
sza — może podlegać okresowo nadzwyczaj silnym zakłóceniom 
spowodowanym przez przyciąganie Jowisza. Z tych wielkich per- 
turbacyi ruchu Thule, znając dokładnie jej pierwiastki drogi, 
z czasem można będzie obliczyć masę układu Jowisza, analogicznie 
jak to zrobił Newcomb, badając ruchy (33) Polyhymnii ?), ale 
uczynić to będzie można, jak to słusznie zauważył Bidschof?), 


у S. Newcomb. Astronomical Papers prepared for the use of the 
American Ephemeris vol. V, part V. 

2) Sitzungsberichte der Ak. der Wissenschaften, Wien Bd. 100: 
pg. 937. 
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dopiero wtedy, gdy mieć będziemy przynajmniej 25 lat obserwa- 
cyi tej planety. Pomimo słabego blasku (1.35—15.5) jakim świe- 
ci Thule, wskutek czego może być ona obserwowaną jedynie przy 
pomocy wielkich narzędzi, obserwowano tę planetoidę podczas 
wszystkich opozycyi, które miały miejsce od chwili odkrycia, t. j. 
25/X 1888 r. aż do 10/11 1891 r. Te wszystkie obserwacye zostały 
bardzo szczegółowo opracowane przez Bidschofa!), który na 
ich zasadzie podał wartości pierwiastków drogi tej planetoidy, 
oraz efemerydy dla czasów bliskich opozycyi od czasu odkrycia do 
r. 1894. W pracy swej wspomina Bidschof, że 3 razy w ciągu 
stulecia następuje znaczne zbliżenie pomiędzy Thule i Jowiszem, 
w ruchu planetoidy występować muszą bardzo znaczne pertur- 
bacye. 

Na zasadzie wszystkich obserwacyi, które mu były dostępne, 
wyprowadził Bidschof następujące ostateczne elementy: 


Epoka: 1881 Luty 20.0 cz. Sr. Berl. 
ois 1049-2130 +8 
M = 155 36 48.8 
308 44 42.0 


45 => 

Q-— 75 26 12.1 

w — 9233 18 299 | równ. 1890.0 

i= 2 2 342 \ 


p = 4 43 14.2 
log.a — 0.629 6674 
u = 403”. 1860, 
oraz wyznaczył wielkość Thule dla średniej opozycyi, 13"76. We. 


dług obliczeń Bidschof'a, najkorzystniejsze opozycye bywają 
w lipcu i sierpniu; wtedy planetoida osięga 13" 3, ale zato, dla 


obserwatoryów leżących w większych szerokościach północnych, 


jest ona wtedy dość niekorzystnie położona na niebie, gdyż jej 
deklinacya wówczas waha się pomiędzy — 24.5 a — 15°. 5. Ele- 
menty, które wyżej przytoczyłem zostały wyprowadzone na pod: 
stawie 8 miejsc normalnych. Te miejsca normalne są dość do- 


1) Sitzungsberichte d. Ak. d. Wiss. Wien Bd. 100, pg. 937. 
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brze przedstawione przez ostateczne elementy; największe odchy 
lenie wynosi dla 

cos 6d a. — +3". 1; 
za$ dla 

dà = — 3". 0. 


W r. 1908 A. Wedemeyer!) ogłosił nowe rachunki, odno- 
szące się do planetoidy Thule. Powodem podjęcia na nowo tych 
rachunków było z jednej strony to, iż od chwili odkrycia do roku 
1906 Thule dwa razy obiegła swoją drogę, z drugiej zaś strony, 
dzięki pracy Bidschofa, który wprawdzie uwzględnił w swych 
rachunkach obserwacye do 1891 r. lecz podał efemerydy aż do r. 
1894 i dalej, zdołano podczas wielu opozycyi zaobserwować (279). 
Ze względu na ten nagromadzony materyał, A. Wedemeyer 
zdecydował się, wychodząc z elementów Bidschofa, opracować 
wszystkie obserwacye, jakie tylko były zrobione od r. 1888 dor. 
1906 (ogółem 53), uwzględnić dokładne pozycye gwiazd porów- 
nawczych i na podstawie tego całego materyału ugrupowanego 
w 14 miejsc normalnych wyprowadzić nowe elementy. 

W swoich rachunkach Wedemeyer posługiwał się bardzo 
dogodną metodą Oppolzer'a, obliczania perturbacyi pierwiast- 
ków drogi. Ostateczne elementy, otrzymane po wyrównaniu 
miejsc normalnych, były następujące 2): 


Epoka i oskulacya 1907 grudzień 6.5 cz. $. B 
M = 210° 15' 47". 26 


WII 999v E28: 520159 | 
1910.0 


Q — 75 35 57.14 
¿= 2 22 29.93 | 
p = 4 37 3594 


u — 404" 29204. 


Obserwacyi późniejszych nie uwzględnionych przez Wede- 
meyera wykonano niewiele, mianowicie w 1906 r. dn. 17/X 
w Wiedniu i w 1911 r. w dniach 20i 30 kwietnia fotograficzna 
w Heidelbergu š). Planetoida coraz bardziej przybliżała się do Jo- 


1) Archiv. der deutschen Seewarte XXXI Jahrgang. 
2) Elementy Wedemeyer'a służą za podstawę niniejszych rachunków. 
3) A. N. 4591, 4563. 
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wisza przez co perturbacye jej biegu stawały się coraz znaczniej- 
sze i zgodność pomiędzy elementami Wedemeyera a rzeczy- 
wistością stawała się gorszą. Krytyczne zbliżenie Thule i Jowisza 
następowało w r. 1912. W lipcu 1913 r. Viljev!) ogłosił nowe 
oskulacyjne elementy Thule. Te jednak elementy nie mogą być 
uważane za bardzo dokładne, gdyż obliczone zostały wprost z ele- 
mentów Wedemeyer'a przez dodanie tylko perturbacyi od osku- 
lacyi 1907 r., które obliczył przy pomocy jednej z metod perturba- 
cyi szczególnych bez uwzględnienia obserwacyi z lat 1906 i 1911. 
Viljev podał eiemerydy na podstawie swoich elementów. Dla 
chwili opozycyi miejsce wskazane przez Viljev a znacznie się 
różniło od tego, jakie było podane przez B. J. na podstawie ele- 
mentów Wedemeyera, planetoidy jednak nie zaobserwowano. 
Elementy podane przez Viljev'a są następujące: 


Epoka i oskulacya 1913 czerwiec 17.5 cz. Sr. Berl. 
M = 358" 35' 20". 7 


w — 220 43 38.9 
1 065200] „B 2:6 

p = 3 39 490 

„ == 397". 60002 


log a = 0.6337068. ' 


Aby więc zbadać ruch tej ciekawej planetoidy, podlegającej 
tak znacznem zaburzeniom, wydaje mi się, że jedyna droga jaka 
nam pozostaje, jest podać, przynajmniej w przybliżeniu, wyrażenia 
analityczne perturbacyi ogólnych aby umożliwić na długi czas 
obliczanie efemeryd, pozwalających obserwatorom odnalezienie 
tego ciała niebieskiego. Metoda Gyldena-Brendel'a wydaje 
mi się tutaj najodpowiedniejszą ze względu na to, że odrazu, przy 
pierwszem obliczeniu perturbacyi, pozwala uwzględnić wpływ 
wyrazów zależnych od wyższych potęg masy Jowisza. Gdy nowy 
materyal obserwacyjny zostanie zdobyty, to wtedy okaże się ko- 
nieczność obliczenia nowych pierwiastków drogi Thule, a rozpo- 
rządzając już całym materyałem obserwacyjnym, obejmującym 


у A.N. 4661. 
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znaczny przeciąg czasu, będzie można dopiero wówczas podjąć 
zagadnienie wyznaczenia; na zasadzie perturbacyi Thule, masy 
układu Jowisza. 

Do naszych rachunków użyliśmy elementów Jowisza według 
„Connaissance des Temps pour 1914“. Masę Jowisza m' wzięłem 
według Newcomb'a, wynosi ona: 


lg m’ = 6.9799061). 


yrefta Neg | 
"n PRIUS 


Ponieważ w niniejszej pracy ograniczymy się do wyszukania 
przerwy (luki) dla planetoid typu š , pozostawiając do następnej 


części wyrażenia liczbowe spółczynników dających rozwinięcie 
perturbacyi Thule, dla tego też przytoczymy tutaj tylko rachun- 
ki, które odnoszą się bezpośrednio do tego zagadnienia. W roz- 
dziale drugim podaliśmy, zgodnie z teoryą, którą do naszego 
przypadku stosujemy, wyrażenia analityczne rozwinięcia funkcyi 
perturbacyjnej, w których występowały pewne arytmetyczne wiel- 
kości, dające się wyrazić przez stosunek a połowy osi drogi Thule 
do połowy osi drogi Jowisza. Te spółczynniki liczbowe były to 
wielkości f i y zaopatrzone odpowiedniemi wskaźnikami. 

Przystąpmy do ich obliczenia. Z elementów Jowisza i Thu- 
le znajdziemy wielkość a: 

lg a= 9.912 6572. 

Ta wielkość mało się różni od jedności i widzimy, że spół- 
czynniki 8 i y, zależne od tej wielkości, nie będą się znajdować 
ani w tablicach Masala"), w których mamy tabułowane wartości 
B i+ dla argumentu lg a od 9.30 aż do 9.60, ani w tablicach Gyl- 
dén'a*), w których jest zawarty lg a od 9.620 do 9.800. 

Musimy więc bezpośrednio z elementów obliczyć wszystkie f. 
Na str. 24 mieliśmy wzór (8): 

1) Masal H.: Tables de l'intégrale | a Lo fe 

о fl—a? sin? en 
ska iakttagelser och undersókningar Anstalda pá Stockholms observatorium, Fje- 
rde Bandet, Stockholm 1891, pg. 51. 


2) Gyldén H.: Astr. Gesell. Publ. 21. Hülfstafeln zur Berechnung der 
Hauptungleichheiten der Kl. Planeten. 


Astronomi- 
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I de 2 jp sin?" pde 
m, 1 
geret. 
Znając a możemy przystąpić do obliczenia wszystkich B° . 
Gyld&n!), Harzer?) obliczali te wielkości w ten sposób, 
że przedewszystkiem znajdowali wszystkie EH dla n=l. Przy po- 
mocy pewnych wielkości pomocniczych f,, À, i 34, określonych 


w sposób następujący: 


di (2n +3)? a? DR PP m) | 
"  Qn-F2)2n-F4) (1a * C 2n+2 a pr? 
] But: 


Po, fn+m—1 
1 — nm Un + mech) | 
Użycie tego wzoru jest następujące: najpierw oblicza sie y, 
z największem n zachodzącem w naszych wzorach, wprost ze 


WZOTU: 
^ ao 
Metis. 
Ie pv fusta 
ARRET m-f 1 


a potem przez koleine stosowanie wzoru redukcyjnego 
| 


1 — fn Ўл 


otrzymuje się wszystkie poprzedzające y. 
Wprowadźmy jeszcze wielkość pomocniczą Pn = An Bn, wte- 


dy mieć będziemy wzór następujący: 
Pn =PB'oPoD, Pa +++ Prai, 


Ma ZZ 


1) Gyldćn H. Orbites absolues I, livre Ш. 
7) Harzer P. Untersuchungen etc. p. 26. 
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który może posłużyć do tego, aby po znalezieniu 8! obliczyć 
wszystkie ß!„. Potem wszystkie inne В, otrzymamy na zasadzie 
następującego wzoru rekurencyjnego: 
D B FR 9 rem 

Trudność w obliczeniu wielkości В polega na tem, że uła- 
mek ciągły, wyrażający %„, jest w niektórych przypadkach mało 
zbieżny. Harzer zauważył, że dla planetoid zwykłych wystar- 
czy przyjąć ostatnie n = 12; ale w naszym przypadku stosunek a 
jest mało różny od 1, wielkości f„ maleją bardzo powoli, nawet 
dla n — 30 jeszcze nie możemy osiągnąć tego, aby dwa nastepuja- 
ce po sobie przybliżenia w rozwinięciu ulamka?), różniły się od 
siebie o wielkości zaniedbywalne, musimy więc, aby módz obli- 
czyć Bi, postąpić inaczej, bez obliczenia wielkości pomocniczych 
An i fa. W tym celu przedewszystkiem obliczymy spółczynnik Bi, 
możemy to zrobić najlepiej przy pomocy metody średniej arytme- 
tyczno - geometrycznej. Na tę metodę poraz pierwszy wskazał 
C. F. Gauss?), jako na praktyczny sposób obliczania całki elip- 
tycznej zupełnej, mamy bowiem: 


T 


dy 


9 2 
jako pierwsze przybliżenie kładziemy 
c=]; y=Vl—a? 
obliczmy z tych dwu wartości średnią arytmetyczną i geome- 


tryczną: 


mises irs y=Vxy, 


następnie z x’ i y’ znowu średnią arytmetyczną i geometryczną: 
/ , Lee 
TV y" AS Voy Ne d 
aż przy pewnej wartości m znajdziemy: 


ym — qc 


1) Gyldén H. Uudersókningar af theorien for himlakropparnas rórel- 
ser. ll, str. 23, (Bihang till sv. Vet.-Ak. Handl. Bd. 6). 

2) Hansen zauważył, że dla n dużego de = | + œ. 

3) Gauss C. F. Werke III, pg. 352. 
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wtenczas: 
== zm) => "C ` 
Ta metoda obliczenia ß', jest bardzo praktyczną i nadzwy- 


czaj prędko prowadzi do celu. 
W naszym przypadku mamy: 


lg x=0; 16 И1—о°:== 9.7600268 
т!Ү —9.8882366 y! = 9.8882666 
lg Bi — 0.1117634. 


Najwyższa wartość wskaźnika n jaka nam jest potrzebna jest 
n=12, obliczyć więc musimy В,! dlan=1..., 12. 

Te spölczynniki otrzymamy najpredzej przy pomocy metody 
podanej przez Ludendorfffa!). Tablice pomocniczą podaną 
przez Ludendorffa rozszerzyłem tak co do rozciągłości, jak 
i stopnia dokładności i przy pomocy jego wzorów znalazłem sto- 
sunkowo łatwo pozostałe 81, (= 1, 2... 12). Dla kontroli drugi 
raz te same spółczynniki [!„ (n=0, 1, 2... 12) zostały obliczone 
przy pomocy kwadratur mechanicznych. Tę metodę do obliczenia 
spółczynników B, poraz pierwszy zastosował Harzer, później zaś 
O. Callandreau. Dokładność w rachunkach jest znaczna i kon- 
trola łatwa, gdyż każdy spółczynnik jest obliczony jako średnia 
arytmetyczna z dwóch sum mało różniących się od siebie; 2e zgo- 
dności tych dwóch rezultatów widzimy, czy nie został popełniony 
jakiś gruby błąd rachunkowy. W naszym przypadku kąt 90° dzie- 
litem na przedziały ро 5° і dopiero wykonałem całkowanie. Po- 
nieważ cały rachunek wykonałem według wzorów podanych przez 
Harzer’a i O. Callandreau, dlatego też w bliższe szczegóły, 
odnoszące się do tej metody, nie wdaję się, odsyłam do rozpraw 
skąd one zostały wzięte 2). 

Wartości 81, obliczone za pomocą kwadratur mechanicznych 
wypadły zupełnie zgodne z wartościami obliczonemi przy pomocy 
powyższej tablicy. Aby otrzymać pozostałe spółczynniki f*,,, 6%, 
B^»; Pro, Bis, które mają posłużyć jako punkt wyjścia do użycia 
wzoru Íl, znów zastosowałem metodę kwadratur mechanicznych, 


a więc 


') Ludendorff. A. N. Bd. 160. 
2) Callandreau. Journal de l'Ecole Polytéchnique, li série, t. VII; 
Harzer, 1. c.; pg. 26. 
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oraz kontrolę przy pomocy tablicy Ludendorff'a; zacząłem od 
obliczenia powyższych spółczynników dlatego, ponieważ wedle 
uwagi Harzer'a!) najkorzystniej obliczyć spółczynniki 8, naj- 
pierw dla najwyższej wartości n, gdyż w ten sposób najlepiej 
zmniejszymy wpływ błędów, pochodzących od zaokrągleń. Gdy 
już spółczynniki 8; zostały znalezione, zastosowałem wzór II 
i w ten sposób otrzymałem wszystkie pozostałe spółczynniki. Aby 
skontrolować cały rachunek obliczyłem kilka dowolnie wziętych 
spółczynników metodą kwadratur. Rezultaty otrzymane różniły się 
tylko o pare jednostek w 7-ej cyfrze logarytmów, co dla naszego 
przypadku, gdzie się ograniczamy do rachunku 6-cio cyfrowego, 
jest najzupełniej obojętne. 

W ten sposób znalezione logarytmy spólczynników f. są па: 
stępujące: 

Logarytmy spółczynników £, 


n/s 1 3 9 7 9 
0 0.111763 0.388660 0.719219 1.092421 1.495190 
] 9887778 0.236762 0.620523 1.027947 1.451216 
2 9.761532 0.153599 0.563644 0.986993 1.420154 
З 9.693387 0.102476 0.524368 0.955978 1.394980 
4 9.642434 0.062590 0.492468 0.929894 1.373181 
9 9.601597 0.029631 0.465363 0.907191 1.353821 
6 9.567459 0.001428 0.441683 0.887002 1.336327 
7 9.538099 9.976724 0.420599 0.868721 1.320326 
8 9.519324 9.954718 0.401561 0.852056 1.305577 
9 9.489344 9.934848 0.384186 0.836699 1.291872 
10 9.468604 9.916727 0.368193 0.822446 1.279062 
11 9.449703 9.900062 0.353368 0.809138 1.267029 
12 9.439337 9.884629 0.339546 0.796650 1.255680 


Ш 


{n.a 


otrzymamy spółczynniki Yu... 


) Harzer Р, |. c. pg. 26. 


225 


Gdy już znaleźliśmy wszystkie wielkości I, to na 
wzoru następującego 
du Y are 


| +20+1 (204-1) 
A ch Bato 


zasadzie 
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Te spólczynniki sa nastepujace: 


Logarytmy spólczynników ze. 


n/s 0 1 2 3 4 
O 0.024421 9.673703 9.700962 9.839429 0.093954 
1 9713092 9583198 9.574342 9.726009 9917951 
2 9.499504 9.364732 9.455100 9.615956 9.819414 
3 9.344016 9.237503 9.340652 9.508494 9,709071 
4 9.205720 9.117201 9.299699 9.402800 9.606979 
5 9.077540 9.001658 9.121202 9.298793 9.505931 
6 8.956059 8.889613 9.014822 9.196093 9.405778 
7 8.839356 8.780260 8.910104 9.094497 9.306403 
8 8.726239 8.673047 8.806768 8.993846 9.207710 
9 8.615916 8.567583 8.704600 8.894015 

10 8.507834 8.463575 8.603436 8.794901 

11 8.401590 8.360800 8.503145 8.696416 

19 8.296881 8.259083 8.403622 8.598585 


lg тоо = 9.859395 . 


Znaleźliśmy w ten sposób spółczynniki f, i qwa, możemy 
obecnie przystąpić do obliczenia wielkości A i B opatrzonych od- 
powiedniemi wskaźnikami u góry i u dołu. Te wielkości podsta- 
wowe 4 i B mogą być obliczone dwoma sposobami. Jeden stano- 
wi zastosowanie wzorów (20) i (24) (str. 26, 27) przy pomocy 
wielkości pomocniczych 9, m, p które sa funkcyami yn, „, oraz wiel- 
kości Pn m рі Qu, m p, Znów będących funkcyami Q, m >. Podo- 
bne postępowanie Brendel zaleca; jest ono istotnie bardzo prak- 
tyczne, gdy chcemy odrazu obliczyć wszystkie potrzebne spół- 
czynniki w rozwinięciu funkcyi perturbacyjnej dla danćj planetoi- 
dy. W naszym jednak przypadku podobny sposób postępowania 
nie zalecałby się, gdyż ograniczając się do wyrażeń zerowego 
stopnia, daleko prościej odrazu obliczyć spółczynniki A i B wprost 
we funkcyi Т. Wzory stosowne dla takiego rachunku zostały po- 
dane przez Masal'a!). Są one wygodne w użyciu o ile nie posu- 


1) H. Masal: Formeln und Tafeln zur Berechnung der Absoluten Stör- 


gen der Planeten, Kongl. Svenska Vetensk. Ak. Handlingar. Bd. 23. 
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wamy sie w rozwinięciach do wyższych stopni, gdyż wtedy stają 
się one skomplikowane. W granicach dokładności przyjętych przez 
nas, wzory Masala są następujące: 


Ann = — 22 Yno 
Bnoo = -|- 2N Yno +4 Tn. 
Än = Yno [002 4-20) — 2n? q] -+ 2n Yn. 
Bho = mol (m+n) 2n? e] 4-14 [— (4m + 6) 4- 4n q]. 
— 8 Yne 
A = R (02 2n) + 2n? g] A qm. 
LL = ho [— (12410) 229] 4-1 [ — (4n 4-6) — 4n q] — 
— 811.2 
Ac. = Yno (2n?--4n)4-4n ts 
Dh, = — Tao 2n (n 4- 1 — 1,4 (8n + 12) — 164,5 
Apo = - Yno(n + 5n? 4-68) — 144 (4n? + 14n) — 8n y,» 
RI, = + quon (n-4- 1) (n-1-2)4-141 6 (n -4- 2)? + 
41.2. 24(n--3)--484,.3. 
gdzie p = n Na zasadzie tych wzorów otrzymali$my nastepu- 


jące wartości logarytmów spółczynników 4 i B. 


n 0 4 8 12 
Алоо 0.108810, 9.930359, 9.677099, 
В, оо 0.075763 0.257352 0.017139 9.738 837 
de 0.991403 1.001115 0.876663 
D. 1.136647, 1.173232. 1.088907, 0.931 817, 

Bäi 1.682384, 1.834337, 1.818637, 
IE 1.906612 1.946349 1.945695 1.972880. 
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Znaleźliśmy więc ostatecznie potrzebne spółczynniki rozwi- 
nięcia funkcyi perturbacyjnej; możemy obecnie przystąpić do wy- 
rachowania przerwy dla planetoidy typu Thule w pierścieniu pla- 
netoidalnym. 

Mieliśmy równanie (37): 


(1) (20, +24) 8, =p'>p"B,+p""B,?, 


które zawiera tylko jedną niewiadomą f, (à, wyraża się przez f) 
przy pomocy związku (61) str. 84 i jest względem tej niewiadomej 
stopnia trzeciego; (jak zobaczymy spółczynnik przy 8, w tym rów- 
naniu zawiara 0). 

Brendel okazał, że w celu znalezienia obu granic przerwy 
należy wyznaczyć ten spółczynnik, przez nadanie odpowiedniej 
wartości na 0, w taki sposób, ażeby to równanie dawało na f, tylko 
dwie, od siebie różne wartości, wówczas dwie wartości 6,, odpo- 
wiadające tym dwum różnym wartościom f,, określają właśnie gra- 
nice przerwy. 

Brendel!) wykazał, że B, w przypadku planetoidy charak- 
terystycznej, nigdy nie może się stać większem co do rzędu wiel- 
kości od pierwiastka trzeciego stopnia z masy zakłócającej; gdy 
mamy do czynienia z planetoidą tego rodzaju, że stosunek ruchu 
średniego dziennego planety zakłócającej do takiegoż ruchu dzien- 
nego planetoidy wyraża się ściśle przez prosty ułamek, wtedy 


6—0, zaś 0, = — 4h; \ 
ale jak to ju2 poprzednio widzielismy 
3 
= B 


Ponieważ f, w let planetoidy charakterystycznej co 
najwyżej inoże być rzędu ra , przeto najmniejsza wartość (à,) 


jest rzędu mniejszego (co do rzędu wielkości) niż ef Jeżeli mia- 
nowicie ta wartość graniczna byłaby przekroczoną, to wtedy nasze 
rozwinięcia stałyby się rozbieżne”), ruch planetoidy, któraby mia- 


ła takie ó, byłby niestały. Dla tych krytycznych wartości 2, : 


! Brendel. A. N. 3346. 
2 Brendel, l.c; Brendel. Theorie d. kl. PI. I, Kap. VII; 
Gyldćn Acta Math.: t. IX. 
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a więc if, występuje (w znaczeniu Brendlowskim) przerwa (la- 
kuna) w pierścieniu planetoidalnym; przerwa ta przychodzi do 
skutku w ten sposób, że 2, zrywa ciągłość; wartości jego po obu 
stronach przerwy są przeciwnego znaku. Ponieważ 0,, jak to wy- 
i 
nika z poprzedniego, nigdy nie może być zerem, przeto H. = а 
1 

nigdy nie może się równać jakiemuś prostemu ułamkowi. Nato- 


miast ó może się stać równem zeru: ma to miejsce wtedy, gdy u 
jest równe prostemu ułamkowi. dla. tego też we wszystkich na- 
szych dotychczasowych wzorach, gdzie zachodzi mały mianownik 
?,, nie wolno zastąpić go przez ô, jakkolwiek 8 może być bliz- 
kie 9, . 

Stosownie do tego, cośmy wyżej powiedzieli. zajmijmy się 
przedewszystkiem wyznaczeniem wartości D. W tym celu weźmy 
na uwagę przytoczone przed chwilą równanie (1). Ponieważ Š, 
nie jest rzędu zero, to 2,78, jest rzędu wyższego, niż pozostałe wy- 
razy wtym wzorze, możemy je więc opuścić. Zamiast ó, pod- 
stawmy 

0 — 4p y 0 — ôu p, >, 


wtedy otrzymamy: 
268, — 12 p =р' 4-р” 8, + p"8,?, 
co możemy napisać: 


REN Aa p pM 


| f | 2 RE 
czyli 
(3) Pr FB JB PELO, 
gdzie dla krótkości położyliśmy: 
"iu ei "Y 2 28. ¿À p' 
0) SNIEG is Gino ln czaso mou 


Z (4) widzimy, że tylko spółczynniki ż oraz k są niezalażne 
od 6. Równanie (2) musimy tak przekształcić, aby sie pozbyć wy- 
razu B,?^; w tym celu połóżmy: 

і 
(а) B, = b, — 3 
wtedy znajdziemy: 
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(5) р tab, te=0, 
gdzie oznaczylismy: 
NEE, 
(6) d=J — 3 
oraz © 
бега S 1% - Z + К. 


Równanie (4) będzie posiadało pierwiastek podwójny wtedy, 
gdy wyróżnik równania (5) jest równy 0. Po obu stronach przer- 
wy, b, ma wartości przeciwnego znaku; ponieważ B, —b, jest małe 
w porównaniu 2 /,, przeto także i wartości ß, będą przeciwnego 
znaku. 

Wyróżnik ten jest 


4 6 
pr | 


z teoryi równań trzeciego stopnia wiemy, że gdy D>0, wtedy 
mieć będziemy tylko jeden rzeczywisty pierwiastek rówńania (5), 
a więc dla ß, otrzymamy tylko jedną wartość; gdy D=0, wtedy 
otrzymujemy trzy rozwiązania, z których dwa są równe i wartość f, 
znaleziona z tego warunku daje nam właśnie krytyczną wartość à;; 
wreśzcie gdy D < 0 mieć będziemy wszystkie trzy pierwiastki rze- 
czywiste, ale tylko jeden, jak to możemy się przekonać przez pod- 
stawienie, spełniać będzie nasze równania. 

Aby znaleźć ten krytyczny warunek, z otrzymanych poprze- 
dnio rozwinięć spółczynników spółrzędnych G yldénowskich 
obliczmy stałe wielkości 2 oraz k; w naszym przypadku znale- 
zlismy: 

lg i =8.442234_,9; 1g —6.672816. 4, 


te wielkości są stałe; wielkość d zależy od j, a więc od 6, po sze- 


regu prób dokonanych co do wartości j znajdujemy, że warunko- 


wi D —0 czyni zadość ta wartość j, która daje ze wzoru (5) 
lg d=8.12508,_10 . 1 lge=7.773160_vo, 


skad znajdujemy: 
lg фу == lg bL,” = 8.824044 
za$ 
lg 5," — 9.125074, 
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na zasadzie wzoru (a) z łatwością otrzymamy: 
Ig(81)i — 8.759360; lg(b) = 9.154128, 
te wartości ß, są krytyczne. Znając wartość krytyczną j z latwo- 
ścią znajdziemy tej wartości odpowiadającą wartość 6. W naszym 
przypadku znaleźliśmy: 
lg à = 8.559884 ; 


(B): 1 (Br)u 
otrzymamy dwie wartości dla y na zasadzie następującego wzoru, 
dokładniejszego niż ten, który poprzednio stosowaliśmy !): 


l 
ee ur Pira. 
© "= 1—60) 


znając 


Wartości te w przypadku naszym są następujące: 
lg yı = 7.696539; lg yu = 8.495577. 


Teraz ze wzoru (2), rozdziału poprzedzającego, znajdziemy 
dwie wartości 0,, odpowiadające tem wartościom y, mianowicie: 


Ig(ö,)ı =8,334772; 1g(8, )ı = 8.750883, 
dalej ze wzoru 
en 
93-8, 
znajdziemy obie wartości prawdziwego ruchu średniego, które 
ograniczają szukaną przerwę dla typu */,. Te wartości są: 


(n) = 401". 732; (n )u=391". 485 


nı 


wreszcie stałą n, dla wartości krytycznych ö,, znajdziemy ze wzoru 
następującego: 


skąd: 
n= 403". 723, 
Rezultaty otrzymane możemy więc wysłowić w sposób na- 
stępujący: 
Dla planetoid typu ?/, wartość B, musi spełniać 
warunek następujący: 


' Brendel A. N. 3346 wzór № 38. 
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lg B, <8.759360 lub lg 9,2 9.154128, 
zaś 
lg 0,<8.334772 lub lg 2,28 750883, 


zatem prawdziwy ruch dzienny n, dla typu Thule 
musi być 
n, > 401”. 732 lub n, < 391". 485. 


Pomiędzy temi wartościami m, znajduje się przerwa w zna- 
czeniu Brendelowskim w pierścieniu planetoidalnym; planetoida, 
która posiadałaby ruch średni dzienny prawdziwy (w definicyi 
Brendelowskiej), zawarty pomiędzy powyższemi granicznemi war- 
toSciami n, poruszałaby się po drodze niestałej. 

[Musimy tu zwrócić uwagę, że ze względu na definicyę n, ta 
wielkość nie jest bezpośrednio porównywalną z wartością ruchu 
średniego dziennego oskulacyjnego, jaki znamy z elementów pla- 
netoidy]. 

W ten sposób otrzymaliśmy wartości krytyczne prawdziwe- 
go ruchu średniego dziennego, warunkującego przerwę dla typu */,; 
taka przerwa dotąd była wyznaczoną dla Hekuby i Hildy (zrobił 
to Brendel?). Prace Schwarzschild'a dla Hekuby‘) i Dzie- 
wulskiego?°) dla Hildy, wykonane na zupełnie odmiennej dro- 
dze, mianowicie metodą Poincarć go wynajdywania rozwiązań 
okresowych, potwierdzają rezultaty otrzymane przez Brende Га 
jego własną metodą, co daje nam najlepsze świadectwo o war- 
tości tej metody, stosunkowo łatwej w zastosowaniach prak- 
tycznych. 


—— s—assoəsL 


') Brendel. A. N 3396. Buchholz H. Untersuchungen etc. 
2) Schwarzschild. A. N. Bd. 160. 
35 Dziewulski. А. N. Bd. 181. 
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dur le mouvement des petites planètes du type de Thule (2). 


Détermination des termes élémentaires et de la lacune 
pour le type T'hule. 


Le nombre croissant des petites planétes a suscité l'appari- 
tion des méthodes spéciales en mécanique céleste qui permettent 
de calculer plus ou moins rapidement les perturbations. 

Gyldén et ses nombreux ćleves ont divisć les petites pla- 
netes en groupes, appelés „types“ ou la relation de commensura- 
bilité entre le mouvement moyen diurne de la petite planéte et ce- 
lui de Jupiter joue un róle tres important. On sait que si ce rap- 
port peut sexprimer approximativement par un rapport simple 
des deux nombres entiers, les perturbations, causées par la planéte 
principale, deviennent de plus en plus grandes á mesure que le 
rapport de ces deux moyens mouvements diurnes s'approche de 
la commensurabilité exacte. Jupiter produit dans l’essaim de 
petites planetes en général des grandes perturbations gráce á sa 
masse énorme. Les perturbations Joviennes deviennent déja trés 


1 
sensibles pour les planćtoides du grouppe 3 (moyen mouvement 


(moyen mouvement diurne ca 


| 
diurne ca 900”) et du grouppe D 


600"). 

Si la commensurabilité devient trés approchée, les perturba- 
tions croissent trés vite et les méthodes ordinaires du calcul des 
perturbations deviennent insuffisantes ou trés difficiles et pénibles 
dans les applications, surtout s'il s'agit des planétoides-types. 


104 J. KRABSOWSKI. 


C'est à Gyldén que nous devons une des méthodes du 
calcul approché des pérturbations. Cette méthode, remaniée et 
perfectionnée par ses éminents éléves (M. M. Backlund, Bren- 
del, Harzer, Wolf etc.), peut étre méme appliquée dans les 
cas les plus difficiles oü la commensurabilité des moyens mouve- 
ments de la planété principale et de la petite planéte est à peu 
prés exacte. 

Cette méthode a fait ses preuves; elle a été déjà plusieurs 
fois appliquée, avec bon résultat, aux différents cas difficiles du 
calcul des perturbations des petites planétes appartenant aux grou- 
pes d'Hécube, de Hestia ou Hilda. 

M. M. Brendel et Kramer nous ont donné des tables 
étendues, qui nous facilitent notablement le calcul des perturba- 
tions dans le cas des types Hestia et Hécube et permettent rapide- 
ment de déterminer leurs positions géocentriques. Les calculs de 
M. Harzer ont été sur une autre voie brillamment controlés par 
M. Simonin pour Hecube. Le cas difficile du type de Hilda E | 
a été commencé d’être étudié par M. Buchholz, qui a déterminé 
les termes caractéristiques, sans donner l'expréssion méme des per- 
turbations. 

La petite planéte Thule (279) appartient au type ?/, des péti- 
tes planètes et jusqu'à présent c'est l'unique réprésentant de ce 
type. Thule peut s'approcher de l'orbite de Jupiter beaucoup 
plus que Hilda. M. Bidschof trouve que des grandes proximi- 
tés entre Jupiter et Thule arrivent jusqu'à trois fois dant l'éspace 
d'un siécle. Pendant un pareil rapprochement le mouvement de 
Thule est fortement perturbé par la puissante action de Jupiter, 
qui produit des variations trés notables dans son orbite. Ainsi 
pendant la proximité qui a eu lieu en 1912, les calculs de M. Vil. 
jev nous montrent que le moyen mouvement de Thule a subit de 
trés fortes variations. 

M. Kobb, en se basant sur les travaux de M. Darwin sur 
les orbits périodiques, a cru pouvoir admettre que l'orbite de Thu. 
le peut devenir instable et s entrecouper méme avec celle de Ju 
piter. 

Jusqu'à présent le mouvement de cette petite planéte, appar- 
tenant à un type planétoidal si intéréssant, n'a pas été étudié 
d'une maniére générale. A ce que Je sais, M. M. Bidschof et 
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Wedemeyer ont seulement calculé et corrigé les éléments de 
Thule. M.Viljev de son cóté n'a fait que le calcul des pertur- 
bations de son mouvement, limité à un court espace de temps en 
se basant sur les éléments de M. Wedemeyer. 

Il m'a semblé intéressant d'essayer d'appliquer à Thule les 
méthodes de Gyldén, modifiées et perfectionnées avec tant de 
succés par M. Brendel. L'orbite de départ de nos calculs a été 
celle de M. Wedemeyer. 

Dans le présent mémoire je donne un aperqu succint de la 
méthode de Gyldén, modifiée par M. Brendel, je détermine pour 
le type Thule les termes caractéristiques ainsi que la lacune (dans 
le sens de M. Brendel) dans l'anneau planétoidal, qui jusqu'à 
présent n'a pas été determinée pour ce type. 

Pour ne pas étre trop long je me permets de renvoyer le lec- 
teur aux formules contenues dans le texte polonais, aux quelles, 
du reste, je fais les renvois donnant les numéros des formules et 
les pages correspondantes du texte. | 


1. Dans le premier chapitre je me suis attaché d'exposer, 
aussi brièvement qu'il m'a été possible, les idées générales sur les- 
quelles repose la théorie de Gyldén ainsi que les équations les 
plus importantes de la théorie qui nous intéresse, 

On sait que c'est la méthode de Hansen qui a sérvi à Gyl- 
dén du point de départ. Le équations fondamentales sont: I pg. 
4; Il pg. 8; et III pg. 9, texte pol. Dans ces équations v désigne 
la longitude vraie de la planéte et 7 son rayon vecteur. 

Les équations de Hansen en coordonnées polaires, Gyl- 
dén transforme, en introduisant des nouvelles variables p, y S 
et v, qu'il appele ,longitude vraie“ et qui est une variable in- 
dépendante. Les équations nouvelles ainsi obtenues sont: IV, 
[VPis pg. 11, V, Vbis pg. 12. Nous introduisons maintenant les no- 
tations de M. Brendel!) et nous rappelons la définition des ter- 
mes élémentaires et caractéristiques à longue et courte 
période. Parmi ces termes nous distinguons, suivant M. Harzer’) 
ceux qui appartiennent au genre 4, B, C, D et ceux que M. Bren- 


l Brendel. Theorie d. kl. Planeten |. 
3 Harzer P. Untersuchungen etc.: pg. A 
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del appelle „rein erster Ordnung“. Avec M. Brendel nous 
designerons par les lettres latines les termes du premier ordre, et 
par les lettres grecques ceux qui contiennent, comme diviseur, la 
petite quantité 8. 

Dans le présent mémoire nous nous bornons à ne considérer 
que le mouvement de la petite planéte qui s'effectue dans le plan 
méme de l'orbite: c'est pour celà que nous laissons de cóté les per- 
turbations de l'inclinaison 

En nous servant des rélations connues, données par M. Bren- 
del, entre les quantités Gyldéniennes 7, II, R nous obtenons: 


e GELEET 
l+ncosv 


(1) (r) = 


ой (r) c'est le rayon vecteur „absolu*. 
Aprés quelques transformations on parvient à l'équation sui- 
vante: 


(dzi) 

l +T cos v 
Cette équation obtenue par Gylden est tres importante, 
car le second membre ne contient aucun élément à longue période, 


tandis que le premier membre de cette équation est formé par une 
Íonction pareille. 


Pour achever toutes les transformations nécessaires on intro- 
duit encore le „temps reduit“ (7) qui est donné par l'équation 
suivante: 


(3) n(t)-- A=v-- XB, sin nv ou 
L=nt--N 


=1+ En В, cos nv. 


(2) 


et la quantité & est definie par: 
dE d -Ba cos zl PE . ab, sin nll 
"o. - — Sin nv — A ——_ COS NV, 

dv dv dv 


alors (12) devient: 
du M— 9% , dE 


—— 


(4) da. ACOS ЕИО 


Si W nous donne la différence entre l'angle déterminé par le 
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temps vrai ní et celui qui correspond au „temps reduit“ 
n(t), nous aurons encore la relation suivante: 


(5) W —nt—n(t) 


qui enfin nous donnera l'équation différentielle VI pg. 20 ou ent- 
rent aussi les quantités R et S. 

Pour les applications de la méthode les relations IV, IVPis, V, 
Vbis et VI sont les plus importantes. Nous allons maintenant tâ- 
cher d'intégrer ces équations. Ces équations contiennent les quanti- 
tés P et Q qui sont des dérivées partielles par rapport à v et à r de 
la fonction perturbatrice Q que nous devons développer en série 
suivant les procédés, donnés par Gyldén. 

2. Onsait que Gyldén pose: 


a 


(6) A = h, + 2 R, cos H+-2R, cos 2 H Lt ... 
ou z | 

dt A POWO AZ 
© dek еа RO rear Sin? ç 


en substituant dans cette expression les coefficients numériques 
: a 
connus Tee, ne dépendant que de a (on se rappele que a= le 


rapport des demi-axes de la planète troublée et troublante) on 
trouve: 


a! Ke, | | 
(8) Jr, — q (1 19 A)? | 1n.0 Хал A + In.2 № — Tu.3 M. ( E 


Ce développement converge encore dans le cas de Thule. 
En introduisant dans (6) l'expression trouvée, nous obtenons 
| m "Da i 
le développement en série seulement de la partie ^ de la fonction 
perturbatrice, qui, comme on le sait d'ailleurs, étant complete, 
s'exprime au moyen de la relation suivante: 
a L 

(9) a (Q9) — m' A р = ZONY H | | 

Si nous dćveloppons cette expression en sćrie, en introdui- 
sant les coefficients numćriques connus, nous obtenons: 
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(10) a(Q)— 2m’ > Q (n, s, s"), v Lëlz 1? cos nH 
am f, 8,8, yv, y 
(11) 0, =m' >: Q(n.s.s') v.v psp! y» nv. 


na, 8 v. v! 


Ces coefficients sont des fonctions numériques de a et peu- 
vent s'exprimer au moyen des т. !) Cette relation est donnée pg. 
20 texte pol. (20). 

Les quantités P et Q peuvent s'exprimer, aprés quelques 
transformations, au moyen des relations suivantes: 


(12) | P=2 Pan (^ p^ n? 3?" cos n(v — v) 


Q —2 X Q inss’) vw D° p'a. uh na sin n (v — v) 


(oü * signifie que pour n=0 le facteur 2 doit étre omis). 

Les coefficients de ces expressions sont données par les for- 
mules pg. 27 — 30 texte pol., qui n'ont pas besoin d'aucune expli- 
cation. 

Dans les expressions (12) nous écrivons p en séparant la par- 
tie élémentaire du genre B et les termes caractéristiques: 


(13) p=(p) R 


oü (p) désigne les termes élémentaires du genre Bet P contient 
les termes ordinaires et caractéristiques. En ne considérant dans 
p que les termes caractéristiques seuls et les substituant dans P 
et Q, nous obtenons l'expression (39) pg. 31 texte pol. pour P et 
une expression analogue pour Q. 

Pour avoir les expressions définitives des quantités P et Q 
il est nécessaire d'effectuer encore une transformation, parce que 
dans ces quantités entrent les S, R, W, qui de leur cóté contiennent 
les termes des genres différents (ordinaires, à courte et longue pé- 
riodes, élémentaires etc.). En effectuant des transformations ana- 
logues aux précédentes (pg. 32 texte pol.) on obtient les expres- 
sions définitives pour P et Q (pg. 32, 33 texte pol.) 2). 


') Brendel, Th. d. kl. Plan. I, pg. 92; Harzer, Untersuch., pg. 49; 
Kramer I, Th. d. kl. Р]. pg. 9. 
° Brendel, c. pg. 70; Harzer: 1. c. pg. 30. 
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3. Nous allons déterminer maintenant les termes caractéri- 
stiques et élémentaires dans les cas de Thule (279). 

Aux pg. 12 et 14 texte pol. nous avons donné les équations 
différentielles qui déterminent les quantités 2 et p; l'équation dif- 
férentielle en W leur est analogue: 


(haj: W=0 + VW. sin (A 0— В,). 


Dans toutes ces équations a„, da, c, sont des petites quan- 
tités, A, est un petit diviseur. Quelquefois cependant, aprés avoir 
effectué l'intégration, il arrive que malgré que les quantités a,, Dn, 
c, sont petites, les termes renfermant le diviseur A, peuvent être 
ou peu différents de zéro ou trés grands!), car, aprés cette opéra- 
tion, les termes considérés peuvent avoir comme facteur ~ ou 

"m 
e Si A, est peu différent de zéro, l'argument de la fonction 
cosinus est un tres petit angle, ce qui nous donne pour le terme 
considéré une période tres longue. 

Si A, diffère peu de 1 alors la période du terme, renfermant 
un pareil A est peu différente de la période de la variation de v, 
par conséquent elle differe peu de la durée de révolution de la pla- 
néte considérée. Des pareilles expressions nous donneront les ter- 
mes á courte période. 

Par l'intégration, dans les expressions se rapportant à 5, se- 
ront agrandis seulement les termes à longue période, tandis que 
ceux à courte période vont démeurer sans changement; inverse- 
ment dans les expressions de p aprés l'intégration augmenteront 
seulement les termes à courte période et ceux qui auront une pé- 
riode longue ne seront pas changés. 

M. Brendel a remarqué, que dans le cas d'une commen- 
surabilité trés approchée entre les moyens mouvements de la pla- 
néte troublée et troublante, il arrive que les termes du second ordre 
sont encore trés importants. L'idée de Gyldén répose sur ce 
principe que dans les expressions approchées, obtenues pour R 
et Q, nous ordonnons les termes non suivant leur ordre ou degré 
1) Brendel, Om amvändningen etc., pg.6. Ludendorff, Die Jupi- 
terstörungen d. kl. PI. von Hecuba-Typus, pg. 4. 
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mais selon leur valeur véritable, qu'ils obtiennent aprés l'intégra- 
tion. C'est pour celà que nous devons d'abord pour Thule déter- 
miner dans les développements de S, p, R. W les termes à longue 
et courte période. 

Nous commencerons par séparer avec M. Brendel dans 
S, р, W, R les parties à longue et à courte période ainsi que les 
termes ordinaires, en souscrivant sous chaque lettre un des indices 
d, К, e. Nous aurons donc ainsi 


S = Su-F- S, S, R= Rit Е, + R. 
p=(p) 4 R W-(v4-W.-4-K 
K=W,+W 


où 


et yv c’est le terme séculaire. 

Dans notre cas spécial le rapport du moyen nouvement 
diurne n' de Jupiter au méme mouvement n, de la planéte Thule 
(279) est: 

п! 3 — 0, 
"nes Ba == 0.7398 = 17 


(б, représente une petite quantité) d'où : 


Généralement l'argumeut des fonctions trigonométriques qui 
entrent dans nos développements a la forme suivante: 


nw=n(l —p,) v—n B—ny Wi. 


Donc, dans les termes à courte période, le coefficient de v sera 
peu différent de 1 et dans ceux à longue période ce même argu- 
ment différera peu de 0. Maintenant, dans les expressions préce- 
dentes pour P et Q, il faut déterminer les termes de ces genres dif- 
férents; les expressions de ces deux quantités, données aux pg. 42— 
43 texte pol. nous donnent tous les termes particuliers qui cara- 
ctérisent le type planétoïdal considéré. Nous devons encore trans- 
former ces expressions pour pouvoir calculer les différentes par- 
ties mentionnées plus haut, des expressions de S, k et W. 
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Comme on a remarquć précédemment, ces quantités sont 
formées avec des termes à longues et à courtes périodes, qui d’ail- 
leurs dépendent des mémes arguments que dans les expressions 
pour P et Q. 


En introduisant la classification de M. Harzer des diffe- 
rents termes particuliers et faisant des transformations nécessaires, 
nous obtenons les expressions (10 pg. 40 texte pol.) oü les diffé- 
rents termes de chaque catégorie de M. Harzer sont spécifiés 
par la majuscule correspondante. 


Dans (10) pag. 46 texte pol. cinq expréssions ne contiennent 
point la quantité à, qui, comme on le sait, est seule liée avec le 
mouvement d'une petite planéte déterminée, ces expressions se- 
ront donc appliquables au mouvement d'une petite planéte quel- 
conque!), par contre les termes contenus dans (10) 1. c. doivent 
étre calculés à nouveau dans chaque cas traité spécialement. 
C'est pour celà que Gyldén a appelé ces termes caractéri- 
stiques, comme ceux qui nous déterminent dans chaque cas 
spécial les conditions du mouvement d'une certaine petite planète. 
Dans les termes de la catégorie (C) et (D) nous trouvons dans 
l'argument la quantité W;, qui est une fonction variant peu avec 
le temps; c'est d'ailleurs la partie à longue periode de W, qui est, 
suivant la nomenclature de M. Brendel, „la reduction du 
temps". Cette quantité nous rend difficile l'intégration de nos 
équations différentielles. Dans les termes de la catégorie (4) l'in- 
tégration introduit un petit diviseur 0, de l'ordre de la masse per- 
turbatrice. Dans (D), (C) nous avons ó en dénominateur après 
avoir integré l'expression pour S, qui comme on le sait, contient 
la masse avec l'exposant moindre que l'unité. Ces termes devien- 
nent zéro avec la masse. Nous devons aussi remarquer que l'in- 
tégration de l'expression pour S fait augmenter les termes du genre 
(A) et (C) tandis que dans p, sont augmentés seulement les terines 
de la catégorie (D) et (D). Nous pourrons donc constater que Sz 


, 


sera de l'ordre de m’ et $9, Ri. Ra, Wa, Wa, K,— de l'ordre ^ 


et que les termes ordinaires qui entrent dans S, R, K seront de 
l'ordre de la masse perturbatrice. 


' Brendel. Theorie d. kl. Plan I, pg. 106. 


112 J. KRASSOWSKI. 


Nous trouverons donc les expressions definitives des diffe- 
rentes parties de 5, E, K, qui sont données aux pages 48 et 49 du 
texte pol. Nous devons encore substituer les series trouvćes 
dans les expressions pour R et Q. Les expressions ainsi obtenues 
sont données aux pages 50—60 du texte pol. G, désigne la somme 
des termes ordinaires; Pa, qu et ga sont des coefficients numéri- 
ques qu'on peut calculer sans aucune difficulté. 

4. Nous allons effectuer maintenant l'intégration des équa- 
tions différentielles données, (IV pg. 11, V pg. 12, VI pg. 20 texte 
pol.). Plus haut nous avons déjà séparé les différents terines sin 
guliers qui entrent dans les expressions des mouvements des peti- 
tes planétes. M. M. Brendel et Kramer ont divisé les petites 
planétes suivant leur moyen mouvement diurne en certaines clas- 
ses. Suivant cette classification la planéte Thule appartient à la 
premiére classe!). Le mouvement des petites planétes apparte- 
nant à cette classe est caractérisé par des termes de degré 0, qui 
entrent dans À et K, qui sont grands par rapport à la masse de 
Jupiter. Dans cette classe on aura aussi dans le développement 
de À un terme important qui, comme l'a démontré M. Brendel, 
appartiendra à la catégorie D de M. Harzer c. à d. caractéri- 
stique à courte période. 

Comme Thule peut se trouver en grande proximité avec Ju- 
piter, les séries développées suivant le rapport des axes ou des 
rayons vecteurs, deviennent lentement convergentes, il faut donc 
dans les limites d'exactitude?) que nous nous sommes assignées, 
dans les développements des coordonnées et dans les équat. IV, 
V, VI, omettre la troisieme puissance de la masse de Jupiter, ainsi 
que dans, les expressions du degré O, les quantités qui seront du 
troisieme ordre par rapport A la masse perturbatrice, divisée par 
une petite quantité s, parce qu'elles ont encore en facteur une 
petite quantité 4 ou n’ ou les carrés de ces quantités. Comme 
nous poussons nos dćveloppements jusqu'a la troisiéme puissance 
de la masse perturbatrice, nous devons A l’exemple de M. Kra- 
mer (dans le cas de Hécube) compléter les développements de 
M. Brendel. Nous obtenons ainsi les équations (1) pg. 63, (2), 


') Brendel, Th. d. kl. Pl. I pg. 92; Kramer I, pg. 20, 144. 

2) L'exactitude que nous nous sommes assignés dans la reprćsenta- 
tion du mouvement de Thule cst celle que M. Brendel admet dans sa 
„Theor. d. kl. Plan.*. 
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(3), (4) pg. 64. (5), (6) pg. 65, (9) pg. 66. (12), (13) pg. 67 et les re- 
lations (14)—(18) pgs, 68—70. Nous avons donc trouvé, dans les 
limites d'éxactitude désirée, tous les coéfficients de P et Q. Ces 
coéfficients sont des grandeurs numériques que l'on peut calcu- 
ler sans beaucoup de difficultés. 

Maintenant nous pouvons aborder l'intégration en se ser- 
vant de la méthode de M. Brendel, que M. M. Ludendorff, 
Buchholz et Kramer ont appliqué avec tant de succés. Nous 


: 45 
commencerons par intégrer l'équation qui nous donne de En se 


limitant aux termes du degré 0 nous aurons = Q——4qsindv 
et comme w=(l—u,)v— B — y Wa, contient aussi des termes 
périodiques d'un ordre supérieur à zéro (nous savons aussi 
que le terme séculaire lui appartient) nous trouverons donc après 
l'intégration de S des termes qui seront au moins de l'ordre 1, 
parceque dans $, il ne peut être des termes du genre A et C. 

En intégrant (1) nous trouvons: (2) S,=a,--a, cos 4w--S,.00 
(S,oa désigne les termes ordinaires), aprés quelques transforma- 
tions nous obtenons l'équation (2) pg. 73 (texte pol.). 

(3) So = 4, + (9 +0", + 0" 8,2) cos 420 + Snoo- 

Dans (3) nous pouvons calculer tous les termes, sauf a, que 
nous déterminerons ultérieurement. 

Intégrons maintenant l'équation différentielle qui détermine 
la coordonnée À. Comme nous omettons pour le moment les 
quantités du premier ordre, nous devons remarquer Ou alors (p)—90, 
parceque cette quantité contient seulement les termes à courte pé- 
riode du genre 5, qui sont au moins du premier degré. L’equa- 
tion (1) pg. 63 texte pol. va nous donner !): 

da? Т 4 OR 
(4) jr | R = 25,—P, — 0 zo car p — (p) +R. 


Nous pouvons écrire aussi: 


P, = p, -+ p, cos 410. 


1) Les indices O indiquent que les fonctions correspondantes sont de 
l'ordre O. 
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S, nous est connue par (3), il faut donc calculer Q, ES - | 
0 

l'équation 6 pag. 34 nous donne l'équation (25) pag. 73 texte pol. 

d'oü, aprés quelques transformations nous obtenons l'équation sui- 

vante: 


а R 
(9) "E + R-2a, + p' + {Dh 4- (9), ё, (Ds + PB? + 
+ (D) But, + 0D, 1 cos 4 ww 


où les constantes autres que a, sont données par (33) et (34) pg. 


74 —75 texte pol. 
L'équation (12) pg. 48 texte pol. nous donne: 


d? P 4 
(6) dio + R=b,—[1—(1--0,)*]8, cos 4w, 
= b, — (28, 40,5) f, cos 4w 
ou 


b, =2a, De et 20, 4-0, =— IM), + Mh HCD) Ti + 
+ (0), 8,2 + O95 But + (le gd 
Au lieu de y, écrivons 
_ 28 
| 1--8,' 
| alors apr&s quelques transformations nous trouvons: 
= шым d Ai pp 
(7) T 24,4 Po PR Z 28, +- Е —1— cos 4w, 
où les coefficients sont donnés par les (29) (37) — (40) pg. 66 tex. p. | 
Nous avons donc déterminé tous les termes singuliers de la | 
fonction £; nous devons encore calculer entre les termes ordinai- 
res ceux qui seront grands. Nous pouvons écrire: | 
d? R 
(8) duż + R = ZW, oo COS NW. | 
ой Чоо == 258400 — Dann, 
En appliquant à l'intégration de l'équation (8) la méthode de | 
| M. Brendel, nous obtenons: 
(9) Р, 0; sin V — gą cos 0 
oü les g, et g, sont donnćs par les relations (43) pg. 78 tex. pol. 
Apres avoir ćffectuć la substitution dans (9), des valeurs trou- 
vees pour les constantes, nous arrivons A la relation (44) pg. 78 
tex. pol. qui est: 


| | 
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ЖӨН Е m p; А 
(9) Боо = — e 
1 


pen à) y 
Um 


а 


où, pour n + 4, 8, оо est donné par (21) pg. 71 tex. pol. 


ll nous reste encore à intégrer l'équation différentielle qui 
nous donne la relation entre W et v. En nous limitant aux termes 
de l'ordre 0 nous aurons: 


(10) 24 = Š, — 2R, —2R,8, + ЗР, + ARA — ARS. 
et les relations (47) pg. 79 tex. pol. En substituant ces relations 
dans l’equation (10) nous trouvons enfin: 

dW 
(11) T =c + + (a, — 28, — 28, F 6058, + 

2 a, 8,2 — 38,9) cos 4w — (G, b, —3 8,4 
+ $ ao Bi — 6 bo B1?) cos 8w + ($a4,6,* — 81°) cos 12w, : 

+= B, et c= ao — A,B, + š 0,8, — 60,5," — 20,. 
c, contient m’ en facteur; AB, contient m’ divisé par la petite quan- 
tité ð}. 

Nous allons intégrer (11) en se rappelant que 


(12) W=w+W.+K, 


où 


Mais y c'est la partie séculaire, donc dans W; et a fortiori 
dans K les termes à longue période du degré O n'existant pas, 
nous aurons: 


AW —.dK 
GR dy TT dv 
En comparant (11) et (13) nous trouvons 
(14) YZ sas 


En ayant ćgard a la relation (13) pg. 49 tex. pol. et differen- 
tiant (10) pg. 49 tex. pol., nous trouverons l'équation (49) pg. 80 
tex. pol. 

En substituant dans (13) cette derniére équation, on trouve 
l'équation 50 pg. 80, d'où nous pourrons conclure les relations 
entre les coéfficients qui y entrent. Puisque Thule suivant M. 
Brendel appartient aux petites planétes critiques de la premiere 
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classe, il est évident que nous aurons c, =0 et y==0. Dans notre 
cas spécial 4 — y D et l'équation (13) deviendra 
(14) W=yv-+-K. 


Si nous connaissons dćja la constante ge, nous pouvons, en 
se rapportant aux formules (51) et (53) pg. 81 tex. pol., calculer les 
membres singuliers de l'équation (14), en se limitant toujours aux 
termes du degré zéro. 

En éffectuant maintenant sur la fonction W les mémes trans- 
formations que nous avons appliquées à R, nous obtiendrons sans 
peine. 


(15) W, — 9509 — 2 Bano 


i (1--0,) 

Déterminons à présent les constantes d'intégration. Ces con- 
stantes se reduiront au nombre de quatre, puisque nous ćtudions 
pour le moment le mouvement de la petite planéte dans le plan 
méme de son orbite. Les constantes à déterminer seront: k le mo- 
dule de l'excentricité,  — la longitude du périhélie, a — le demi- 
axe de l'orbite et n la constante du mouvement; a, joue dans nos 
formules le róle de constante surnuméraire, que nous pouvons 
pour celà traiter comme une constante arbitraire. 

La constante a, nous déterminons au moyen de d,, en sup- 
posant que dans le développement de R la partie constante doit 


; (nz 1). 


etre effectivement du premier ordre (rein erster Ordnung). La fon- 


ction W = И.К et K est formée au moyen des fonctions À et S, 
pour celà la partie constante de W se composera des parties con- 
stantes de À et de S. Nous aurons donc: 
(16) сот — 20. 

Mais M. Brendel?) a démontré, que dans le cas d'une pe- 
tite planéte caractéristique nous avons la relation 

Clo Z 2b, = () = CC. +7. 

Si la planćte considérée est critique (le cas de Thule) y ne 

peut étre zero et dans ce cas l'équation suivante: 


a, — 20, =c ++ 


1) Brendel, l. c. pg. 92, 100. 
2 Brendel, l. c. kap. V, sparsim. 


I 
I 


- - 
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a lieu aussi— c'est pourquoi nous trouvons: 
(17) TEN 


On voit donc que dans le cas qui nous interesse, l'équation 
53 pg. 81 texte pol. et l'équat. (17) vont nous déterminer ła con- 
stante Ou 


La fonction W nous pouvons écrire: 
(18) W = (coH Y+ 19) v + les termes périodiques 
où (= T+ o 


mais x, est au moins du second degré, alors y = y. 
On a trouvć la relation suivante: 


w, =(1—p)v—B—uW, 


et 

(19) 10, = 10 — y K 

donc v=w +yK=(1—u)v — B—wW 4-yK, 
mais WC Wan 


et w=(1 —- u.) -|- les termes period. [w =p(l--7)]. 
De ces relations nous trouvons ainsi la partie sćculaire de W: 
pars sec w= (1 -|- pa) v. 


Rappellons— nous le définition des quantités y et р... 
On avait les relations suivantes: 


=p; et ==. 


(n c'est une constante d'intégration „la constante du mouvement" 
selon M. Brendel, n, c'est le mouvement diurne vrai). Dans 
les relations précédentes introduisons à et à, (qui sont des petites 
quantités) alors nous trouvons: 


On voit que dans nos expressions, dont nous avons besoin, in- 
tervient maintenant 0, seul, qui dans le cas méme de commensu- 
rabilité exacte est toujours beaucoup plus grand que la masse per- 
turbatrice, tandis que la quantité 2, dans le cas d'une commensu- 
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rabilité exacte, peut devenir 0. L'idée de M. Brendel est que 
dans les expressions des perturbations apparaissent, comme petits 
diviseurs, non les ë mais le à. M. Brendel a démontré aussi 
que le rapport de commensurabilitć exacte peut avoir lieu seule- 
ment entre n еї n’ et non entre n, et n’. Dans ce cas aussi nous 
aurons ô= Q. 

Nous aurons maintenant les relations suivantes !): 


py = (T + Co +1) 
6, — 0,0--04-[-6 — 30,— 37. 


Comme nous étudions le mouvement d'une petite planéte 
critique pour laquelle: 


et 


Со) рышы 
nous aurons: 
ó =ô — 4 y. 


D'une maniére tout à fait analogue nous pouvons écrire: 


/ 
Ar = flo = - 5 et à, 0, — 41%. 

Quand 9, devient trés petit ou presque zéro, nous aurons 
alors à faire à la libration qui, dans la méthode de Gyldén nest 
pas difficile à étudier. La quantité ną nous donne le mouvement 
moyen diurne dans le cas de la libration; nous devons aussi con- 
naitre exactement cette quantité pour pouvoir déterminer la pé- 
riode de révolution de la petite planéte. 

Nous devons remarquer encore que p, et p, se distinguent 
seulement par la présence de la constante у. Cette quantité est 
de second degré et comme telle pour le moment elle n'entre pas 
dans nos calculs; <= est toujours positive et, dans le cas des 
planétes caractéristiques de premiére classe (auxquelles appar- 
tient Thule), étant méme du degré zéro, elle doit étre conservée 
dans nos calculs. 

9. Passons maintenant aux applications numériques de nos 
formules. La petite planéte Thule (279) fut découverte par M. Pa- 
lisa à Vienne le 25 octobre 1888; au moment de la découverte 
elle était de grandeur 13.5. Les premiers éléments calculés par 


— 
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М. M. Lange et Halm ont montré que cette nouvelle planète 
est un astre trés lointain, qui se meut à peu prés à la distance de 4 
unités astronomiques du Soleil. Une orbite plus éxacte a été cal- 
culée à nouveau par M. Bidschoff!), qui a utilisé dans ses cal- 
culs toutes les observations de Thule de 25/X 1888 à 10JII 1891. 
M. Btdschof, dans son travail, trouve que Thule se rapproche 
beaucoup de Jupiter (en moyenne trois fois pendant un siecle). 


Les éléments calculés par M. Bildschof sont donnés pg. 86 
texte pol. 

Aprés M. Bidschof, M. Wedemeyer?) à calculé une nou- 
velle orbite, en utilisant toutes les observations et en leur appli- 
quant les perturbations, calculées par la méthode d'Oppolzer. 
Les éléments de M. Wedemeyer se trouvent pg. 87 texte pol. 
Les observations de Thule, effectuées aprés la correction de son 
orbite par M. Wedemeyer, sont trés peu nombreuses, notam- 
ment en 1906 Thule füt observée le 7/X Vienne et en 1911 le 20 
et 30 avril photographiée à Heidelberg. 

La petite planéte s'approchait de plus en plus de Jupiter, ses 
perturbations augmentaient, et sans doute l'accord entre les po- 
sitions de Thule effectives et calculées devenait de moins en 
moins satisfaisant. La, proximité avec Jupiter arrivait en 1912. 
M. Viljev?) à effectué à nouveau les calculs des perturbations 
et les a appliqué aux éléments de Thule de M. Wedemeyer; il 
a donné ainsi des nouveaux éléments osculateurs (pg. 88 tex. pol.) 
et une éphémeride, déduite de ces éléments. Thule ne fut point 
retrouvée en 1913. 

Pour nos calculs nous avons utilisé les élements de M. We- 
demeyer (B. J. für 1914) et les éléments de Jupiter donné par la 
„Conn. des Temps“ pour 1914, en prenant comme valeur de la 


masse ml = —— . Dans ce travail nous nous limitons pour 
le moment à la recherche de la lacune (dans le sens de M. Bren- 
ү 3 
del) pour les petites planetes du type 4” c'est pour cela que nous 
donnerons ici seulement les calculs qui s'y rattachent. Dans le 8 2 
') Sitzungsberichte d. Ak. d. Wiss. Wien, Bd. 100, pg. 937. 


2 Wedemeyer, Archiv. d. deutschen Seewarte XXXI Jhgg. 
3) Astr. Nachr. 4661. 
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nous avons donné, conformément à la théorie, que nous appli- 
quons à notre cas spécial, les expressions analytiques du dćvelop- 
pement de la fonction perturbatrice. Dans ces développements 
entre la quantité a, qui s'exprime au moyen du rapport du demi- 
axe de l'orbite de Thule à celui de Jupiter. Les coéfficients nu- 
mériques de ces développements sont exprimés au moyen des quan- 
tités В et y. 

Calculons ces coéfficients. Avec les éléments de Jupiter et 
de Thule nous trouvons: 


log « — 9.9156572. 


Comme cette quantité differe peu de l'unité, nous voyons 
que les coéfficients 8, y ne se trouvent point ni dans les tables de 
M. Masal?) ni dans celle de Gyldén?). 

Nous devons donc les calculer spécialement en se basant sur 
leurs expressions analytiques. Nous calculons d'abord B, au mo- 
yen des formules de la moyenne arithinético-geométrique; nous 
trouvons: 

log 3%) = 0.1117634. 


Aprés avoir calculé ß,, nous calculons les autres coéfficients 
3, pour m =1, 2...12. Ces coéfficients nous obtiendrons le plus 
facilement en employant la méthode de M. Ludendorff?). 


La table des valeurs auxiliaires, donnée par M. Luden- 
dorff fut augmentée et au moyen de cette table ainsi transtormée, 
jai obtenu les f,. Comme contrôle de nos calculs, j'ai obtenu 
les mémes coéfficients en se servant des quadratures mécaniques. 

L'accord des valeurs des coéfficients ainsi obtenus à été com- 
plet. Les coéfficients 83,4, (5,5, Dia, @°,, ont été calculés au mo- 
yen des quadratures mécani ques, et controlés par les tables de M. 
Ludendorff. Pour les coéfficients qui nous restent encore à de- 
terminer, j'ai appliqué la relation suivante: 


84-2 _ (8 958-F2 
Ps TE Br Baa | 


!') Masal H. Astronomiska iakttagelser och undersökningar anstälda pa 
Stockholms Observatorium. Fjerde Bandet. Stockholm, 1891. 

2) Gyldén Н. Astr. Ges. Publ. 21. 

з) Ludendorff H. A. N. Bd. 160. 
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Les résultats obtenus furent controlés aussi par la quadrature 
mécanique. Ainsi nous avons obtenu les logarithmes de tous les 
„ qui sont données pg. 93 texte pol. 


Maintenant en appliquant la relation: 


*. 1 o De 5... (20 —!1) an 2o EF" 1 
T DIARIO n--1 


nous avons obtenu les logarithmes de 1,,, qui sont donnés pg. 94 
tex. pol. Nous pouvons déjà calculer les A et B (avec des indices 
différents) qui entrent dans nos formules. Ce calcul peut étre 
effectué de deux maniéres différentes. On peut appliquer les re- 
lations (20) et (24) pg. 27 et 29 tex. pol. et faire lecalcul de A et B 
au moyen des quantités Q, nx qui sont des fonctions de y... et des 
Pamp EL emp, lOnctions а ано ао - 

Cette facon de procéder est recommandée par M. Brendel, 
elle est trés pratique, quand on doit calculer tous les coéfficients 
nécessaires pour le calcul des perturbations, mais dans notre cas 
spécial où nous nous limitons aux termes du degré zéro, c'est 
beaucoup plus pratique de calculer les A et B directement en fonc- 
tion de 14, en se servant des tables trés pratiques de M. Masal’). 

Le résultat de notre calcul est donné pg. 95 tex. pol. 

Nous avons donc trouvé tous les coéfficients nécessaires 
dans le développement de la fonction perturbatrice, nous pouvons 
maintenant calculer pour Thule la lacune correspondante dans 
l'anneau planétoidal. 

Nous avons eu l'équation suivante, où p, est l'inconnue: 


(20, --0,?)8, — p! +p"8, +p""p,? 


(ou 0, s'exprime au moyen de f, (67) pg. 84 texte pol.). 

M. Brendel a montré que pour trouver les limites de la la- 
cune il faut déterminer 8, de cette maniére, qu'en assignant une 
certaine valeur pour 6, cette équation nous donne seulement deux 
différentes valeurs pour 8,, qui correspondent aux deux valeurs 
de ß,. Ces valeurs particulieres determinent les limites de la la- 
cune. En se basant sur la théorie de M. Brendel, nous pouvons 
déterminer ó, de telle manière que cette fonction devienne discon- 


1) Masal H. Formeln und Taícln zur Berechnung der absoluten Stö- 
rungen der Planeten. Kongl. Sv. Vet. Ak. Handl. Bd. 23. 
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tinue, alors suivant M Brendel A ces valeurs de 8, correspon- 
dent les ß,, qui déterminent la lacune. 


| ^ п! н í 
Puisque à, ne peut être zéro donc p=, пе peut étre égale 
1 


a une fraction simple, tandis que à peut être égale à zéro, quand р 
devient une fraction simple. 

C'est pourquoi dans toutes nos formules où nous trouvons 
en dénominateur ö,, nous ne pouvons le remplacer par à, méme 
dans le cas lorsque ces deux quantités différent trés peu l'une de 
l'autre 

Nous avons maintenant à calculer g,. Siö, ne peut être de 
l'ordre 0, alors 0,2}, est d'un ordre plus élévé que les autres ter- 
mes qui se trouvent dans notre formule. Au lieu de à, écrivons 


ô — 4p y =à — 68,7, 
nous aurons: 


208, — 125. =р' р", -- p" p,? 


qb LÉI 20 zen) ane 
РЗ, EI BF 19, Br ani" = (), 


ou bien: 
Pr +28, JB, tk=0  !) 
OU 
aS HN ERIS ne EL 
WTYK HOT A TD! 


Résolvons l'équation (1) en admettant qu'elle doit avoir une 
racine double; nous trouvons pour f, les valeurs suivantes: 


log(8,) — 8.759360 i log(B,)u — 91542123, 


Les deux valeurs de ß, ainsi obtenues sont justement les va- 
leurs critiques. Nous pouvons aussi sans difficultć calculer 6. 
Nous trouvons: 

log 0 — 8.559884 ; 


Connaissant (äi et (В,)и nous obtiendrons deux valeurs pour 
y en appliquant la formule suivante, donnée par M. Brendel": 


1 


ee ln 
жай PTE 
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Les valeurs de 1 sont: 
log yı = 7.696539; log yı = 8.495577 


et les valeurs correspondantes de 6,. 
log(ö,)ı = 8.334772; 1ор(ё,) = 8.750883, 


Maintenant de la relation 


nous trouvons les deux valeurs du moyen mouvement diurne vrai, 
qui limitent la lacune pour le type à. Ces valeurs sont: 


(ni) — 401". 732; | (nj) — 391". 485 
La constante n pour les valeurs critiques 8, sera: 


_ ¿de 
кас зон 
d'où: 
n — 403". 723. 

Les résultats, que nous avons trouvés peuvent s'exprimer de 
la maniére suivante: 

Pour les petites planétes du type 3 la quan- 
tité B, doit satisfaire aux conditions suivantes: 


log В, <8.759360 ou log f,>9.154128,, 


et 
log 2, <8.334772 ou log 0, 28.750883, 


d'où le vrai mouvement diurnen, pour le type Thule 
doit ćtre 
n, 401". 732 ou т, < 391". 485. 


Entre ces deux valeurs. de n, se trouve la lacune, dans le 
sens de M. Brendel, dans l'anneau planétoidal; une petite pla- 
néte avec un moyen mouvement diurne vrai (comme il est défini 
par M. Brendel) contenu entre ces limites de n, aurait une or- 
bite instable. 

(Il taut remarquer qu'on ne peut comparer n, immédiatement 
avec le moyen mouvement diurne osculateur, calculé avec les élé- 
ments de l'orbite de la petite planéte). 


124 J. KRASSOWSKIT 


Nous avons donc obt enu les valeurs critiques du vrai moyen 
mouvement diurne qui nous dćterminent la lacune dans le type 2. 
Cette lacune jusqu'à présent n'a point été déterminée. 

Cet élément ne fut déterminé par M. Brendel que dans le 
cas d'Hécube, par M. Buchholz dans le cas de Hilda et par 
M. M. Schwarzschild et Dziewulski dans le cas de Hécube 
et Hilda. Les deux derniers auteurs employérent la méthode de 
Poincaré et recherchaient les solutions périodiques; les résul- 
tats quils ont obtenus, sont en parfait accord avec ceux de M. 
Brendel et confirment ainsi la haute valeur de la nouvelle mé- 
thode. 


‘À 
A 
Pa | 


DO NABYCIA WE WSZYSTKICH KSIĘGARNIACH 


NASTĘPUJĄCE DZIEŁA, 


+ 

et 

* z WYDANE 7 ZAPOMOGI 
KASY POMOCY DLA OSÓB PRACUJĄCYCH NA POLU NAUKOWEM 

imienia D-RA MED. JÓZEFA MIANOWSKIEGO, 


LUB OFIAROWANE NA RZECZ KASY. 


ch, 
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Wł. Wojtowicz 1916, XV + 492. . . . aL-— 
Kowalczyk. O sposobach obliczania przeszkód biegu ciał niebieskich 
. 7 


1901, XII, 624 + 3 tabl. 50 
Loria Gino. Przeszłość i stan obecny najważniejszych teoryi geome- 

trycznych. Przekt. uzup. dod. S. Dicksteina 1889, VII + 112. . — 60 
Majewski W. Geometrya praktyczna. Podręcznik dla rzemieślników. 

1903, 301 (karton; . . t 75 
Malanowicz J. Kreślenie geometryczne i jego praktyczne zastosowania. 

1907, XI 4- 176 + 40 tabl. i 346 rys. (karton) . — 60 
—  Kreslenie i zdobienie geometryczne. Kurs pierwszy, wyd. P» 1912, 

УШ + 110 + II + 45 tab. i 225 rys. (karton) . . . — 60 


b. — Rzuty geometryczne. (Geometrya wykreślna dostosowana do po- 
F trzeb przemysłu i rzemiosł}. Podręcznik dla szkół techn. i zawod. 

| ОССЕ А esas — 75 
Pieri Mario. Geometrya elementarna oparta. na pojęciach . „Punktu i Kuliś 

z oryginału włoskiego przełożył Stefan Kwietniewski 1915, 152 . 1 — 
h Poradnik dla samouków. Wskazówki metodyczne dla studyujących 
poszczególne nauki Wydawaictwo A. Heflicha i St. Michalskie- 
go. Wydanie nowe. Tom 1-52у w opracowaniu J. Łukasiewicza, 
| Z. Janiszewskiego, St. Kwietniewskiego, St. Mazurkiewicza, Wł. 
8 Sierpińskiego i St. Zaremby 1915, XXXIX + 618, z 314 fig. w tek- 

B ście i I tablicą. . . ro A 2:40 
d Prace matematyczno - fizyczne, wydawane przy współudziale Wład. 


^ Natansona, J. Puzyny, M. Smoluchowskiego, S. Zaremby, K. Żó- 
HN rawskiego przez S. Dicksteina. Tom XXVII, 1916, 277, rys. 10 . — — 


Rozmarynowicz T. Matematyczne podstawy ubezpieczeń na wypadek 
niczdolności do pracy w zastosowaniu do urządzenia kas emeryt. 
Wyd. Bolesław Danielewicz, 1886, IV -- 53. . . 

Schiir F. Podręcznik geometryi analitycznej. Przełoż, T Lopuszañski 
1901, X + 246 + 1 nlb., rys 76. . . 

Sierpiński Wł. dr. Teorya liczb niewymiernych. Wykłady uniwers. 
1910, 4 nlb + 149. . ua 

— Teorya liczb (kurs uniwersytecki) 1914, XIV R 412 . 

Sochocki. Rozwiązywanie równań liczebnych. 1884, XII + 212 

Sprawozdania z piśmiennictwa naukowego polskiego w dziedzinie nauk 
matematyczno - przyrodniczych: 

I ~ (1882, 1883, III + 187, 
II — (1883) 1885, III -+- 222, 
HI — (1884) 1886, 294 -+ 10, 
IV — (1885) 1887, 300 -+- VI, każdy tom 

Steiner J. Konstrukcye geometryczne wykonane za pomocą linii prostej 
i stałego koła; przełożył Stefan Kwietniewski 1915, VIII + 69, 
rysunk. 25 

Szczepański Józef. Stopień wyższy matematyki elementarnej 1 poczat- 
ki rachunku nieskończonostkowego. Podręcznik dla klas WMASSYCH 
szkół średnich i dla samouków. 1906, 450 

Wektor. Czasopismo matem. fiz. Red. Władysław Wojtowicz w War- 
szawie. 5, — z przesyłką 6, — zeszyt pojedyńczy . 

Witwiński Romuald. O układach odwracalnych powierzchni potrójnie 
ortogonalnych. Odbitka z tomu XXVIl-go Prac mat. plizy CZy ilas 
1915492, гуѕ.. Оз 

Zaremba Stanisław. Wstęp do analizy. Część [- sza. ` Pojęcie defen 
matematycznego oraz inne wiadomości pomocnicze 1915, 124 . 


NAUKI PRZYRODNICZE. 

Berdau Feliks dr. Flora Tatr, Pienin i Beskidu Zachodniego, 1890, VI + 
827 + 55. 

Braun Juljan. Badania ` w dziedzinie azotowych związków organicznych 
i ich pochodnych (1900 — 1908), 1908, VII + 238 . 

Chmielewski Z. Podręcznik analizy chemiczno-rolniczej 1905, 169 . 

Dyakowski B. Zarys metodyki elementarnego kursu historyi naturalnej. 
Wyd. W. Jezierski. 1909, 38. 

Dziedzicki H. Atlas organów rozrodczych (hypopyglum) typów Winnertz'a 
i gatunków znajdujących A w jego zbiorze. Mycetophilidów, 
1915, 16, tabl. XXI . 

Dzieje myśli. Tom I zesz. I. o rozwoju metod badań naukowych. Wie- 
dza ludów pierwotnych. Dzieje astronomii. Rys rozwoju fizyki. 
W opracow. Wł. Heinricha, Ludwika Krzywickiego. Stanisława 
Kramsztyka i Ludwika Bruncra, 1907. XXXI + 296, rys. 82 

— Tom | zesz. 2. Rozwój historyczny pojęć chemicznych. Szkic 
ewolucyi pojęć w mineralogii. Zarys rozwoju matematyki: a) roz- 
wój matematyki do końca XVI w., b) zarys rozwoju geometryi 
w starożytności, wiekach średnich i w epoce odrodzenia, c) rozwój 
matematyki od początku w XVII. W oprac. Leona Marchlewskiego, 
Józefa Siomy, M. Feldbluma, WŁ. "UIAA Stef. Kwietniew- 
skiego, 1911, 279. z 33 ilustracyami 4 

— Тот Il zesz. 1. Historya ogólnej nauki o ziemi (geografii — geolo- 
gii). Dzieje nauk biologicznych. Dzieje antropologii. Dopełnienie 
do historyi fizyki. W oprac. Wacława Natkowskiego, Józefa Nus- 
bauma, Ludwika Krzywickiego i L. Brunera. 1907, 471, 40 illu- 
stracyi w tekście, 2 tablice ' PRI ', 

— Тот Il zesz. 2. Dzieje psychologii. Dzieje językoznawstwa. 
W ae S. Lorii i J. Baudouina de Courtenay. Warszawa, 1909, 
str. AA 

Faraday. Dzicje świecy, przeł, M. i St. Kalinowsey, Warszawa, 1914, 
XXIII + 105, rys 35. . 
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Filipowicz Kazimierz dr. Wiadomości początkowe z botaniki (podług 
dzieła d-ra Le Maout: ,Leçons élémentaires de botanique“) z 194 
drzew. w tekście, 1884, II + 225 -|- lI. (karton.) : 


Grzybowski J. prof. Przeglądowa mapa geologiczna ziem wta 
z tekstem objaśniającym z trzema przekrojami, a redakcya prof. 
J. Morozewicza, wydał Z. Weyberg, 1912, 139, 1 mapa kolorowa . 


Guenther Konrad. Zagadnienia życia w świetle darwinizmu. Z upoważ- 
nienia autora Pal Ad. SCH i Бааша Kulwieć. 1906, 
ХІХ +- 425 . 


Holleman А. F. prof. Podręcznik TO Maa nen sed 
przeł., według 7 wyd. niem. PBA) K. z a M 2. 1910, 
X + 410 0+ I 

Jedrzejewicz 1. Kosmografia. Wyd. 2, оргас. przez d-ra M. Ernsta, 
z 246 fig. w tekście i 11 tabl. 1907, XVI= 442770. : 

Joubert J. Zasady elektryczności Z czwartego wydania PR S 
przełożył Maryan Grotowski. 1915, XV + 507 z 354 rys. w tekście 

Klein P. Metereologia ogólna. Przełożył R. Merecki. Warszawa. 1916, 
99e V ар ЛУИ ыл TES SR ve nr кү; W. ons AE 

рога ејсгук January. Stosunki Póty cane jeziora Switezi. 1915, 

VII + 437. Sprostowania str. 7. ar. CES. An^. 

Kontkiewicz S. Krótki podrecznik mineralogii, 1907, V -|- 226 +- 3 tabl. 
(Karton) : 

Kozłowski Wł. M. Zasady przyrodoznawstwa w świetle teoryi pozna- 
nia, 1905, 311 

Kulwieć Kazimierz. Chrząszcze polskie. Klucz do określania owadów 
tęgopokrywych, dla użytku er amatorów i ОО? 
19072257 

Loth E. Wskazówki do badań antropol na czlowieku żywym 1914 

Malinowski Edmund dr. Świat roślin. O kształtach roślin, powstawaniu 
gatunków, krążeniu soków w roślinach. 1912, VI + 2 nib 145 z 
2 nlb., 108 rys., 2 tabl. barwne. 


Męczkowska T., Rychterówna St. Zbiór cal i AOS AEC Z ps 
rody martwej (202 do$wiadczen z 112 rysunkami), 1915, 156 . 
Mendel G. SS nad mieszańcami roślin, przełoż. W. Wolska. 1915, 
lH + 6 š 

Merczyng H. Teorya prądu elektrycznego. Zarys zasadniczych praw 
ustalonego i nieustalonego prądu elektrycznego i towarzyszących 
mu zakłóceń magnetycznych. Podstawy PR MAŚCI Mod) 
światła. 1905, IX + 92. . . . A 


Merecki R. Klimatologia ziem ED. WEE 1914, 313. DE 
Miłobędzki Tadeusz. Szkoła analizy jakościowej. 1910, VIII—271 (Karton) 
Mohn H. Zasady meteorologii, przełożył St. SE M 1888, XVI + 
218 + VI, 46 rys. i 25 tablic. . | 
Neumayr M. prof. Dzieje ziemi, w opr. prof. d- -ra Wiktora Uhliga: 

I. Geologia ogólna. Wyd. 2 pod red. J. Morozewicza, opra- 
cował K. Koziorowski, z dopełn. M. Limanowskiego. 1912, XX+837, 
mapa barwna, 16 tablic, 300 rys. w tekšcie . 

li. Geologia opisowa, przeł. z 2 niem. wyd. J. Lewiński 
i K. Koziorowski; dopełnienia poczynili: K. Bohdanowicz i J. Grzy- 
bowski. Wydał J. Morozewicz. 1908. XVI + 674 M 343 He 
w tekście, 2 mapy barwne, 9 tabl. (1 kolor.) . 


Nusbaum Józef dr. prof. Zasady anatomii porównawczej. 
. Wiadomości wstępne i anatomia porównawcza zwierząt 
bezkregowych; 1899, III + 744 + XXI, rys 212 + tabl 5. 
ll. Anatomia porównawcza zwierząt kręgowych. 1903, X + 
992, 134 rys. OM T у 
Nusbaum J. dr. Zootomia praktyczna. Wyd. staraniem d-ra Jana Tura, 
zem OUdTzewosytami 641008 eV 111 263,20 KEE ©... ae 
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Nusbaum Hilarowicz Józef. Szlakami nauki ojczystej. Wyd. Jan Tur. 
1916, lll + 280 + 11 potretów. ез РДЫ] A EMS SEM 
Pamietnik Fizyograficzny, wydany staraniem E. Dziewulskiego i B. Zna- 
towicza: 
Tom II. Dział I Meteorologia i hydrografia. IL Geologia 
z chemią. Ill. Botanika i zoologia. IV. Antropologia. V. Miscelanea. 
1883 536 + 2 + 2 + 13 tab., rys lit, 21 drzeworytow w tekście; 
V. Dział I, LI, IIl, IV, V. 1885, 4 nlb. 113 + 76 + 223 + 
74 4- M1 4- 4. 
VIII. Dział I, IT, HI IV, V. 1888, 2 nlb. + XIX + 19 + 
155 + 389 + 17 + 33 + 4 nlb.; 27 tabl rys. lit. i drzew. w tekście; 
Wydawcy: A. Slösarski i Br. Znatowicz. 
IX. Dział I. Il, IM. IV. 1889 2 nib.--XIX-235-45+-11-- 
295 + 77 4-IV, 24 tabl. rys. litogr. i drzewor. 
X. Dział 1, 11, 10, IV. 1890. 2 nib. + XXI + 202 + 75 + 
437 + 2 nlb. + 20 -- 1I + II. 29 tabl. rys. lit. i drzew. w tekście. 


XI. Dział L IL HL 1891, 8 4-18 + 186 + 162 + 133-|-1I-|- . 


H 14 tabl. rys. litogr. i drzewor. w tekscie. 
XH. Dział I, IT, Ш, IV. 1892. 17 + 214 + 235 + 23 + Il + 
II + 12 tabl. rys. litogr. i drzewor. w tekscie. 
XIII. Dział I, 11, HI. 1895, 19 + 152 4-231 + 1-- 1 +7 ta- 
blic rys. lit. 
XIV. Dział L II, Ш. 1896, 23 + 151 + 30 + 228 +1 + I+ 
7 tabl. rys. lit. 
Wydawcy: W. Wróblewski i Br Znatowicz. 
XV. Dział I, II, HI. 1898, 19 + 183 + 285 + 39 +1 + 1 + 
4 mapy + 3 tabl. lit. 
XVI. Dział I, II, Ш. 1900, 13+139+13 + 44-208. 
XVII. Dział 1, II, Ш, IV. 1902, 16-+ 134 + 144 + 104 + 
22 + 1 + 1-+ 1 mapa i tabl. litogr. 
XVIII. Dział 1, H, Ш, IV, V. 1904, 61 + 193 + 147 + 
104 + 24 + 42 +1 + I. 
XIX. Dział I,II, Ш, IV. 1907, 79 + 183 + 59 + 82 + 7 + 1 +-1. 
Pamiętnik Fizyograficzny, wyd star. E. Dziewulskiego i B. Znatowicza. 
XX. Meteorologia i Miscelanea 1910, XLI + 203 + 46, 
Wydawcy: K. Kulwieć i K. Stołyhwo. 
XXI. Dział I, Il, Ш, IV, V. 1913, IX + XV + 155 + 30 +- 
25 + 117 + 48 + 41 +4 mapy + 19 rys. + 24 tabl. fot. 
ХХІІ. I, Il, Ш, IV, V. 1914 IX + XV + 155 + 30 +- 25 + 
117 + 48 + 41 + 4 mapy + 19 rys. + 24 tabl. fot. 
XXIII. Dział I, 11, HI, IV, V. 1916, 262 + XXVIII tabl. + 
142-+ 1 tabl. + 246 + 1 mapa + 16 rys. + 16 -- 8, każdy tom 
Pogorzelski W. Badania teoretyczne ilości ciepła, otrzymywanych na 
d ziemskiej, z uwzględnieniem strat promieniowania w atmos- 
Be 41.3217, TE Wagę ET: 
Pol G. Słownik łacińsko - polski nazw gatunkowych roślin (12 + 17), 
1904, 59 M* rawa Solec TA RT вун NI 
Pawłowski Stanisław. Ze studyów nad zlodowaceniem Czarnohory. 
15, 61, XI tabl. . So. LA oe ui el yb ty v IE 
Pożaryski Mieczysław. Podstawy naukowe elektrotechniki łącznie z za- 
sadami pomiarów. 1915 X + 415, z 427 rysunkami w tekście . 
Roszkowski W. i Żebrowska A. O budowie pochewek prącia u błotnia- 
rekMALOTS POSZLI "tabla UAT PA pa LCR 
Routh E. I. Statyka teoretyczna z licznymi przykładami z drugiego wy- 
dania angielskiego przełożył Zygmunt Straszewicz. 1016.3 5 
458, Hy SEITE) PIETER RAR Kal M Nous 
Rydzewski Bronistaw. Pröba charakterystyki paleobotanicznej Dabrow- 
skiego Zagłębia Węglowego. 1915, 86, tabl. 5 . PASA ARA -- 
Siemiradzki I. Gąbczaki jurajskie ziem polskich (Paleontologia ziem 
polskich pod red. J. Lewińskiego N2 1), 1913, 49 + tabl. VIII 
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Silberstein Ludwik. Elektryczność i magnetyzm. |. 1908, VIII + 366 . 
I. 1910, 304 . . 
JM... cz 15319139 372. . . 
Słownik Geograficzny Król. Pol. i innych krajów słowiańskich. Kompiet 
Słupski Światopełk Zygmunt. Atlas ziem polskich Tom I, część I Wiel- 
kie Księstwo Poznańskie, 46 map i planów. 1915 . . 
Strasburger E. dr., lost L. dr., Schenk K. dr., Karsten 6. dr. Podręcznik 
botaniki dla szkół wyższych. Z XI wyd. niem. przełożyli Jadwiga 
i Karol Szteinbokowie. Zeszyt I. 1913, 160. SECH ]. 1914, 
161 — 320. Zeszyt lll: 1915. 
Świat i człowiek. Zeszyt I, wyd. 2. Pojęcie rozwoju. Wszechświat 
i jego rozwój. Rozwój ziemi, oprac. 1. Waserberg, S. Kramsztyk, 
W. Nałkowski, 1908, XVI + 215 + 82 ilustr. + 3 tablice kolor. . 
Zeszyt II, wyd. 2. Rozwój życia organicznego. Genealogia 
roślin. Genealogia zwierząt. Pochodzenie człowieka. Rozwój czło- 
wieka, opr. J. Nusbaum, Z. Wóycicki, J. Eismont, K. alot wo, 
L. Krzywicki, 1912, 321 -- 73 ilustr. + 1 tabl. 
Zeszyt III, wyd. 2. Rozwój kultury. Rozwój mowy. Roz- 
wój stosunków gospodarczych. W opr. L. Krzywickiego i K. Appe- 
Ja. Warsz. 1912, str. 356 + 65 ilustr. . 
Zeszyt IV, wyd. 2. Rozwój społeczny. Rozwój psychiczny. 
Rozwój w dziejach sztuki. Znaczenie rozwoju. W opr. L. Krzy- 
wickiego, M. Borowskiego, Wł. Tatarkiewicza i F. a wi? 
Warszawa, 1913, str. 355 +5 ilustr. 
Szokalski W. T. Początek i rozwój umysłowości w przyrodzie. ` 1885, 
VIII + 468 +. 
Tenenbaum Szymon. Fauna koleopterologiczna wysp Balearskich. 1915. 
150 + IV . 
Tombeck D. i Gouard E. Chemia przemysłowa, przełożył J. Haraba- 
szewski. 1915. XI + 422. 1 
Tur Jan. Nowe badania nad rozwojem układu nerwowego potworów 
platyneurycznych. 1915, 128. 
Warming E. Zbiorowiska roślinne. Zarys ekologicznej geografii roślin. 
Z wydania niem. E. Knoblaucha przeł. z upow. autora E. Strumpf 
i Trzebiński. 1900, XV + 450 . . 
Witkowski Aug. prof. uniwers. Jagiellońskiego. Zasady fizyki. Tom l, 
wyd. 8. (Fizyka ogólna. Dynamiczne wlasnošci materyi. Akustyka). 
1908, XV + 536 + 205 fig. ӯ 
Tom I, wydanie 4-е, 1915. XX z 535, rys. 205 . 
Tom Il, wyd. 2. (Ciepło. Fizyka cząsteczkowa. Promienio- 
wanie.). 1908, X + 651+ 285 fig. + 2 tabl. kol. . 
Tom Ill. (Elektryczność i NEO 1914, IX + 1 nbl. + 
655 + 326 fig. : 
№. К. Rzeki i jeziora, tekst objaśniający do mapy hydrograf. dawnej 
słowiańszczyzny, część północno - zachodnia. 1883, II+ 125 + 1 nbl. 
Wóycicki Zygmunt. Obrazy roślinności Królestwa Polskiego. Zeszyt 1. 
Roślinność niziny RS AA 1911, 12 nib + tabl. 10 + 20 
str. nlb. objaśnień. . 
Zeszyt ll. Roślinność wyżyny Kielecko - Sandomierskiej. 
1912, 36 + 10 tabl. . . 
Zeszyt III. Roślinność wyżyny Kielecko - Sandomierskiej. 
1912, 32 + 10 tabl. . . 
Zeszyt IV. Roślinność Bolesławiai Olkusza. 1913.34--10tabl. 
Zeszyt V. Roślinność Ojcowa. 1913, 39 + 10 tabl.. 
Zeszyt VI. Roślinność Ojcowa. 1913, 26 + 10 tabl. . 
Zeszyt VII. Roślinność okolic Częstochowy i Olsztyna. 1914. 
30 + 10 tabl. . . : 
,Ziemia" Tygodnik ilustrowany. 'Redaktor i wyd. K. Kulwieé. Rocznie 
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Redaktor i Wydawca 


Jan Tur. 


Adres Redakcyi: Kaliksta N 8 (w lokalu Towarzystwa Nau- 
kowego Warszawskiego). 


Cena kop. 50. 


